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The first line of this example is an initia! line. The second and third lines are continuation 
lines. The statement label is 200, The statement is equivalent to: 


CARD COLUMN 


NOTE 


Where brackets appear in FORTRAN statement formats or other 
formats, they indicate that the data enclosed is optional. 


1.2.3 Comments 


Comments are a convenient method for the programmer to describe the program, A 
comment line cannot be followed by a continuation line. It can only be followed by another 
comment line or by the initial line of a statement. 


Comment lines have the character C or * in column 1. Columns 2 through 72 may be used 
in any desired format. A line with blanks in columns 1 through 72 is also treated as a comment 
line. 


1.2.4 Source Statement Ordering 


Figure 1-1 shows the order in which source statements of each program must be written. 
Within each group the statements may be written in any sequence. A DATA statement must 
appear after any declaration (COMMON, DIMENSION, or type) affecting the variables to be 
initialized. An INTRINSIC, EXTERNAL, or RECURSIVE statement containing function names 
in its declaration list must not appear before any type statements for those function names, 
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PROGRAM, FUNCTION, SUBROUTINE, OR BLOCK DATA STATEMENT 
IMPLICIT STATEMENTS 
OTHER SPECIFICATION STATEMENTS 


EXTERNAL, INTRINSIC, AND 
RECURSIVE STATEMENTS 


Figure 1-1. Ordering of FORTRAN Source Statements 


1.2.5 Program Size Restrictions 


The programmer must observe certain program size restrictions when developing a FORTRAN 
program (due to the nature of the compiler constructions). The number of each of the 
following items within each main program, Subroutine, or function must be less than 1024. 


COMMENT 


FORMAT COPY 


STATEMENTS 


STATEMENTS 


° Scalar variables 

° Array variables 

° Common variables 

° Equivalenced variable names 

° Statement numbers 

° Names in explicit type statements 


° Unique EXTENDED INTEGER, REAL, DOUBLE PRECISION, and COMPLEX constants 


° Unique INTEGER constants 
° Unique LOGICAL constants 


© Unique subprograms called 
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° Arithmetic statement function definitions 


Characters used in CHARACTER constants. If a CHARACTER constant is used 
more than once, only one copy of the constant requires storage allocation. 


1.3 Copy Statement 


The COPY statement enables parts of the program to be stored in more than one file 


and allows many FORTRAN programs to use a single source. This statement has the 
following format. 


FORMAT 
COPY name 


In the preceding format, name is a fully-qualified file access name or synonym (under 
TAXO). The name must be completely contained on one line (not extended across 
continuation lines) and may not contain blanks. A TAXO program may contain any number 
of COPY statements, but they should not be nested deeper than five. The file is rewound 
before the COPY operation begins. 


NOTE 
COPY statements may NOT be nested at all under any TAXO operating system. 


The contents of the named file are inserted into the source program so that the first record 
of the file will be the next line after the COPY statement. Thus, COPY statements may be 
labeled and referenced the same as CONTINUE statements, 


In the following example, the TAXO file. COPY.DATA contains FORTRAN source. This source 
can be included in another FORTRAN source program with a COPY statement as shown. 


EXAMPLE 


COPY .COPY.DATA 


END 


1.4 Fortran Character Set 


The FORTRAN character set consists of 26 letters, 10 digits, and 12 special characters. 


—ee 
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1.4.1 Letters 
A letter is one of the 26 characters: 


ABCDEFGHIJKLMNOPQRSTUVWXYZ 


1.4.2 Digits 


A digit is one of the 10 characters: 


0123456789 


1.4.3 Alphanumeric Characters 


An alphanumeric character is a letter or a digit. 


14h Special Characters 


A special character is one of the 12 characters: 
Character Name of Character 


= Equal 

+ Plus 

- Minus 
Blank 

¥ Asterisk 

/ Slash 

( Left Parenthesis 

) Right Parenthesis 

’ Comma 

. Period 

: Apostrophe 

> Greater Than 


Other characters may appear in a FORTRAN statement only as part of a CHARACTER 
constant or within a FORMAT statement. Any printable character may appear in a comment. 
Refer to Appendix E for additional characters. 


1.5 Programming Technique 


The FORTRAN language was the first of the higher-level languages. No formalization of 
computer languages existed when the first FORTRAN compiler was implemented in 1954; 
no concept of "structured" programming had evolved; and the FORTRAN syntax reflected 
(and still reflects) the architecture of the original target computer. 


Some of the more effective means of structuring programs can be achieved by imposing 
standards of coding more strict than those of the compiler. A well-constructed program 
debugs faster than a poorly-constructed one, an aspect of programming that is rapidly 
becoming more important with rising software costs. 
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This section includes some recommended programming techniques. 


° 


Make liberal use of comments, It is important for the programmer to document 
his program thoroughly so that a user less familiar with it can read through and 
follow it easily. Even an experienced programmer will be able to read through a 
program more easily if it includes comments. One programmer's coding steps 
may not be clear to other programmers. 


- FORTRAN is not one of the more readable languages and is frequently used 
in highly technical, difficult program modules. Comments help the writer of the 


program to think through his problem more clearly and speed up the debugging 
process during initial development. 


- Comments can be read more easily if they are clearly delineated from program 
code, for example by a standard indentation or a border of asterisks, 


Plan statement number assignments carefully. A coherent scheme of statement 
number assignments makes the program steps easier to follow. Do not number 
successive branch targets 1, 2, 3, .... Generally, during the normal coding and \y 


debugging process, a need for more branch points will arise after the first draft 
of the module. 


- Allow increments of at least 10 number assignments. Statements may then be 
inserted later and be assigned numbers that fit into the existing sequence. If 
statement numbers are related to their location in the listing, it is easier for 
other programmers to locate the three branch points of an IF statement or locate 
any other statement that is not executed in program sequence. 


- Another possible convention is to use different digit positions in the statement 
number to indicate the function of the statement, as in the following scheme. 


Statement Number Meaning 

000n0 Branch target 

00n00 DO loop end 

0n000 FORMAT statement 
nooo00d READ/WRITE statement 


- An intelligently planned scheme of statement numbering can be valuable as 
a documentation aid. 


Avoid branching into and out of DO loops. Branches of this type are rarely 
necessary; a program with many of them can be difficult to read and modify. If 
the programmer avoids them, his programs will be more easily understood. 
Whenever feasible, use subroutine CALL statements and function references in 
place of GO TO statements. These substitutions should be made for the following 
reasons. 

- The program will be more readable. 


~- The chances of endless loops in the program are reduced substantially. 


—— SST 
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- This practice encourages the programmer to approach the problem as a 
set of smaller tasks instead of one large program module. 


- Debugging time will be reduced. 


© The programmer may use certain features of FORTRAN-78 as documentation aids, 
if he is willing to sacrifice compatibility with some other compilers. The two features 


are as follows. 


The free format 'F' compiler option may be invoked; therefore, DO loop nests, 
IF blocks, and other logical groupings may be indented to separate them visually 
and improve readability. 


- Identifiers are not limited to six characters in length. This allows variable names 
to be meaningful and readable when a longer word may be spelled out. However, 
because the compiler only looks at the first six characters of the label, duplica- 
tions in these characters must be avoided, 


© Leave as much nonrepetitive program code as possible outside the range of a DO. 
By excluding initialization or other setup code from a DO, execution time may be 
reduced substantially. 


° Use the conditional compile option and D comment records to print program status 
information as a debug aid. 


1.6 Conflicts between FORTRAN-66 and FORTRAN-78 


In April of 1978 the American National Standards Institute released a standard for the FORTRAN 
language that replaced the previous standard released in 1966. The new standard, known 

as X3.9-1978, provides many extensions to the previous FORTRAN standard, which makes 
FORTRAN a more valuable general purpose programming language. FORTRAN-78 has been 
changed to incorporate the entire ANSI X3.9-1978 subset specification as well as several 
features from the full language. These changes cause some conflicts with programs written 

for FORTRAN versions prior to FORTRAN-78. The following describes the conflicts which exist. 


° The concept of Hollerith data has been replaced by the character data type. 
Character data can no longer be stored in program variables unless they are 
explicitly declared to be of type CHARACTER, This restriction applies to assignment 
statements, DATA statements, I/O transmission under the A format specification, 
internal READs and WRITEs, and arrays being used as format references. 


° The way in which DO statements are executed has been changed, If the test for 
completion of the DO loop is true initially, the body of the DO loop is not executed, 
Modification of the loop control variable or any of the loop parameters will not have 
an effect on the number of iterations executed. 


° The evaluation of exponents in an expression has changed; evaluation now goes 
from right to left rather than left to right. 


| 
| 
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2 DATA SPECIFICATIONS 


2.1 General 


This section discusses the handling of data in FORTRAN, an specifies the data conventions 
that are or are not allowed in a program. The topics discussed include: 


° Representation of data in FORTRAN statements 

° Types of data 

° FORTRAN statements that specify or manipulate data 
° Mathematical functions 


The types of data used with FORTRAN-78 are integer, real, double-precision, complex, 
logical, character, and fixed. Data of these classes may be handled as individual data 
elements or grouped in multidimensional arrays. FORTRAN statements define the data cha- 
racteristics and data arrangements for each program. If declarations are not stated, the 
compiler defaults to predefined characteristics, 


2.2 Definitions of Terms 


The terms identifier, variable, and constant, which are used to specify and define data, must 
be understood before any data characteristics are discussed, 


2.2.1 Identifiers 


An identifier is a symbolic name that represents a variable, function, or subprogram ina 
program, The identifier may be any alphanumeric combination as long as the first character 

is a letter. The compiler does not restrict the number of characters that may be used in an 
identifier; however, when the identifier is used in the program, the compiler recognizes only 

the first six characters of the identifier. Therefore, no two identifiers can contain the same 
characters in the first six character positions. So if MODIFICATION is used as an identifier, 
MODIFIER2 cannot be used. The first six characters in both are MODIFI. Blanks are not 
significant within identifiers, therefore, ROOT OF EQUATION in the following is a valid identifier. 


ROOT OF EQUATION = SQRT (FIRST COEFFICIENT) 


2.2.2 Variables 


Variables are data whose values can change during the execution of a program. Because 
the value is dynamic, variables are assigned unique, fixed identifiers for reference within 
the program. In a program, the value of the variable is the most recent value stored in that 


variable's storage area, For example, the following three statements define the relationship 
of the variable identifiers, A and B. 
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EXAMPLES 
A = 8.0 
B=A+2.0 


A=A+ 1.0 


The first statement above assigns a value of 8.0 to A. The second statement defines the 
relationship of B to A. The value of A is used to calculate the value of B. The third statement 
assigns a new value to A by adding 1.0 to the old value and storing the new value in the 
same storage location. A variable identifier may refer to data in any of the seven categories 
of data listed in paragraph 2.1, 


2.2.3 Constants 


Constants are data that do not change in value during the execution of a program. They 
are not assigned identifiers and are referred to by their explicit values. For example, the 
integer seven is represented in a program as the number 7. 


2.3 Type Statements 


The type statement is one example of the FORTRAN statement known as a declaration. 
FORTRAN declarations are used to supply descriptive information about the program rather 
than to specify computation or other action. Such descriptive information is concerned 
primarily with the interpretation of source program identifiers and object program storage 
allocation. Data may be one of seven types. Each type may be explicitly declared in a type 
statement. Integers and real numbers may also have their variable names implicitly declared. 


The seven type statements used to explicitly specify identifier types appearing in the program 
are INTEGER, REAL, DOUBLE PRECISION, COMPLEX, LOGICAL, CHARACTER, and FIXED. 
The following rules apply to these type statements. 


© The name of the statement type is followed by identifiers (used for variable names) 


separated by commas. The FIXED statement includes a scale factor in parenthesesyy 


before the identifier. The CHARACTER statement may include a length specification 
following the keyword CHARACTER, as well as a separate length specification 
following each identifier. 


° Type declaration statements must precede any nondeciarative statements, 
arithmetic function definition statements, or DATA statements in the program. 
Arithmetic function definition statements and DATA statements are discussed in 
Section 6. 


° Each type statement must appear in the program before the first use of any 
identifiers named in the statement. 


° An identifier may appear in only one type statement. 


° Type statements may be used to declare arrays that are not dimensioned in 
either a DIMENSION or a COMMON statement. 


a 
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Integers and real numbers may be declared implicitly using the following rules. 


© If the variable name starts with the letter I, J, K, L, M, or N, it is an 
integer variable. 


© If the variable name starts with any other letter, it is a real variable. 
The rules for integers and real numbers may be altered through the use of an IMPLICIT 
statement. If the programmer wishes to use a variable name for integer or real number data 
that does not conform to the rules above, Such as a rea! variable name beginning with the 
letter |, he may do so. Assignment of variable names to data types in this manner is 
accomplished through the IMPLICIT statement. 
A statement similar to a type statement is available for declaring that a subprogram name 
is an argument in a subprogram call. This statement, the EXTERNAL statement, is described 
in Section 6. The INTRINSIC statement declares that an intrinsic function name is an argument 
in a subprogram call. The INTRINSIC statement is also described in Section 6. 
2.3.1 INTEGER Statement 
FORMATS 

INTEGER identifier, identifier,... 

or 
INTEGER*2 identifier, identifier,... 
and 

INTEGER*4 identifier, identifier,... 
The first two forms are equivalent and are used to declare that the listed identifiers are 
integer type with each datum occupying two bytes, The third form, INTEGER*4, declares 
extended integer type data occupying four bytes each, 
EXAMPLES 

INTEGER ALPHA,P 


INTEGER*4 LAMBDA 
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2.3.2 REAL Statement 
FORMATS 
REAL identifier, identifier,... 
or 
REAL*4 identifier, identifier,... 
and 
REAL*8 identifier, identifier,... 
The first two forms are equivalent and are used to declare that the listed identifiers are real 


type with each datum occupying four bytes in floating-point format (refer to paragraph 2.4.2). 


The last form, REAL*8, declares eight bytes of floating-point data, REAL*8 is equivalent 
to DOUBLE PRECISION. 


EXAMPLES 
REAL LOGX, MASS(10,4) 


REAL*8 A 


2.3.3 DOUBLE PRECISION Statement 
FORMAT 
DOUBLE PRECISION identifier, identifier,... 


This statement declares that the listed identifiers are double-precision type. Double-precision 
data occupies eight bytes in floating-point format (refer to 2.4.3). 


EXAMPLE 


DOUBLE PRECISION RATE,Y,FLOW 


2.3.4 COMPLEX Statement 
FORMAT 


COMPLEX identifier, identifier,... 
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This statement declares that the identifiers are complex type. Each datum occupies eight 
bytes, two floating-point numbers representing the real and imaginary parts. 


EXAMPLE 

COMPLEX ZETA,W,ROOT 
2.3.5 LOGICAL Statement 
FORMAT 


LOGICAL identifier, identifier,... 


This statement declares that the listed identifiers are logical type. Each datum occupies two 
bytes where zero represents .FALSE. and nonzero represents .TRUE.. 


EXAMPLE 


LOGICAL BOOL,P,Q,ANSWER 


2.3.6 CHARACTER Statement 
FORMAT 

CHARACTER #*len identifier *len, identifier *len,... 
Len is an unsigned integer constant in the range 1 through 255 that defines the number of 
characters associated with an identifier. If a length value is given following the CHARACTER 
keyword, it specifies a default length for each identifier that does not have its own length 
specifier. If no length follows the CHARACTER keyword, all identifiers without their own length 


specifier are assumed to be of length one. CHARACTER variables occupy one byte of storage 
for each character of their length. 


EXAMPLE 
CHARACTER *10 NAME,BUFFER*80,COLUMN* 1 


This statement declares NAME to be a character variable having a length of 10, BUFFER 
has a length of 80, and COLUMN has a length of 1. 
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2.3.7 FIXED Statement 


FORMAT 
FIXED (scale) identifier, identifier,... 
Scale is a signed or unsigned integer constant within the range of -31S scale $31. 
The FIXED statement declares that the listed identifiers are fixed type with the stated binary 
scale factor. Each datum occupies two bytes. 
EXAMPLE 
FIXED(3) WEEKLY,NET 
This statement declares that the two variables WEEKLY and NET are fixed type with three 
bits to the left of the binary point (the value is multiplied by 27>). 
EXAMPLE 
FIXED(-2 1)HUGE(10) 
This statement declares an array with the binary point of each element assumed to be 21 
bits to the right of the word; i.e., the value is multiplied by 221, 


Examples of fixed integers are given in paragraph 2.4.7. 


2.3.8 IMPLICIT Statement 


The IMPLICIT statement establishes the implicit type of an identifier if the programmer wishes 
it to be different from the default type. For the default, identifiers beginning with letters |, J, K, 
L, M, and N are integer, and others are real, 


FORMAT 


IMPLICIT type,(A pArpeee) pt yPOrlA 3A yyooe) pees 


Type represents one of the following: INTEGER, INTEGER*2, REAL, REAL*4, REAL*8, DOUBLE 
PRECISION, COMPLEX, LOGICAL, CHARACTER*ien, or FIXED (scale); and A4,Az,... represent 
single alphabetic characters or a range of characters (in alphabetic sequence) denoted by 

the first and last character of the range separated by a minus sign (e.g., A-D). 


This statement applies only to variables not mentioned in an explicit type statement. Variables 
whose first character is listed following a type identifier assume the type appearing before 
the fist. Unlisted variables assume their defauli types. 
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EXAMPLE 
IMPLICIT INTEGER (A-B,X), CHARACTER*50 (C),LOGICAL (L) 


This statement would cause the following implicit declarations to be in effect. 
© |dentifiers beginning with A, B, |, J, K, M, N, and X are integer. 
° |dentifiers beginning with C are character length 50. 
° Identifiers beginning with L are logical. 


© |dentifiers beginning with D, E, F, G, H, O, P, Q, R, S, T, U, V, W, Y, and Z are real. 


2.4 Types of Data 


The compiler recognizes seven types of data as input for processing, These types and the 
number of words of machine storage they occupy are as follows, 


Data Type Number of Storage Bytes 
Integer 2or4 

Real number 4 
Double-precision number 8 

Complex number 8 

Logical quantity 2 

Character string n (where n is string length) 
Fixed 2 


The following paragraphs describe each of these data types. Data type declarations (described 
in paragraph 2.3) are used to specify a particular data type for use in a program. 


2.4.1 Integers 


Integers are whole numbers whose least significant digit represents the one's position. Integers 
may be positive, negative, or zero. If unsigned, they are assumed to be positive. If not explicitly 
declared, an integer variable name must start with the letter I, J, K, L, M, or N. 


Integers cannot contain a decimal point. The compiler assumes that the decimal point is located 
to the right of the right-most digit position. A comma may not be used in a number (to separate 
the thousands digit from the hundreds digit, for example) Only digits are valid characters 

in an integer. 
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Simple integers occupy two, eight-bit bytes; extended integers occupy four, eight-bit bytes, 
In either case the high-order bit is reserved for the sign. Therefore, the range of integers 
in FORTRAN~78 is -32768 to +32767 (-25 to +2 15-1) for simple integers and -2147483648 
to +2147483647 (-231 to +2311) for extended integers. 
The following numbers are valid integers for input to the compiler. 

3298 

-3298 

1234567890 

0 

0001 


The following numbers are not valid integers for input to the compiler. 


3,287 Commas cannot be used in integer data. 
42. Decimal points cannot be used in integer data. 
3234567890 This number exceeds the upper bound of 2147483647 for positive 


integer values. 


Hexadecimal constants are treated as integers or extended integers within expressions. They 
may be written in two forms: 


>D D2... 


| 
In the second form, n is an unsigned integer constant that indicates the number of hexadecimal 
digits in the hexadecimal constant. Z or > indicates hexadecimal, and D represents a hexadecimal 
digit (0-9, A-F). Hexadecimal constants are stored internally, right-justified in one word if there 
are four or less hexadecimal digits, or in two words. 
ww 


The following are examples of valid hexadecimal constants. 
>1B 


| 

| 

| 

>2A | 
3Z1AB | 
| 

| 

| 

| 
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2.4.2 Real Numbers 


Real numbers contain up to seven significant digits and must have a decimal point, Real 
numbers may have an Msi in the input and can represent any values within the approx- 
imate range 10-” to 10 » including zero. A real Number may be represented in either of 


two forms: 


© Digits and decimal point, which may be preceded by a plus or minus sign. 

° Power-of-ten form =— a number times 10 raised to an integer power. 
The power-of-ten form is as follows. 

zzzEn 
The letters zzz represent a numerical value of any number of digits (only the first seven are 
significant) that may or may not contain a decimal point. The letter E indicates that the number 
following it, represented by n, is a power of 10. The power must be an integer constant: 
positive, negative, or zero. The number must be within the specified range for real numbers, 
An unsigned number is assumed to be positive. Commas cannot appear in the input data, 
If not explicitly declared, a real variable name must start with a letter of the alphabet other 


than I, J, K, L, M, or N. 


The following numbers are correct representations of real numbers. 


1387. 
1.23456 

150000. 

1.37E2 (1.37 X 102) 

137E-2 (137 X 1072 or 1.37) 


-4.53716E2  (-4.53716 X 102) 
3E3 (3 X 103) 


The following numbers are not correctly formulated real numbers for input to the compiler, 


387 Numbers not in exponential form require a decimal point. 

3.E The exponent value has been omitted after the letter E. 

3,871.E2 Commas cannot be used in real number data input. 

5E2.5 The exponent is noninteger. 

4.72E76 The number is larger than the upper bound of the permissible range. 


2/10 


The real numbers described in the preceding paragraphs are single-precision. Single-precision 
real numbers are stored in memory in four, eight-bit bytes as illustrated in figure 2-1. Before 
being stored in memory, however, the number is transformed to a normalized hexadecimal 
fraction, a corresponding hexadecimal exponent that is a power of 16 and a sign bit. 


The fraction portion of the number is normalized; that is, it is shifted to the left to eliminate 

leading zeros between the radix point and the first significant bit of the fraction. Normalization 

is by hexadecimal digits; therefore, the number is shifted 4 bits at a time until the left-most 

4 bits are not all zeros, Each bit position shift in the normalization process produces a correspond- 
ing change in the exponent portion of the number to maintain the correct magnitude of the 
number. When completely normalized, the hexadecimal fraction is stored in the second, third, 

and fourth bytes of the four-byte field. The radix point for the fraction is assumed to be 

positioned immediately to the left of the mantissa (at the start of the hexadecimal fraction). 


RADIX POINT 


Oo 1 7 pe 15 
MOST SIGNIFICANT 8 BITS 
ills 5 : pasaprkimeian OF NUMBER NORMALIZED 
HE XADECIMAL 
FRACTION 
2ND WORD 


LEAST SIGNIFICANT 16 BITS OF NUMBER © 


Figure 2-1. Memory Representation of Real Numbers 


The exponent portion of the number is biased by 40 4¢ (excess 64 notation), so that an exponent 
for the number 16° is represented in memory by 404¢. Positive exponents are therefore 
represented by numbers greater than 404g, and negative exponents are represented by numbers 
less than 404¢. For example, 1678 is represented in the exponent field by a value of 384¢- 

The exponent may be any value from 004¢ to 7F 4g, Using the 404¢ bias value, these numbers 
represent exponent values from -404¢ to +3F 4¢ (164 to 1663), The seven exponent bits 

are stored in bits 1 through 7 of the first byte of the four—byte field. 


_— 
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Bit O of the first byte is used as a sign bit. When this bit is zero, the number is positive; when 
this bit is one, the number is negative. 


EXAMPLE 
Hexadecimal Contents of Memory 
Base Ten Number Byte 1 Byte 2 Byte 3 Byte 4 
pins emecnn an renrrmnins nettee C 
1.0 41 10 00 00 
0.5 _ 40 80 00 00 
312500.5 45 4C 4B 48 
100.0 42. 64 00 00 
03125 3F 80 00 00 
-1.0 C1 10 00 00 


2.4.3 Double-Precision Numbers 


Double-precision numbers are similar to real numbers, except that they occupy four more 
memory bytes and provide up to 16 significant digits (56 bits) instead of the seven available 
with real numbers. pouble-precision numbers, like real numbers, are restricted to the range 
of values from 107 to 10 , including zero. The numbers may be positive, negative, or Zero; 
if unsigned, they are assumed to be positive. The numbers must be written in an exponential 
form similar to real numbers. 


zzzDn 
In the preceding form, zzz represents any number of digits (only the first 16 are significant) 
that may or may not include a decimal point. The letter D indicates that the number following 
it is a power of ten; it must be included in the data format. The power, represented by n, 
may be positive, negative, or zero, and must be within the specified range for double-precision 
numbers. Commas cannot appear in the input data. 
The following numbers are correct representations of double-precision numbers. 

3D-16 (3 X 107%) 

487.318264D0 (487.3 18264) 


The following numbers are not correctly formulated double-precision numbers for input to 
the compiler. 


3.8D + 327 The exponent exceeds the upper bound of 75, 
4,134,904D5 Commas cannot be used for double-precision data input 
8.333333D-2.7 A nonintegral exponent of 10 is not allowed. 
8.33333333333 The letter D is not included, 
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Double-precision real numbers are stored in memory in eight, eight-bit bytes as illustrated 

in figure 2-2. The most significant bit of the first byte is a sign bit for the data value field: 

0 if the number is positive and 1 if it is negative. Bits 1 through 7 of the first byte are the 
exponent. The exponent follows the same form as for real number exponents. The remaining 
seven bytes contain the double-precision mantissa with the radix point positioned just to 
the left of the mantissa. 


Oo 1 7 6 1s 
oom [lamer RE 
aTH WORD 


16 LEAST SIGNIFICANT BITS OF NUMBER 


Figure 2-2, Memory Representation of Double-Precision Numbers 


2.4.4 Complex Numbers 


A complex number has two components, a real part and an imaginary part. The imaginary 
part is a real number multiplied byi= -1. 


In FORTRAN-78, complex numbers are written as an ordered pair of real numbers in this form: 


(c C2) 


The variable c, represents the real part of the complex number and co represents the coefficient 
of i for the imaginary part of the complex number. Since c; and cg are both real numbers, 

they must conform to the format for real numbers described in paragraph 2.4.2. Each of these 
numbers may be signed or unsigned. The parentheses and comma are required in the complex 
number representation. The following are correctly formatted complex numbers for use with 

the compiler. 


~ 
(1.25,-8.6) Represents the expression 1.25 - 8.6i 
(0.,0.) Represents the expression 0 + Oi 
(2.5E-7,-3.) Represents the expression 2.5 X 10-7 - 3i 


2.4.5 Character Constants 


A character constant consists of a string of printable characters enclosed in single quotes. 
The maximum number of characters allowed is 255, In order to denote a single quote within 
a character constant, the user needs to have two successive single quotes for every single 
quote that is to appear. Character constants require one byte of storage for each character 
in their length (two successive single quotes count as one character). 
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EXAMPLE 
"LENGTH SHOULDN"T EXCEED 255' 


2.4.6 Logical Constants 


A logical constant may assume one of two logical states, true or false, These states are 
expressed as follows for use with the compiler. 


»TRUE. 


-FALSE, 


The enclosing periods are part of the constant and must be used regardless of context. 

These values can be assigned to variables just as numerical values can, provided the variables 
have been previously declared to be logical variables (refer to the declarations of data types 
in paragraph 2.3). A logical constant occupies two bytes (16 bits) in memory. That word is 
either all zeros (.FALSE.) or any nonzero value (.TRUE.). The following expression is valid 

if the variable A has been declared to be a logical variable. 


A =.TRUE. 


The use of the relational operators and logical operators in logical expressions is described 
in Section 3. 


2.4.7 Fixed-Point Integers 


A fixed-point integer is an integer data type having a binary scale factor associated with 
it. This scaled integer feature provides a mechanism for performing simple floating-point- 
like operations using integer arithmetic, which is much faster than floating-point arithmetic. 
The feature is also useful in scaling and normalization for accuracy optimization. 


Fixed-point integers are represented in source text by variable names declared in a FIXED 
type statement or by constants of the form: 


tiQts 


In the preceding form, i is the integer portion which consists of a string of digits in the range 
~32768 to 32767 with no periods allowed and s is the scale factor portion which consists 
of one or two digits in the range -31 to 31. 


The integer portion is represented internally as a 16-bit, two's complement binary number, 
The scale factor represents the position of the integer binary point. A scale factor of zero 
indicates that the binary point is just right of the least significant bit of the integer portion. 
This orientation of the binary point is identical to that of the usual integer data type. A positive 
scale factor (s) indicates the binary point is s bits to the left of the least significant bit; a 
negative scale factor, s bits to the right. 
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EXAMPLES 


19Q0 is represented as: 


efefefofofofofofefefe]+fofof || 
si, a 


BINARY POINT 


o 
~ 


19Q3 is represented as: 


oa, ee 


BINARY POINT 


and 19Q-3 is represented as: 


VIRTUAL BIT 
POSITIONS 


poe eh 8) a] | 9) pt sey 
= 


BINARY POINT 


NOTE | 


The binary point is not stored internally as a part of the integer portion. Its position 
is maintained separately by the FORTRAN compiler. 


The following are examples of valid fixed-point integers. 


512Q0 

256Q-1 All of these have a decimal value of 512. 
4096Q+3 

13Q0 is equivalent to a decimal value of 13. 
13Q4 is equivalent to a decimal value of 0.8125. 
13Q-2 is equivalent to a decimal value of 52. 


Refer to Appendix D for a more detailed discussion of fixed-point integer use. 
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2.5 Array Data 


A programmer must often perform calculations using large sets of related data. It is generally 
convenient to handle such data as a unit, and FORTRAN allows a single variable name to 

be assigned to such a set of data. To be used in a program, the set of data consists of a 
specified number of individual data items. The discussion of arrays, subscripts, and the DIMENSION 
statement in the following paragraphs describes the handling of grouped and related data in 
FORTRAN. 


2.5.1 Arrays and Subscripts 


A set of related data is called an array, reflecting the fact that the data can be arranged into 
lists, or tables of rows and columns, Table 2-1 shows a group of data presented in a two- 
dimensional array. Each number in the table represents the average precipitation in inches 

of a particular city during a particular month, and is an element of the array. If the array is 
given the variable name PRECIP, any element in the array can be identified by a row number 
and a column number following the name, for example PRECIP (3,5) = 4.8, the average 
precipitation in Dallas (row) in May (column). The numbers 3 and 5 are called subscripts. 

The number of supscripts following the variable name is the same as the number of dimensions 
in the array. Commas separate the subscripts. 


Any element in the array in table 2-1 could be denoted as PRECIP (I,J), where | equals an 
integer value from 1 to 6 for the cities and J equals an integer value from 1 to 12 for the 
months, Arrays are often used in programs where a computation is to be performed on each 
element in turn, 


An example of a one-dimensional array is the data for the average precipitation in Atlanta 

in each month, If this array is named P'1, it has as its elements P1(1) = 4.4, P1(2) = 4.5, and 
so forth through P1(12) = 4.4, 

The number of subscripts permissible in a FORTRAN-78 array is limited to three. A sufficiently 
large area of storage must be available to accommodate the array. 


Table 2-1. Example of an Array: Precipitation in Selected United States Cities 


(1) (2) (3) 4) BY (9) (10 (ATP (12) 


City Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. 
(1) Atlanta 4.4 4.5 5.4 4.5 3.2 3.8 4.7 3.6 3.3 2.4 3.0 4.4 
(2) Chicago 1.9 16 27 3.0 3.7 4.1 3.4 3.2 rf 2.8 ae 1.9 
(3) Dallas 2.3 2.6 2.9 40 48 3.2 1.9 1.9 2.8 2.7 Re 2.7 


(4) Los Angeles 3.1 3.3 2.3 12 O02 1 00 00 O02 04 +t 2.9 
(5) New York 33 28 40 34 37 33 3.7 £44 3.9 3.1 3.4 3.3 


(6) Seaitle 5.7 42 38 2.4 Bey 4 16 08 1.0 2.1 40 54 63 
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2.5.2 Storage of Array Elements 


Arrays are stored in storage locations in the following manner. The first subscript is incremented 
through its possible values; then the second subscript is incremented by one and the first 
subscript is again incremented through its possible values. When the second subscript has 

been incremented through its possible values in this manner, the third subscript is incremented 
by one and the first and second subscripts are again incremented through their ranges of 
values. 


For example, the 3 by 3 by 3 array variable A(I,J,K) would be stored in this order, 
A(1,1,1);A(2,1,1);A(3, 1, 1);A(1,2, 1);A(2,2, 1);A(3,2, 1);A(1,3, 1); 
A(2,3,1);A(3,3,1);A(1,1,2);A(2,1,2);A(, 1,2);A(1,2,2);A(2,2,2); 


A(3,2,2);A(1,3,2);A(2,3,2);A(3,3,2);A(1, 1,3);A(2, 1,3);A(3,1,3); 
A(1,2,3);A(2,2,3);A(3,2,3);A(1,3,3);A(2,3,3);A(3,3,3) - 


The identifiers for arrays follow the same rules as identifiers for undimensioned variables. 
For example, these are valid array names: 

ARRAY 

E 

ST64A 
These are not valid array names: 


2X2 An identifier may not begin with a numeral. 


RUN# Characters other than letters and numerals may not appear in an identifier. 


2.5.3 DIMENSION Statement 
The DIMENSION statement is used to declare that identifiers are array identifiers and to specify 
the number and bounds of the array subscripts. The information supplied in a DIMENSION state- 
ment is required for the allocation of memory for arrays. Any number of arrays may be decla_) 
in a single DIMENSION statement. 
FORMAT 

DIMENSION S qS 2:00 


S is an array specification. 


Each array specification gives the array identifier and the maximum values that each of its 
subscripts may assume, thus: 


identifier(max , [,maxg [,max3]]) 
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The maxima must be integer constants in a main program or integer constants or variables 

in a subprogram. An array may have one, two, or three dimensions, The dimension information 
for an array must be given before the occurrence of that array variable in the program, 
Subscripts must not exceed their maximum values declared in the DIMENSION statement. 


EXAMPLE 
DIMENSION EDGE(10,8) 


This statement specifies that EDGE is a two-dimensional array. The first subscript may vary 
from 1 to 10 inclusive, and the second from 1 to 8 inclusive, Zero and negative subscripts 
are not permitted. 


EXAMPLE 

DIMENSION PLACE (3,3,3),HI(2,4),K(256) 
Arrays may also be declared in COMMON or type declaration statements in the same way. 
EXAMPLES 

COMMON X(10,4),Y,Z 

INTEGER A(7,32),B 

DOUBLE PRECISION K(6, 10) 


Each array variable appearing in the program must represent an array element declared in 

a DIMENSION statement, unless the dimension information is given in another statement. When 
the dimension information is provided in a COMMON or type declaration statement, it must 
not appear in a DIMENSION statement. 


Note that each element of an integer or logical array requires two bytes of storage. Each 
element of a real array requires four bytes of storage, double-precision elements require eight 
bytes of storage, and complex require eight bytes of storage. Character array elements require 
between 1 and 255 bytes of storage depending on their declared length. 


Within a subprogram, the maximum array subscript values may be integer variables instead 

of constants, provided that the array name and variable dimensions are dummy arguments 
(refer to Section 6) of the subprogram. The actual array name and values of the dummy 
variables are given by the calling program when the subprogram is called, The maximum array 
subscript values are not changed by any subsequent modification of the dimension variables, 
In addition, arrays which are dummy arguments of a subprogram may have the last dimension 
factor specified by an asterisk(*), This type of dummy array is known as an assumed-size 
dummy array. Subscripts for the assumed-~size dimension may take on any positive integer 
value, as long as the total size of the array is not exceeded, 
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EXAMPLE 
DIMENSION BETA(L,M,G),B(20),C(N,*) 


The identifiers BETA, L, M, G, C, and N must all be dummy arguments. 


2.5.4 Array Bounds Checking 


A common source of program errors is the referencing of an array outside its declared bounds, 
This results in a reference to, and possible alteration of, memory adjacent to that array. The 
FORTRAN-~78 compiler has an option which may be selected to ensure that array subscripts 
remain within their declared bounds. The 'B' compiler option causes the compiler to check 

all constant subscripts at compile time and issue syntax errors for illegal references. It also 
causes the compiler to generate object code that checks variable array subscripts at runtime. 
If any subscript calculation falls out of bounds, a fatal error is given and execution terminates. 
The bounds for each subscript in an array reference are 1 < subscript value =n, where n 


is the maximum value declared for that dimension. In the case of assumed-size dummy arrays. 


the last subscript may have any positive nonzero value, 


The BOUNDS OFF and BOUNDS ON statements are provided to control the bounds checking 
feature. The syntax for these statements follows. 


FORMATS 


BOUNDS OFF 


(any statements) 


BOUNDS ON 


When the 'B' option has been selected, these statements are used to selectively turn the 
bounds checking feature off and, if desired, on again. If the 'B' option is not in effect, these 
statements have no effect. Since the bounds checking option creates longer object modules, 
the BOUNDS OFF and BOUNDS ON statements can be used to control the amount of overhead 
incurred. 


2.6 COMMON and EQUIVALENCE Statements 


The COMMON and EQUIVALENCE statements provide flexibility in the naming of variables 
and the assignment of storage locations to them. 


The COMMON statement is used to make an identifier in a program and in separate subprograms 


refer to the same variable, Generally, variables that have the same name but are in different 
programs or subprograms are distinct variables. The COMMON siatement assigns these distinct 
variables to the same storage location. 
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The EQUIVALENCE statement causes two or more variables in a single program or subprogram 
to be placed in the same storage location. This statement can be used for these purposes: 

© Two or more identifiers may be defined to have identical meanings. 

° A storage location can be used to contain two or more variables that are never 

needed at the same time. 

2.6.1 COMMON Statement 
A COMMON statement has the following form. 
FORMAT 

COMMON A,B(2,3),C,D(4) 
This statement causes A to be stored in the first location of the common area, the array B 
to be stored in the next six locations (in the order specified in paragraph 2.5.2), then C in 
the following location, followed by the array D in four locations, 
One COMMON statement by itself does not cause multiple assignments of an area of storage. 
One or more additional COMMON statements, each assigning storage in a specified order 
to variables, accomplishes this. For example, if one program contains the above example 
and a subprogram contains this statement: 


COMMON G(6),F,E(2,2),H 


then the common area shared by the two programs or subporgrams will contain the following 
data. Each rectangle represents a storage location. 


B(1,1) G(2) 
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The preceding illustration is an example of blank common, where all COMMON statements 
refer to a single unnamed common area of storage. It is possible to have a number of distinct 
common areas or blocks, each named with an identifier. A named common block is known 
as labeled common. 


The common area may be divided into separate blocks identified by block names, A block 
is specified in this way: 


COMMON identifier /identifier 4, identifiers,...,identifier,, 
The identifier enclosed in slashes is the block name. The subscripted identifiers that follow 
are the names of the variables or arrays assigned to the block. These elements are placed 


in the block in the order in which they appear in the block specification. 


Quantities placed in a common block by means of equivalences (refer to paragraph 2.6.3) 
may cause the end of the common block to be extended, 


EXAMPLES 
COMMON/R/X,Y,Z 
DIMENSION A(4) 
EQUIVALENCE(A,Y) 


These statements cause the common block R to extend from X to A(4), arranged as follows. 


Equivalence, which causes extension of the start of a common block, is not allowed, For 
example, the sequence 


COMMON/R/X,Y,Z 
DIMENSION A(4) 
EQUIVALENCE(X,A(3)) 


is not permitted since it requires block R to be arranged. 
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A(1) and A(2) extend the start of block R. 
Redundant COMMON entries are not allowed. For example, the following is invalid. 
COMMON A,B,C,A 
A variable written with subscripting information in a COMMON or type statement must not 
be mentioned in a DIMENSION statement. A variable in a DIMENSION statement, however, 
may be mentioned without subscript dimensions in a COMMON or type statement. 
The block list of the COMMON statement consists of a sequence of one or more block 
specifications. The following statement indicates that the elements X,Y, and T, in that order, 
are to be placed in block R and that U, V, W, and Z are to be placed in block C, 
EXAMPLE 
COMMON/R/X,Y,T/C/U,V,W,Z 


Block entries are concatenated throughout the program unit, beginning with the first COMMON 
statement. For example, the statements 


COMMON/D/ALPHA/R/A,B/C/S 
COMMON/C/X,Y/R/U,V,W 


have the same effect as the statement 
COMMON/D/ALPHA/R/A,B,U,V,W/C/S,X,Y 
One block of common storage may be left unlabeled, in effect making it blank common, Blank 
common is indicated by two consecutive slashes. The following statement indicates that B, 
C, and D are placed in blank common. 
COMMON/R/X,Y//B,C,D 
The slashes may be omitted when blank common is the first block of the statement. 
EXAMPLE 
COMMON B,C,D 
Array names appearing in COMMON statements may have dimension information appended, 


as in a DIMENSION statement. The following statement specifies the dimensions of the array T 
while entering T in blank common. 
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EXAMPLE 
COMMON ALPHA,T(15,10,5),GAMMA 
NOTE 


Dummy arguments for SUBROUTINE or FUNCTION statements 
cannot appear in COMMON statements. 


A single COMMON statement may contain variable names, array names, and dimensioned 
array names (but nor array elements). The following are valid, 


EXAMPLE 


DIMENSION B(5,15) 
COMMON A,B,C(9,9,9) 


Variables or arrays that appear in the main program or in a subprogram may be made to share 
the same storage locations with variables or arrays in other subprograms, by use of the COMQyN 
statement. For example, if one program module contains the statement 
COMMON TABLE,A,B,C 
and a second program module contains the statement 
COMMON LIST 
the variable names TABLE and LIST refer to the same storage locations, 
If the main program contains the statement 
COMMON A,B,C 
and a subprogram contains the statement 
COMMON X,Y,Z,XX,YY,ZZ 
and A, B, and C are equal in length to X, Y, and Z, respectively, then A and X refer to the ww 
same storage locations, as'do B and Y, and C and Z. XX, YY, and ZZ are unique; that is, 
they do not share a storage location with any other variables, 
Within a specific program or subprogram, variables and arrays are assigned storage locations 
in the sequence in which their names appear in a COMMON statement. Subsequent sequentia! 


storage assignments within the same program or subprogram are made with additional COMMON 
statements, 


—— ee OOS 
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A dummy variable can be used in a COMMON statement to establish shared locations for 

variables that would otherwise occupy different locations. For example, the variable Z of 

the previous example will share storage with S if the following statement is used: 
COMMON Q,R,S 


Q and R are dummy names that are not used elsewhere in the program. 


2.6.2 SAVE Statement 
The purpose of the SAVE statement is to guarantee that entities in labeled common retain 
their definition status upon completion of a given subprogram, This statement has no actual 
effect, but is included for conformance to the ANSI standard. The format of this statement 
follows. 
FORMAT 

SAVE /name 4/,/name2/,...,/Namep/ 


The names which appear in the list of the SAVE statement must be the names of labeled com- 
mon blocks which have been declared in that program module, 


EXAMPLE 


SAVE/BLOCK 1/,/COM/ 


2.6.3 EQUIVALENCE Statement 
The EQUIVALENCE statement allows more than one identifier to represent the same quantity. 
FORMAT 

EQUIVALENCE (r ql Qpeee gl) (S S2reeeySrp) peony lt qrtapeeesth ) 
The variables r, s, and t represent the equivalenced variables, and n, m, and k represent 
the number of variables equivalenced. Variables r, through r, are equivalent, that is, they 
represent the same quantity. Variables s; through s,, all represent a second quantity. 
The references of an EQUIVALENCE statement may be variables or array identifiers or array 


element references. The subscripts of an array element must be integer constants. The number 
of subscripts must be equal to the array dimension or.must be one, 
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Since entire arrays are shifted to satisfy the equivalence, only the relative positions of the 
references are important. In the following example, both statements accomplish the same 
result. 


EXAMPLES 
EQUIVALENCE (BETA(1),ALPHA(7)) 
EQUIVALENCE (BETA(2),ALPHA(8)) 
Note that the relation of equivalence is transitive, e.g., the two statements 


EQUIVALENCE (A,B),(B,C) 
EQUIVALENCE (A,B,C) | 


have the same effect. Redundant equivalence statements are allowed, 
EXAMPLES | 
EQUIVALENCE (X,Y) 


EQUIVALENCE (Y,X) 


2.6.4 EQUIVALENCE and COMMON Interaction 


Identifiers may appear in both COMMON and EQUIVALENCE statements, provided no two 
quantities in common are set equivalent to one another. 
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3 OPERATORS AND EXPRESSIONS 


3.1 General 


The data used in FORTRAN computation — constants, variables, and functions —= may 
be combined or altered using mathematical operators. The operators indicate the manner 
in which a group of values is to be changed to obtain a new value. 


A sequence of mathematical or logical elements, separated by operators and possibly by 
parentheses according to specific rules, is called an expression, The mathematical elements 
are constants, variables, and functios. The parentheses indicate how elements are grouped 
within the expression. 


FORTRAN expressions belong to two general classes: arithmetic and logical. The mathematical 
data values in the expression and the types of operators that relate them determine the class 
of a given expression. 


The seven data types described in Section 2 may be divided into those which have numeric 
values, those which have character string values, and those which have logical values, 
Integers, real numbers, double-precision numbers, fixed, and complex numbers are all numeric 
data. Character data consists of strings of characters from the FORTRAN character set, Logical 
data is distinct from the others because it may only assume the value true or false. 


Data values may be combined or transformed by means of operators, yielding a result which 
is a new data value. The FORTRAN operators may be any of three types. 


° Arithmetic operators, which specify operations on numeric quantities. 


° Relational operators, which express a logical relation between numeric values or 
between character strings. The result of the operation is a logical value. 


° Logical operators, which specify operations on logical quantities. 


The class of an expression is determined by the results of the operations within it. If the results 
are numeric values, the expression is arithmetic; if they are logical values, the expression 

is logical. Arithmetic expressions contain arithmetic operators, while logical expressions contain 
relational and logical operators. 


3.2 Arithmetic Expressions 


The following paragraphs describe the arithmetic operators and the rules for forming and 
evaluating arithmetic expressions. 
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3.2.1 Arithmetic Operations 
There are five arithmetic operations, represented in FORTRAN by the following symbols, 


Addition 
Subtraction 
Multiplication 
Division 
Exponentiation 


Net + 


Two of these symbols are also used as algebraic signs to indicate positive or negative 
values. 


+ Unary plus 
- Unary minus 


These algebraic signs are handled as operators by the FORTRAN compiler. They are read 
as positive and negative signs rather than addition and subtraction operators if they are not 
immediately preceded by an arithmetic quantity. 


If the precedence of operations is not given explicitly by parentheses, the order of precedence 
is taken to be the following. 


Operator Operation 

ek Exponentiation (highest precedence) 
* and / Multiplication and division 

+ and - Addition, unary plus, subtraction, and 


unary minus (lowest precedence) 


When the operators are of equal precedence, they are understood to be grouped from the 

left, except for exponentiation, which is performed from right to left. Thus, the expression 
A+B+C+#D is calculated as ((A + B) +C) + D. This convention takes care of such ambiguous 
expressions as 


A**B**C 
X/Y/Z 


In FORTRAN-78, these expressions are calculated as 


A**(B**C) 
(X/Y)/Z 
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3.2.2 Formation of Arithmetic Expressions 


An arithmetic expression consists of numerical elements, arithmetic operators, and paren- 
theses. It has a single value which is the result of the calculations specified by the quantities 


and operators in the expression. 
The rules for forming arithmetic expressions are the following. 


° An expression may consist of a single element (constant, variable, or function 
reference). (Refer to Section 6.) 


2.71828 
Z(N) 
TAN(THETA) 


© Compound expressions may be formed by using operators to combine basic 
elements. 


X+3 
TOTAL/POINTS 
TAN(PI*M) 


° Any expression may be enclosed in parentheses and considered to be a 
a single element. 


(X+Y)/2 
(ZETA) 
COS (SIN(PI*M)+X) 

° Expressions may be preceded by a unary plus or unary minus, 
+X 
-(ALPHA*BETA) 
-SQRT (-GAMMA) 

© No two operators may appear in sequence. For instance, 
X*-Y 
is improper. Use of parentheses yields the correct form. 
X*(-Y) 


° The precedence of operations, if not given explicitly with parentheses, 
is according to the rules in paragraph 3.2.1. For example, the expression 


A*B+C/D**E 
is taken to be 


(A*B)+(C/(D**E)) 
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° Sequences of operations are evaluated according ro the rules in paragraphs 3.2.1 


and 3.2.3. 


° Parentheses may not be used to imply multiplication, The asterisk operator must 


always be used for this purpose. 
Therefore, the algebraic expression 


(ab)(c%) 
must be written as 


(A*B) * (C**D) 
Hexadecimal constants may be used in the algebraic expressions, For example 


b= >73C+4 


3.2.3 Evaluation of Arithmetic Expressions 
The value of an arithmetic expression may be of the integer, fixed, real, double-precision, 
or complex type. If more than one of these types appears in an expression, the expression's 
data type is that of the highest rank of any element in the expression. The arithmetic data 
types are ranked as follows. 

° Complex (highest rank) 

° Double-precision 

° Real 

° Fixed 

° Extended integer 

© Integer (lowest rank) 
The type of a subscripted variable, however, is determined by the type of the variable identi. 
Each operation within an expression is evaluated in the type of the operator's highest ranking 
operand, Thus, the evaluation of an expression is not changed to a higher rank until necessary. 
For example, if | and J are integers, R is REAL*4, DP is a REAL*8 number, and C is a complex 
variable, then the expression 

(I/ J+R)*DP*C 
is evaluated as 


CMPLX((DBLE(REAL(I/J)+R)*DP),0)*C 


The variable CMPLX is a function that converts double-precision to complex, DBLE is a function 
that converts real to double-precision, and REAL is a function that converts integer to real. 
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integer and fixed espressions are evaluated using binary integer arithmetic throughout, giving 
an integer value as the result. In integer arithmetic, fracional parts arising in division are 
truncated, not rounded. For example 


7/3 yields 2 
6/7 yields 0 


All other calculations use binary, floating-point arithmetic. 


Conversions to a higher rank are performed as follows. 


° An integer quantity becomes the integer part of a fixed or real quantity. The 
fractional part is zero. 


° A fixed quantity is converted to real so that the part to the left of the binary point 
is the integer part and the part to the right of the binary point is the fractional part. 


° Areal quantity becomes the most significant part of a double-precision real quantity. 
The least significant part is zero. 


° Areal quantity becomes the real part of a complex quantity. The imaginary part 
is zero. 


° A double-precision quantity is converted to single-precision and becomes the real 
part of a complex quantity. The imaginary part is zero. 


In exponentiation (**), the types of the base and exponent are restricted as follows. 
° Complex exponents are not allowed, 
° Only integer exponents may be used with complex bases. 
° Fixed bases are not allowed. 

In exponentiation (**), the evaluation is as follows. 


° If the base is exponentiated by an integer expression, the result is the same type 
as the base. 


° An INTEGER*2 or INTEGER*4 expression may be exponentiated by an INTEGER*2 
or INTEGER*4 expression. The result of two INTEGER*2 expressions is INTEGER*2; 
otherwise, the result is INTEGER*4, 


° Areal or double-precision base may have a real or double-precision exponent; 
only if both parts are real will the result be real, 


° With a complex base and integer exponent, the result is complex. 
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Operations are performed by evaluating exponentiation first, then multiplication and division, 
then addition, subtraction, unary plus and unary minus. Operations of the same priority are 
evaluated left to right, except for exponentiation. The following example illustrates how this 
process is carried out for the five arithmetic operations. 


EXAMPLE 


2**°3/4+5°6-7 
8/4 + 5°6-7 
2+5°6-7 
2+30-7 
32-7 
25 


Expressions within parentheses are evaluated first. Nested groups of parentheses are evaluated 

by starting with the innermost sets of parentheses and working outward. Nesting of parentheses 
to any level is permitted in FORTRAN. The following example illustrates how this type of exp\y- 
sion is evaluated. 


EXAMPLE 


4*(5-((17 + 37)/3**(1 + 2))) + ((21-9)/(10-6) + 13) 
4*(5-(54/3* *(1 + 2))) + ((21-9)/(10-6) + 43) 
4*(5-(54/3* *3)) + ((21-9)/(10-6) + 13) 
4*(5-(54/27)) + ((21-9)/(10-6) + 13) 

4* (5-2) + ((21-9)/(10-6) + 13) 
4°3 + ((21-9)/(10-6) + 13) 

4*3 + (12/(10-6) + 13) 
4*3+(12/4 + 13) 
4*3+(3+ 13) 
4°*3+16 
12+°16 
28 


When several variables or numbers in sequence are to be multiplied and/or divided, they 
are calculated as if grouped from the left. For example, these two are equivalent: 


A*B/C*D 
((A*B)/C)*D 


Similarly, a string of variables or numbers to be added and/or subtracted are calculated as 
if grouped from the left. These two are also equivalent. 


A-B+C-D 
((A-B)+C)-D 


> a ee a ~ = — 
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3.3 Logical Expressions 


Four basic elements are used in FORTRAN logical computations: constants, variables, func- 
tional references, and relations. These basic elements all represent logical quantitites. Logical 
identifiers are written in the same way as numerical identifiers; however, they must always 
be explicitly declared to be of the logical types. 


A logical quantity may have either of two values: true or false. Logical quantities occupy 
one word. 


The following paragraphs discuss relational operations, logical operations, and the formation 
and evaluation of logical expressions. 


3.3.1 Arithmetic Relational Operations 


Arithmetic relations are constructed from arithmetic expressions of the integer, fixed, real 
or double-precision type. These relations make use of the relational operators: 


-GT. Greater than 

GE. Greater than or equal to 

LT. Less than 

-LE. Less than or equal to 

-EQ. Equal to 

-NE. Not equal to 
The enclosing periods are part of the operator and must be present. 
Complex operands are allowed for .EQ. and .NE. only, Relational operators have lower 
precedence than arithmetic operators. Two expressions of integer, fixed, real, or double- 
precision type, separated by a relational operator, form a relation, The following example 
is a relation. 
EXAMPLE 

X+2 .LE. 3*Y 
The entire relation constitutes a basic logical element. The value of such an element is .TRUE. 
if the relation expressed is .TRUE.; otherwise, the value is .FALSE.. In the example above, 


the element has the value . TRUE. if X is 2 and Y is 2, and the value .FALSE. if X is 2 and 
Yis 1. 
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The following are examples of valid relational expressions, 
EXAMPLES 


Q1.GT. Q2 
Q3**3 .EQ. (6*Q4+2) 
0.4*FLN .LE. 1.83D4 


Q1, Q2, Q3, Q4, and FLN are arithmetic quantitites. In the following examples, A is an 
arithmetic variable and S, T1, and T2 are logical variables. The following are examples of 
illegal relational expressions, 


EXAMPLES 
A.GT.(4.345,7.6 14) Complex numbers are illegal in logical expressions except .EQ. 
and .NE.. 
S.GE.(T 1+T2) Logical quantities may not be joined by relational operators. 


w 
3.3.2 Character Relational Operations 


The relational operators may be used with operands of type character to produce lexical 
comparisons between text strings. Comparisons are performed character by character from 

left to right. If the operands are of unequal length, the shorter string is assumed to be extended 
on the right with blanks. The collating sequence used for comparisons is the same as the 
internal ASCII collating sequence for the TA 1600 computer. The following are examples of 
valid character relations. 


EXAMPLES | 
A .EQ. 'HELLO' | 
B .NE. "YES' 


"STRINGA' .LT. 'STRINGB' 
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3.3.3 Logical Operations 


There are three logical operators, which are arranged according to precedence in the evaluation 
of an expression. 


NOT. Logical NOT (highest precedence) 
-AND. Logical AND 
-OR. Logical OR (lowest precedence) 
The enclosing periods are part of the operators and must be present. 
The logical operations have specific meanings: 
© The expression .NOT.P is true if P is false and false if P is true, 
© The expression P.AND.Q is true if P and Q are both true; otherwise, it is false, 
° The expression P.OR.Q is false if P and Q are both false; otherwise, it is true. 


In a FORTRAN program, P and Q would have been explicitly defined as logical variables, 


3.3.4 Formation of Logical Expressions 


Logical expressions are formed according to the following rules, 
° A logical expression may consist of a single logical element. For example: 


-TRUE. 
BOOL(N) 
X.GE.3.14159 


° Single elements may be combined through use of the logical operators AND. and 
-OR. to form compound expressions, such as 


TVAL.AND.INDEX 
BOOL(M).OR.K.EQ.LIMIT 


° Any logical expression may be enclosed in parentheses and regarded as an element. 


(T.OR.S).AND.(R.OR.Q) 
(BOOL(M)) 
PARITY.OR.((2.GT. Y.OR.X. GE. Y). AND.NEVER) 
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° Any logical expression may be preceded by the operator .NOT. as in 


eNOT.T 
eNOT.X+7,.GT.Y+Z 
BOOL(K).AND..NOT.(TVAL.OR.R) 


The only situation in which two logical operators may occur in sequence is when 
the second operator is .NOT.. An example is 


P.AND..NOT.Q 


° Expressions whose values are either true or false may be combined using logical 
operators. 


3.4 Evaluation of Logical Expressions 

When logical expressions are evaluated, the relational operations have precedence over the 
logical operations; in other words, the relational operations are evaluated before the logical \) 
operations. 


Among the logical operators, .NOT. is evaluated first; then .AND.; .OR, is evaluated last. 
Therefore, the expression 


T.AND..NOT.S.OR..NOT.P.AND.R 
is interpreted 


(T.AND.(.NOT.S)).OR.((.NOT.P).AND.R) 


3.5 Summary of Rules vor Expressions 


The following paragraphs summarize the hierarchy of operator precedence and provide some | 
general comments about permissible arithmetic and logical statements. 


The hierarchy of precedence of arithmetic, relational, and logical operators in FORTRAN is 
summarized in table 3-1. The first item in the table has the highest precedence if functions |) 
and expressions within parentheses have been evaluated, Note that arithmetic operators 


are evaluated before relational and logical operators. All operators of the same precedence | 
are evaluated from left to right, except for exponentiation, | 


For example, the logical expression 
-NOT.ZETA**2+Y*MASS.GT.K-2.0R.PARITY.AND.X.EQ.Y 
is interpreted 


(NOT.(((ZETA**2)+(Y *MASS)).GT.(K-2))). OR.(PARITY.AND.(X.EQ. Y)) 
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Table 3-1 Precedence of Operators 


Operator Operation 

ne Exponentiation 

* and / Multiplication and division 

+ and - Addition, positive sign, 
negative sign, and subtrac- 
tion 

.GT. 

GE 

aT. Relational operations 

LE; 

EQ. 

.NE. 

-NOT. - Logical NOT 

-AND. Logical AND 

OR. Logical OR 


3.6 Assignment Statements 


Assignment statements assign values to variables. FORTRAN has two types of assignment 
statements: arithmetic and logical, and character. 


3.6.1 Arithmetic and Logical Assignment Statements 


In arithmetic and logical assignment statements, a variable name appears to the left of an 
equal sign and an expression appears to the right. The value of the variable is defined by 
evaluating the expression. These assignment statements have the following form. 


FORMAT 
variable = expression 
Expression may be any arithmetic or logical expression. 


The statement defines the value of the variable by evaluating the expression to the right 
of the equal sign. 


Note that an assignment statement is not a mathematical equation, for the "=" sign is not 
used in its usual mathematical sense. In a FORTRAN statement, the "=" means "is replaced 
by the value of" rather than "is equal to the value of." For example, the statement 


N=N+1 


is interpreted to mean "N is replaced by the value N + 1", i.e., one is added to the value 
of N. N is assigned a new value which is greater than its previous value by 1. 


EXAMPLES 


Y = 2*Y 


P = .TRUE, 


3/12 


X(N) = NXZETA(ALPHA*M/PI) 


Type of Variable 


(Left Side of 
Statement) 


INTEGER*2 
INTEGER*2 


INTEGER*2 
INTEGER*2 
INTEGER*2 


INTEGER*2 


INTEGER*4 


INTEGER*4 


INTEGER*4 


INTEGER*4 


Assignment Statement Type Conversion 


Type of Expression 
(Right Side of 


Statement) 


INTEGER*2 
INTEGER*4 


FIXED 


REAL 


DOUBLE PRECISION 


COMPLEX 


INTEGER*2 


INTEGER*4 


FIXED 


REAL 


If the type of the expression on the right of the assignment statement differs from the type 
of the variable on the left, the expression on the right takes on the type of the variable on 
the left. Table 3-2 summarizes the conversion for all combinations of expressions and variables. 


Conversion 


None 
Take least significant word. 


Take least significant 16 bits 
of the integer part. 


Truncate fractional part; 
convert to integer form, 


Truncate fractional part; 
convert to integer form. 


Truncate fractional part of 
real part; convert to integer 
form. 


Extend sign bit to fill most 
significant word. 


None 


Take least significant 16 bits 
of the integer part. 


Truncate fractional part; 
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Table 3-2. Assignment Statement Type Conversion (Continued) 


Type of Variable 
(Left Side of 
Statement) 


INTEGER*4 


FIXED 


FIXED 


FIXED 


FIXED 


FIXED 


FIXED 


REAL 


REAL 


REAL 


REAL 


REAL 


REAL 


DOUBLE PRECISION 
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Type of Expression 
(Right Side of 
Statement) 


COMPLEX 


INTEGER*2 


INTEGER*4 


FIXED 


REAL 


DOUBLE PRECISION 


COMPLEX 


INTEGER*2 


INTEGER*4 


FIXED 


REAL 


DOUBLE PRECISION 


COMPLEX 


INTEGER*2 


Conversion 


Truncate fractional part of 
real part; convert to integer 
form. 


Integer becomes integral 
part; shift according to 
scale. 


Integer becomes integral 
part; shift according to 
scale. 


None, except possible shift. 


Convert both integral and 
fractional parts to integer 
form; shift according to 
scale. 


Convert Doth integral and 
fractional parts to integer 
form; shift according to 
scale, 


Convert both integral and 
fractional parts of real part 

to integer form; shift according 
to scale, 


Convert to floating-point 
form. 


Convert to floating-point 
form. 


Convert to floating-point 
form. 


None 


Take first two words of 
extended floating-point. 


Take real part. 


Convert to extended 
floating-point form. 


Table 3-2. 


Type of Variable 
(Left Side of 
Statement) 


DOUBLE PRECISION 


DOUBLE PRECISION 


DOUBLE PRECISION 


DOUBLE PRECISION 


DOUBLE PRECISION 


COMPLEX 


COMPLEX 


COMPLEX 


COMPLEX 


COMPLEX 


COMPLEX 
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Assignment Statement Type Conversion (Continued) 


Type of Expression 
(Right Side of 
Statement) 


INTEGER*4 


FIXED 


REAL 


DOUBLE PRECISION 


COMPLEX 


INTEGER*2 


INTEGER*4 


FIXED 


REAL 


DOUBLE PRECISION 


COMPLEX 


Conversion 


Convert to extended 
floating-point form. 


Convert to extended 
floating-point form. 


Extend single-precision 
floating-point to extended 
form. 


None 


Extend real part to extended 
floating-point form. 


Convert to real form; assign 
to real part; assign Zero to 
imaginary part. 


Convert to real form; 
assign to real part; assign 
zero to imaginary part. 


Convert to real form; assign 
to real part; assign zero to 
imaginary part. 


Assign to real part; assign 
zero to imaginary part, 


Take first two words of 
extended floating-point; 
assign to real part; assign 
zero to imaginary part. 


None 
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3.6.2 Character Assignment Statements 


Character assignment statements are used to store text strings in variables of type character, 
The form of a character assignment statement follows. 


FORMAT 
character variable = character variable or constant 


The variable on the left side of the equal sign may be a character scalar variable or a character 
array element reference. The right side of the equal sign may contain a character scalar, 
character array element, or a character constant. If the item on the right side has fewer 
characters than the assignment variable holds, the right-most remaining character positions 

are blank padded. If the item on the right has more characters than the assignment variable 
can hold, the right-most character positions are truncated before assignment. 


EXAMPLE 
A = 'GREEN' 


BUFFER = LINE 
TABLE(I) = '*' 
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4 CONTROL 


4.1 General 


Statements in a FORTRAN program are normally executed sequentially. This sequence 
may be altered by the use of control statements, which allow the programmer to make 
decisions during the course of his or her program and to control the order of program 
execution. 


FORTRAN statements may be given numbers to be referenced by control statements. A 
statement number is written as an unsigned, nonzero integer of five or fewer digits. 

Leading zeros are ignored. Although statement numbers are written as integers, they 
represent labels rather than numerical quantities, Statement numbers are used for program 
control, not numerical calculation, and must be unique within any separately compiled program 
unit. For example, no two statements may have the same numer in a program; however, 

a program may contain a statement numbered j and may call a subprogram which contains 

a statement numbered j. 

4.2 GO TO Statement 


The following paragraphs discuss the various forms of the GO TO statement and a statement 
used in conjunction with it, the ASSIGN statement. The statements included are as follows. 


° Unconditional GO TO statement 
° Computed GO TO statement 
° ASSIGN statement 


° Assigned GO TO statement 


4.2.1. Unconditional GO TO Statement 


FORMAT 


GO TO} 


This statement transfers control from one point of a program to another by indicating the 
number of the statement (j) to be executed next. Each time the GO TO statement is reached, 
the program transfers control to statement number j and proceeds sequentially from there 
until another control statement is reached. 
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EXAMPLE 


GO TO 10 
& A=B+C 
ic B=B+1 


Statement 10 is executed immediately after the GO TO statement. 


4.2.2 Computed GO TO Statement 
FORMAT 

GO TO (j qh 2reesir) [,] i 
This statement allows control to be transferred to any one of n different statements where 
Jqlzeeeegip, are the statement numbers, and i is a nonsubscripted integer variable which has 
a value in the range 1 through n. The computed GO TO causes control to be transferred to 
the statement whose label is j,, where subscript i is the value of the integer variable i, | f 
i is outside the range 1 through n, then the next sequential statement following the GO TO 
statement is executed. The comma separating i from the statement list is optional. 
EXAMPLE 

GO TO (89,14,92,2 1),ICE 


Statement 89 is executed next if the value of ICE is 1. If the value of ICE is 3, statement 92 
is executed, 


4.2.3 ASSIGN and Assigned GO TO Statements 
FORMATS 


ASSIGN i TO variable 


GO TO variable 


GO TO variable, (j ysl) 3e=esin) 
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Variable is a nonsubscripted integer variable, i is an executable statement number, and 
jyigeeessin are executable statement numbers, 


NOTE 


Another use of the ASSIGN statement is to assign a FORMAT statement 
label to an integer variable and then use that variable name in place of 
a format reference in a formatted I/O statement. Section 5 explains this 
feature in more detail. 


The second form of the assigned GO TO is allowed for coding compatibility with other 
versions of FORTRAN only. It is important to note that the list of statement numbers is 

not used to limit transfer of control, i.e., the statement number does not have to be included 
in the list for the program to work. However, the list may be used as a handy reference of 
all possible transfer points, and if a listed statement number is not used as a label some- 
where in the program, an error message is printed by the compiler. 


This statement transfers control to the statement whose number was assigned to the 
variable. The assignment must take place in a previously executed ASSIGN statement. 


The variable is a control variable, having a label as a value, not a numerical quantity. At 
the time of execution of an assigned GO TO statement, the current value of the variable 
must have been defined by the previous execution of an ASSIGN statement. The value of 
the integer variable is not the integer statement number; ASSIGN 10 to | is not the same 
as!= 10. 


EXAMPLE 
ASSIGN 40 TO NERROR 
CO TO NERROR 
40 A=A+B 


In the preceding example control is transferred to the statement numbered 40. 


EXAMPLE 


GO TON, (10,25,8) 
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If the current assignment of the integer variable N is statement number 8, then the 
statement numbered 8 is executed. If the current assignment of N is statement number 10, 
then the statement numbered 10 is executed next. If N is assigned statement number 25, 
statement 25 is executed next. If N is assigned statement number 12, statement 12 is 
executed, even though it is not included in the list within parentheses. 


EXAMPLE 


ASSIGN 10 TO ITEM 


13 GO TO ITEM , (8, 12, 25,50, 10) 


8 A=B+C 
10 B=C+D 
ASSIGN 25 TO ITEM 


GO TO 13 


25 C=Exe2 


In the previous example the first time statement 13 is executed, control is transferred to 
statement 10. On the second execution of statement 13, control is transferred to 
statement 25. 


4.3 IF Statement 


The following paragraphs describe the three forms of the IF statement: the arithmetic IF, 
the logical IF, and the block IF. 
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4.3.1 Arithmetic IF Statement 
FORMAT 


IF (e)n pn 23 


The quantity e is an arithmetic expression and ny,no, and n3 are statement numbers, The 
simple GO TO statement causes an unconditional transfer of control to a specified statement. 
The transfer does not depend on any condition of the data, status of the machine, or any 
other changeable hardware or software condition. The unconditional GO TO by itself permits 
little work to be done. It is also necessary to be able to transfer control when some condi- 
tion is met during program execution. The arithmetic IF statement accomplishes this. Wnen 

the value of the expression within parentheses is negative, the statement having statement 
number ny, is executed next; if the value of the expression is zero, statement no is executed 
next ; and if the expression is positive, n3 is executed next. 


EXAMPLE 


IF (ACJ, K)##3-B) 10, 4,30 
4 =B+C. 
30 C=(D+C ) ##2 


10. E=(F"B)/(D+1) 


In this example, if the value of the expression (A(J,K)**3-8) is negative, statement 
number 10 is executed next. If the value of the expression is zero, statement number 4 is 
executed next. If the value of the expression is positive, statement number 30 is executed 
next. 
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4.3.2 Logical IF Statement 
FORMAT 


IF(e)s 


Another tool for transfer of control is the logical IF statement, where e is a logical 
expression and s is any other statement except another logical IF, a block IF, or a DO 
(discussed in paragraph 4.4). The most common form of logical expression in this context 
is one that asks a question about two arithmetic expressions. Such relational expressions 
are written by using the six relational operators: .GT., .GE., .LT., .LE., .EQ., or .NE.. 
(These operators are discussed in Section 2). 


The action of the logical IF is as follows: if the logical expression is true, statement s is 
executed; if the logical expression is false, statement s is not executed. In either case, 
the next statement executed is the one following the logical IF, unless s is a GO TO or an 
arithmetic IF, and the expression is true. 


EXAMPLE 


IF (A. LE.0}G0 TO 25 
C=D+E 


25 WaX*#Z 


If the value of the expression (A.LE.0) is true, indicating that A is less than or equal to 0, 
the statement GO TO 25 is executed next and control is passed to statement number 25. 


If the value of (A.LE.0) is false, indicating that A is greater than 0, the statement GO TO 25 


is ignored and control is passed to the next sequential instruction. 


In the following example, assume that P and Q are logical variables. 


SSS 
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EXAMPLE, 


IF (P. OR. . NOT. G)A=B 
C=Bex2 


In the first statement, if the value of (P.OR..Not.Q) is true, the value of A is replaced by the 
value of B and the second statement is executed next. If the value of (P.OR..NOT.Q) is false, 
the statement A = B is skipped and the second statement is executed. 

4.3.3 Block IF, ELSE IF, ELSE, and END IF Statements 

The block IF statement is used with the END IF statement and optionally the ELSE IF and 

ELSE statements to control the execution sequence of a program. The form of a block IF 
statement follows. 


FORMAT 


IF (e ) THEN 
In the above statement, e is a logical expression. The form of an END IF statement follows, 


FORMAT 


END IF 
The form of an ELSE IF statement follows. 


FORMAT 


ELSE IF (e ) THEN 


In the above statement, e is a logical expression. The form of an ELSE statement follows, 


FORMAT 


ELSE 
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There must be a corresponding END IF statement for every block IF statement in a 
FORTRAN program. The IF level of a given statement s is defined to be the difference 
between the number of block IF statements up to and including s, and the number of 
END IF statements up to and not including s. The IF level of every statement in a program 
must be greater than or equal to Zero. The IF level of the END statement in every pro- 
gram must be equal to zero. The IF level of every ELSE IF, ELSE, and END IF statement 
in a program must be greater than Zero, 


EXAMPLE 
C T=0 
1 IF CI. EG. O} THEN 
2 J=I 
2 IF (J. EG. 0} THEN 
2 STOP 
2 ELSE IF(J. GT. G) THEN 
2 =-J 
2 ELSE 
2 I=0 
2 J=0 
2 ENDIF 
| J=0 
t ENDIF 
0 END 


Execution of a block IF statement causes evaluation of the expression e. If the value of 

e is true, normal execution continues with the statements following the block IF statement 

up to the next ELSE IF or ELSE statement that is at the same IF level as the block IF 

statement. At that point control goes to the END IF statement that corresponds to the block 

IF statement. If the value of e is false, control is transferred to the next ELSE IF, ELSE or END 

IF statement that is at the same IF level as the block IF statement. Execution of an ELSE IF 
statement is identical to the execution of a block IF statement. Execution of an ELSE statement 
causes control to proceed to the statement following the ELSE statement. Execution of an 

END IF statement also causes control to proceed to the statement following the END IF statement. 
Although it is legal to have statement labels on ELSE IF, ELSE, and END IF statements, their \) 
statement labels may not be referenced anywhere in the program. It is also illegal to use 

the statement label of a block IF statement as the terminal statement of a DO loop. 


4.4 DO and CONTINUE Statements 


The following paragraphs describe the DO statement and a statement used in conjunction 
with it, the CONTINUE statement. 
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4.4.1 DO Statement 


FORMAT 


DO nI,] v = pyPo,L p3) 


Most of the computation and information processing involved in a program is essentially 
repetitive. The DO statement serves as a powerful tool in FORTRAN by facilitating the 
definition and control of repetitive processes. For example, the following statements, which 
calculate the squares of 50 numbers, 


I=1 

10 SQCIP=ACID#ACT) 
I=I+1 
IF(I. LE. 50:¢60 TO 10 


can be replaced by 


DO 106 I=1,50 
16 SGC IS=ACIS#ACTI) 


The DO statement simplifies the iterative process controlled by the first, third, and fourth 
statements. 


In the general form of the DO statement, n is the statement label of an executable state- 
ment. The statement identified by n is called the terminal statement of the DO loop and must 
follow the DO statement. The range of the DO loop consists of all of the executable state- 
ments that appear after the DO statement, including the terminal statement of the DO loop. 
The variable v is a nonsubscripted integer, extended integer, real, or double-precision 
variable, called the DO variable; and py, pp, and p3 are integer, extended integer, real, 
or double-precision expressions. 


The expression p4 is called the initial parameter, and po is called the termination parameter. 
The expression p3 is the incrementation parameter and cannot equal Zero. If p3 does not 
appear, it defaults to a value of 1. 


The DO variable v may be used as a subscript or an ordinary variable within the DO loop. 
It may not be changed within the range of the DO loop, Also the DO variable retains its last 
defined value on exit from the DO loop. 
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The execution of a DO statement causes an iteration count to be established. The initial 
parameter m4, the terminal parameter m3, and the incrementation parameter m3 are produced 
by evaluating py,p2, and p3 respectively. If necessary, pyp2, and p3 are converted to the 
same data type as the DO variable. The iteration count is defined as follows. 


MAX( INT( (mo = Mat m3)/m3), 0) 


Note that the iteration count is zero whenever 
m4>mz2 and m3>0, or 


m4< mg andm3< 0. 


The maximum number of iterations is 32767; therefore, the execution of a DO loop will be 
unpredictable if the number of iterations calculated exceeds 32767. 


The DO variable v is defined with the value of the initial parameter m. On entry to the range 

of the DO, the iteration count is decremented and tested to see if its value is less than zero. 

If the iteration count is less than zero, control passes to the statement immediately following 

the terminal statement of the DO loop. If it is not, the range of the DO loop is executed, 

Note that the DO loop will never be executed if the initial value of the iteration count is zero, 

After the statements in the range have been executed, the value of v is incremented by the | 
incrementation parameter m3. Other than normal termination, a DO loop may become inactive | 
if a transfer of control occurs to an area outside the range of the DO loop, or if a RETURN 

or STOP statement is executed within the range of the DO loop. 


EXAMPLE | 
| 

DO 30 I=20,1.-1 

‘ 

| 


30 CHAxB. 


In this example, statements 1 through 5 will be executed 20 times with | ranging from 20 to 1. 
After completion of the DO loop, | will equal 0 and statement 6 is executed, 
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EXAMPLE 
I=1i 
DO 25 A=0,i1,.1 
RESULT (CID=SINCA?} 
25 I=I+i 


=B+C 


In the preceding example, the DO variable A is set to the initial value of 0, and the iteration 
count is computed to be 11. On entry to the range of the DO loop, the iteration count is 
decremented by one. Before the execution of the terminal statement, A is increased by the 
increment, .1, and the terminal statement is executed, After 11 executions of the DO loop, 
A equals 1.1. Since the iteration count is now Zero, control passes out of the DO loop and 
the fifth statement is executed next. 


EXAMPLE 
DO 25 M=11,10,1 
J=MoK 
25 = ARRAY (J) =BRAY(J? 
A=B+C 


In the preceding example, the DO loop will not be executed since the iteration count is 


zero. Control will pass to the statement following the terminal statement and M will have a 
value of 11. 


Several programming considerations for using a DO loop are discussed in the following 
paragraphs, 


The parameters of a DO statement (p yPoP3) may be changed (provided they are not 
constants) within the range of the DO loop. Since an iteration count has already been 
calculated, changing these parameters will not affect the execution of the DO statement. 
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The incrementation parameter p3is checked for a Zero value at compile time and at run 
time. If p3 is a constant, a syntax error will be generated for a Zero value at compile time. 
If the value of p. becomes Zero at run time, the DO statement is not executed, and the DO 
variable is not changed, 


If a DO statement appears within the range of another DO loop, the DO loop specified by 
that DO statement must be contained within the range of the outer DO loop. More than one 
DO loop may have the same terminal statement. A set of DO statements satisfying this rule 
are called nested DOs. 


EXAMPLE 
BDO 50 I=1,4 dr 
ACI}=BC( I) ##2 i Range of 
DO $0 J=1,5 “| Range of { Outer DO 
SC CC StL=ACT} i Inner BO { 
~ ~ 
EXAMPLE 
DO 10 INDEX=L,M ws 
N=INDEX+K i Range of 
DO 15 J=1,100,2 “{ Range of t Outer DO 
is TABLE ( J} =SUM(U, N}-1 t Inner DO H 
ic B(N}=A(N) ~ f 


If a DO statement appears within an IF block, ELSE IF block, or ELSE block, the range of 
the DO loop must be contained entirely within the block. 


If a block IF statement appears within the range of a DO loop, the corresponding END IF 
statement must appear within the range of the DO loop. 


The extended range of a DO statement is defined as those statements that are executed 
between the transfer out of the innermost DO within a nest of DOs and the transfer back 
into the range of this innermost DO. The extended range must be within the program unit 
containing the DO statement. The following restrictions apply. 


© Transfer into the range of a DO loop is permitted only if such a transfer is from 
the extended range of a DO loop. 


° The extended range of a DO statement must not contain another DO statement 
having an extended range. 
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The statement label of a terminal statement may not be used in any GO TO 
or arithmetic IF statement that occurs anywhere but in the range of the most 
deeply nested DO loop. 


In the following example, the transfers specified by the numbers 1, 2, and 3 

are permissible, whereas those specified by 4, 5, and 6 are not. Transfers 1 and 

3 are from a point in an inner DO loop to a point in a DO containing the inner 

DO. Transfer 2 is contained within the inner DO, Transfer 4 is from a point ina 

DO loop to a point in an inner DO loop skipping the initialization code of the 

inner DO which is not allowed. Transfer 5 is from a point in a DO loop to a terminal 
statement shared with a contained DO loop; this transfer is not allowed. Transfer 6 
is into a DO loop from outside its range; it also is not allowed, 


EXAMPLE isle) jate) 
Bo BO 
1 4 
DO 
5 


The terminal statement of a DO loop cannot be an unconditional GO TO, assigned GO TO, 
arithmetic IF, block IF, ELSE IF, END IF, RETURN, STOP, END, or DO statement. If the 
terminal statement is a logical IF statement, it may contain any executable statement except 
a DO, block IF, ELSE IF, END IF, END, or another logical IF statement. 


DO 5 K=i,4 
IFCXCK). GT. YC KD) YOK =X CK} 
CONTINUE 
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The second statement is executed four times regardless of the outcome of the IF statement. 
The third statement is executed when the DO loop becomes inactive regardless of the out- 
come of the IF statement. 


The execution of a subprogram reference within the range of a DO loop is permitted, 


4.4.2 CONTINUE Statement 
FORMAT 


CONTINUE 


The CONTINUE statement causes no instructions to be executed; it serves only as a dummy 
statement. In the normal sequence of a program, the statement following the CONTINUE 
statement is executed next. 


The CONTINUE statement is primarily used as the terminal statement for the DO loop. It 
must be used when the loop would otherwise terminate on an illegal terminal statement as 
specified in paragraph 4.4.1. The CONTINUE statement also acts as a transfer point for IF 
and GO TO statements within the range of a DO loop that are intended to begin another 
iteration of the loop. 


EXAMPLE 
DO 30 I=1,20 
7 IF(ACI)-B(I)) 5,30, 30 
5 A(I)mB(I)+1.0 


BC Id=BCI)-2.0 


60 To 7 
30 CONTINUE 
C=A(3)+B(7) 
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In the preceding example, the CONTINUE statement is used as the last statement in 
the range of the DO loop in order to avoid ending the DO loop with a GO TO statement. 


EXAMPLE 
DO 30 I=1,20 
IF(ACI)-B(I)) 5, 40,40 
5 ACI }=C(I) 


GO TO 30 
40 ACT}=0. 0 
30 CONTINUE 


In the preceding example, the CONTINUE statement provides a branch point enabling the 
programmer to bypass the execution of statement 40. 


The following paragraphs describe the CALL statement and the RETURN statement. 


4.5.1 CALL Statement 
FORMAT 

CALL S(€ 4,2 dye00,2r) 
The CALL statement causes a transfer of control to a subroutine subprogram (described in 
Section 6); s is the symbolic name of the subprogram and a 4,8>,..a, are arguments that 
may be required by the subprogram. The symbolic name of the subprogram is not assigned 
to a type and has no relation to the type of the arguments. If the subprogram requires no 
arguments, 24,8d,.e0.,2, are not required in the CALL statement. 
EXAMPLES 

CALL MULTI (A,B,C) 


CALL PRODUCE (A, 10,50,P) 
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The CALL statement causes control to transfer to the first executable statement of the 
subroutine subprogram. Arguments appearing in a CALL statement may be any constant, 
any subscripted or nonsubscripted variable, any expression, or a subprogram name. A 
subroutine subprogram may use one or more of its arguments to return results to the 
calling program. A constant should not be used as an argument If a value is being returned 
to that argument by the subroutine. The following example illustrates a problem that might 
be encountered when a constant is used as an argument. 


EXAMPLE 
Main Program Subroutine Subprogram 
CALL MIN(6, 3, 25?) SUBROUTINE MIN (N,M, 1} 
- : M=N% I 
RETURN 


20 «g=gt END 


Here the subroutine replaces the constant 3 with the value N * |, which is 150 in this case. 
This results in the variable S being assigned a value of 150 when statement 20 is executed, 
rather than the expected value of 3. 


4.5.2 RETURN Statement 


FORMAT 


RETURN 


This statement returns control from a subprogram to the calling program. Normally, the last 
statement executed in a subprogram is a RETURN statement. It need not be physically the 
last statement of the subprogram. Any number of RETURN statements may be used. RETURN 
may appear only in subprograms. 


4.6 PAUSE, STOP, END, and PROGRAM Statements 
The following paragraphs describe the PAUSE statement, STOP statement and END statement. 


—— LLL enn 
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4.6.1 PAUSE Statement 
FORMAT 

PAUSE [c] 
In this statement, c is an integer constant in the range 0S cS 99999, or else c is a character 
constant. When c is omitted, a character constant of length 0 is assumed, 
This statement is included for compatability with other versions of FORTRAN. The message 
'PAUSE' followed by the constant c is displayed to the FORTRAN error logging device or 
file. If a character constant with length greater than 74 is used, the characters past the 74th 
character position are truncated from the display. 
After displaying the pause message, normal program execution resumes. If a user wishes 
to halt program execution pending operator response he should use the ACCEPT statement 
for a video display terminal, or else a READ statement for a hard-copy input device. 


EXAMPLE 


PAUSE 25 


4.6.2 STOP Statement 
FORMAT 

STOP [ce] 
In this statement, c is an integer constant in the range 0S cS 99999, or else c is a character 
constant. When c is omitted, a character constant of length 0 is assumed, 
The message 'STOP' followed by the constant c is displayed to the FORTRAN error logging 
device or file when this statement is executed. If a character constant of length greater 
than 74 is used, the characters past the 74th character position will be truncated from the 


display. 


After displaying the stop message, program processing stops, and control returns to the 
operating system. 


EXAMPLE 


STOP 'THE MAG TAPE WAS NOT READY' 
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4.6.3 END Statement 
FORMAT 


END 


The END statement is an executable statement that defines the end of a source program 

or source subprogram for the compiler. Physically, it must be the last statement of each 
program or subprogram. If program execution reaches the END statement in a main program, 
execution terminates and control returns to the operating system. If program execution 
reaches the END statement in a subprogram, control returns to the calling routine. The END 
statement must be on a single source line; continuation lines are not allowed. The END state- 
ment may have a statement number. 


4.6.4 PROGRAM Statement 


FORMAT 
PROGRAM N 


If the PROGRAM statement is used, it must be the first statement in the main program (other 
than comments). 


The letter N is an identifier and must follow the same syntax rules that are discribed in para- 
graph 2.2.1. 
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5 INPUT/OUTPUT STATEMENTS 


5.1 General 


The input/output (I/O) statements control the transmission of information between the 
computer and the I/O devices or units, |/O statements are classified as follows: 


° FORMAT Statement. The FORMAT statement is a nonexecutable statement 
that specifies the conversions required between internal and external data 
forms. It describes how the data is arranged for input or output. Data trans- 
ferred without a FORMAT statement is unformatted and reflects the data as it 
appears in computer memory. 


° READ and WRITE Statements. READ and WRITE are statements that specify 
transmission of information between computer memory and various input/output 
devices. 


° Mass Storage File Input/Output Statements. REWIND, BACKSPACE, and END FILE 
provide positioning and file termination for magnetic tape and disk 1/0. 


NOTE 

All devices are opened on first |/O operations, not at initialization, 

A REWIND may be used to open a device prior to the first |/O. 
5.2 FORMAT Statement 
All formatted input or output requires the use of a data format specifying the external format 
of the data and the type of conversion to be used, The data format is given in a FORMAT 
statement or as a character string in a constant or variable, 


FORMAT 


FORMAT (S4 S2,e00,Sp) 


In the above statement, S is a data field specification as defined in the following paragraphs. 
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5.2.1 Numerical Data Specifications 


Conversion of arithmetic data may be one of the following six types: 
° D field specification, The internal form is double-precision, binary floating-point. 
The external form is double-precision, decimal floating-point. The specification 
has the following form. 
Dw.d 


° E field specification. The internal form is binary floating-point. The external form 
is decimal, floating-point. The specification has the following form, 


Ew.d [Ee] 


° F field specification. The internal form is binary floating-point. The external form 
is decimal fixed-point. The specification has the following form. 


Fw.d e 


° G field specification. The internal form is binary floating-point. The external form 
is decimal fixed-point or floating-point. The specification has the following form. 


Gw.d [Ee] 


° | field specification. The internal form is binary integer. The external form is 
decimal integer. The specification has the following form. 


lw 


° Z field specification, The internal form is binary, The external form is hexadecimal 
integer. The specification has the following form. 


Zw 
The letter D, E, F, G, |, or Z designates the conversion type; w is an integer specifying the 
field width, which includes preceding bianks, minus sign, decimal point, and exponent; dis Qy 
an integer specifying the number of decimal places to the right of the decimal point; e is 
an integer specifying the number of digits in the exponent. For example, the statement 


FORMAT (15,F 10.2,D25.8) 


could be used to output the line 


jefe] La[ 5] ef es [oops eafeal safeafeef+rfeafvofaolas aafaa[zqfan aef27]20]20]20]>:|32] 20] 24[ 9] 222] 4 0 eo 
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The field width, w, includes all characters in the number. If the number is too large for 
the field width, the field is filled with asterisks, If the number is smaller than the field width, 
it is rightjustified within the field. 


The type of conversion used must correspond to the type of variable in the input/output 
list. | conversion is used for integer or extended integer variables; E, F, or G conversion is 
used for real variables; the D conversion is used for double-precision variables; and the 
Z conversion is a hexadecimal conversion for any numeric variable type. 


5.2.1.1 |Field Specifications 


The form of the | field specification is shown in the following. 
FORMAT 


lw 


The | field specification is used to input or output an integer value, The form lw specifies 
that w character positions are to be allotted to the number, The number is right-justified 
within the field and unused space remains blank. If a value requires more than w positions, 
an error condition exists and is indicated by w asterisks placed in the field, If the value is 
negative, a minus sign immediately precedes the left-most digit. If the value is positive, 

a blank immediately precedes the left-most digit. A position must be reserved for a minus 
sign, if any; thus, a negative number with a size greater than w-1 digit positions results in 
an error condition. Blanks received as input are interpretes as Zeros unless a BN edit 
descriptor is in effect (refer to paragraph 5.2.3.5). The following examples show how each 
of the values on the left is printed for the format specification 14: 


EXAMPLES 
Value Stored internally Appearance on Printer Line 
459 459 
+459 459 
-360 ~360 
10 10 
10000 eens 


5.2.1.2 F Field Specification 


The form of the F field specification is shown in the following. 
FORMAT 


Fw.d 
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The F indicates conversion between an internal real value and an external number without 
an exponent, The quantity w refers to the length of the total field and d refers to the number 
of digits which are to appear to the right of the decimal point (the fractional portion). The 
total field length must provide positions for a sign, if any, a decimal point, and a digit to the 
left of the decimal point. The sign is printed only if the number is negative. In general, w 
should be greater than or equal to d + 3 for output. 


If a fractional portion requires more positions than are reserved by d, it is truncated from the 
right. If the positions reserved by d are not filled, Zeros fill to the end of the field, The integer 
portion of the field is treated in the same manner as numbers converted by the | field speci- 
fication on both input and output. Blanks received as input are interpreted as zeros unless 

a BN edit descriptor is in effect (refer to paragraph 5.2.3.5). 


The following examples show how each of the values on the left is printed according to the 
specification F7.3. 


EXAMPLES 
Value Stored Internally Appearance on Printed Line 

32.447 32.447 

88.668 -88.668 

-.33 -0.330 

8.6452 8.645 

-3. -3.000 

8.064 8.064 

1325.943 nw 


Asterisks are placed in all error fields. 


For F field specification input, the decimal point in the data field is optional, If a decimal 
point is present that is inconsistent with the format specification, the actual position of the 
decimal point overrides that indicated in the format and governs the conversion. If no de- 
cimal point is present, space should not be reserved for it. The number is given a decimal 
point when converted internally; the format specification is used to place the decimal 

point in the correct position. For example, these two input fields yield the same result inter- 
nally if the specification is F7.3. 


246.135 


246135 
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For F field specification output, fractional numbers are usually printed with a zero 
immediately to the left of the decimal point. If d is zero, a fractional value is printed as a 
sign, a zero, and a decimal point. 


For example, for the specification F4.0 
-.3564 is converted to -0. 
A zero value is printed with a zero preceding the decimal point. 


Numbers in an E field specification format are accepted as data for F field specification 
input. 


5.2.1.3 E Field Specification 


The form of the E field specification is shown in the following, 
FORMAT. 


Ew.d [Ee] 


The E indicates conversion between an internal real value and an external number with 

an exponent. The d indicates the number of digits in the fractional portion, The e indicates 
the number of digits in the exponent and has no effect on input. If the Ee form is not used, 
the default value for e on output is 2 places, For output, the minimum requirement is the d 
field, a decimal point, a digit to the left of the decimal point, a place for a sign (or blank), 
and the exponent. Each of these must have an allotted space on output. If w<d+e+5an 
output error condition will occur. For input, space for all of these positions is not necessary, 
In general, the relationship w= d + 7 should be sufficient. 


The exponent is a signed or unsigned one- or two-digit integer constant preceded by the 
letter E and within the range -79 <exp <76, The integer indicates the power of 10 that is 
used to multiply the number to obtain its true internal value. Blanks received as input are 
interpreted as zeros unless a BN edit descriptor is in effect. 


The following examples show how each of the values on the left is printed according to the 
specification E11.4 or E11.4E2. 


EXAMPLES 
Value Stored Internally Appearance on Printed Line 
4381. 0.4381E 04 
-00000000004381 0.4382E-10 
262610. -0.2626E 06 
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Asterisks are placed in all arror fields. For example, an attempt to print the above values 
in E10.4 will result in a string of 10 asterisks on the printed line, 


For input, the start of the exponential field may be indicated by a + or - rather than an E, 
The following numerical fields would each yield the same result: 


1.25E4, 1.25E+4, 1.25+4 


E conversion will accept F field specification data as input. For example, an input field of 


12500. 


produces the same value internally as the examples immediately above, If the decimal point 
is omitted from the input field, the format specification will place the decimal point in the 
correct position. 


5.2.1.4 G Field Specification 


The form of the G field specification is shown in the following, 
FORMAT 


Gw.d[Ee] 


The G field specification is used for numbers that may be represented in either fixed or floating 
point. The format editor automatically selects one of the two representations for data output. 
The numeric field descriptor Gw.d indicates that the external field occupies w positions with d 
significant digits. The e indicates the number of digits in the exponent and has no effect on 
input. If the Ee form is not used, the default value for e on output is 2. Input processing is 

the same as for the F field specification. 


The method of representation in the external output string is a function of the magnitude | 
of the internal data item. The following tabulation exhibits correspondence between the output 


value N and the equivalent method of conversion that will be effected, ww 

Magnitude Equivalent 

of Conversion | 

Data Item Effected | 

0.1<N<1 F(w-4).d,4X | 

1<N<10 F(w-4).(d-1),4X | 

. . | 
10°2<N<10% F(w-4).1.,4X 
10*"*<N<10? F(w-4).0,4X 

Otherwise 


(N<0.1;N210*) sEw.d 
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In the conversion column, 4X represents four blanks. In the last conversion, s is the current 
scale factor (refer to paragraph 5.2.1.7), and E field conversion is used, Note that the effect 
of the scale factor is suspended when the magnitude of the data item to be converted is 
inside the range that permits effective use of F field conversion. 


When used for data input, values read into variables using a G field specification are converted 
according to the type of the corresponding variable in the input list. 


The following is an example of a G field specification and the resulting data format. 
EXAMPLE 


REAL*8 X(9) 
FORMAT (G13.6) 


External 
Internal Value Representation 
X(1)= .01234567 .123457D-01 
X(2) = -0.12345678 ~.123457 
X(3) = 1.23456789 1.23457 
X(4) = 12.34567890 12.3457 
X(5)= 123.45678901 123.457 
X(6) = ~1234.56789012 -1234.57 
X(7)= 12345.67890123 12345.7 
X(8)= 123456.78901234 123457. 
X(9) = -1234567.89012345 -.123457D +07 


5.2.1.5 D Field Specification 


The format of the D field specification is shown in the following. 
FORMAT 


Dw.d 


In D field conversion, the internal value must be double precision. The numeric field descriptor 
Dw.d indicates that the total field occupies w positions, and the fractional part consists of d 
digits. 


The basic form of the external input field is the same as for real conversions. Any of the 
following are valid inputs in a D field specification, and will produce the same result internally. 
3D1 

3E1 

30. 
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The external output field is the same as for E conversion, except that the character D replaces 
the character E in the exponent. 


5.2.1.6 Z Field Specification 
The form of the Z field specification is shown in the following. 


FORMAT 


Zw 


The Z field specification is used for the transmission of hexadecimal data. The input fields 

are scanned with any leading, imbedded, or trailing blanks treated as Zeroes. In internal 
storage, two hexadecimal digits occupy one byte of storage. If the input data is an odd number 
of digits, padding with zeros occurs on the left. When the input data is too large for the storage 
area, truncation of the high-order digits occurs on output. If the number of characters in 
storage is less than w, the left-most positions are filled with blanks, If the number of characlyss 
is greater than w, the left-most digits are truncated. 


5.2.1.7 Scale FaktorsA scale factor may be specified which applies to all D, E, and F field 
specifications, 

FORMAT 

nP 
In the preceding edit descriptor, n is a signed or unsigned integer specifying the scale factor, 
and P is the identifying character. The scale factor changes the position of the decimal point 
for both input and output operations. 
For F field conversion, the scale factor specifies a power of ten such that. 

external numbers=(internal number)*(power of ten) 7 
For D and E field conversions, the scale factor multiplies the number by a power of ten, but 
the exponent is changed accordingly, leaving the number unchanged except in form. For 
example, if the statement 

FORMAT (F8.3,E 16.5) 
corresponds to the output line 

26.451 -0.41321E-01 


then the statement 


FORMAT (-1P,F8.3,2P,E 16.5) 


n./G 


a 
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Corresponds to the output line 


When no scale factor is given, it is understood to be 0. However, once a scale factor is given, 
it holds for all following D, E, and F field conversions within the same format. The scale factor 
is reset to zero by giving a scale factor of zero. Refer to paragraph 5.2.1.4 for the effect 

of a scale factor on G field conversions. 


Scale factors have no effect on! or Z conversions. 
For input operations, the scale factor may be supplied for all real number inputs, but is ignored 
for all inputs containing an exponent in the external field, When used with input data, the 
direction of the decimal point movement is opposite that for an equivalent output specification, 
For example, if the statement 

FORMAT (F8.3) 
inputs data to be used internally as 

26.451 
then the statement 


inputs data to be used internally as 


264.51 


5.2.2 Logical Data Field Specification 


Logical data can be transmitted in a manner similar to numerical data by use of the logical 
data format. The field specification for the format follows. 


FORMAT 


Lw 


In the above statement, L is the control character and w is an integer specifying the field 
width, 
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Data is transmitted as the value of a logical variable in the input/output list. 


On input, the first nonblank character in the data field must be T or F, and the value of the 
logical variable is stored as true or false, respectively, The remainder of the field is ignored. 
If the entire data field is blank, a value of false is stored, 


On output, w-1 blanks followed by T or F are output, depending on whether the value of 
the logical variable is true or false, respectively. 


5.2.3 Alphanumeric and Literal Data Specifications 


The following paragraphs describe printer carriage control characters, the alphanumeric data 
field specification, and the specifications for the H and X fields. 


5.2.3.1 Carriage Control 


The first character in a record to be printed on a listing device is used for carriage control. 
This character is not printed and is usually specified at the beginning of the format field 
specification for the record as 1Hx or x, where x is a blank, 0, 1, or +. The specification may 
also be a slash (/). The effect of each character is given in table 5-1. 


Table 5-1. FORTRAN Carriage Control Characters 


Carriage 
Control Characters Print Results 
Blank Skip one line before printing. 
0 Skip two lines before printing, © 
1 Skip to first line of next page 
before printing. 
+ Suppress line spacing before 
printing (overprint). 
/ Skip one line before printing 
does not require 1H prefix. 
EXAMPLE 
10 FORMAT(9H1 PAGE ,13/1H0) 


in this statement the character "1" following the 9H prefix is interpreted as carriage control 
and causes the literal message " PAGE " to be printed at the top of the next page. Similarly, 
the character "0" following the 1H prefix is interpreted as carriage control to skip two lines 
before the next line is printed. The slash (/) represents an end-of-record and causes a line 
feed before the 1HO specification is interpreted. The net result is three blank lines between 
the " PAGE " label output and any subsequent lines. 
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5.2.3.2 A Field Specification 


The specification for character data has the following form. 
FORMAT 


Alw] 


In the preceding edit descriptor, A is the field specification character and w is an optional 
field width. When w is omitted, it is assumed to be equal to the length of the character datum 
being transmitted. 


EXAMPLE 


CHARACTER A#25, B#15, C#80 
READ(S, 1034, B,C 
16 FORMAT (3A) 


The previous statements will cause A, B, and C to be read with field widths of 25, 15, and 
80 characters, respectively, from the input file. The effect is the same as if the following 
FORMAT statement were used, 


10 FORMAT (A2S, A1S, ASO} 


Character data can be transmitted in a manner similar to numerical data by use of the form 
Alw]. The alphanumeric characters are transmitted as the value of a variable in an input/ 
output list. The variable must be of type CHARACTER. For example, the sequence 


CHARACTER CODE#2 
READ(2, 5)CODE 
i] FORMAT (A2?) 


causes two characters to be read and placed in memory as the value of the character variable 
CODE, 


The number of characters transmitted is limited by the maximum number of characters which 
can be stored in the space allotted for the variable. If w exceeds the available space, leading 
characters are lost on input and replaced with blanks on output. When w is less than the 
available space, blanks are filled in after the given characters until the maximum is reached, 
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The following examples show the A field specification used for input and output data, 
EXAMPLES 
For Output: 


CHARACTER*4 ITEMP 
DATA ITEMP /ABCD/ 


Format External Field intemal ; 
Representation 
A6 ABCD ABCD 
A4 ABCD ABCD 
A2 AB ABCD 
For Input: 


CHARACTER*4 ITEMP 


internal 
Format Representation External Field 
Aé CDEF ABCDEF 
A4& ABCD ABCD 
A2 AB AB 


5.2.3.3 H Field Specification 


The form of the H field specification is shown in the following. 
FORMAT 

nH 
The H field specification specifies that an alphanumeric string is a Hollerith data field, H is 
the control character and n is the number of characters in the string, counting blanks, In 
all examples, b represents a blank. 

FORMAT (17HbDPROGRAMDCOMPLETE) 
can be used to output 

PROGRAM COMPLETE 


on the output listing. On input the externa! characters are stored in the format itself only 
if there is no list in the READ statement. 


f.Man 
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An alternate form of this field is to specify the alphanumeric string within apostrophes. The 
FORMAT statement above could have been written this way. 


FORMAT (‘PROGRAM COMPLETE") 


The input and output results are the same as the first FORMAT statement in this paragraph. 
Where an apostrophe appears in a character string of this nature, it must be represented by 
two apostrophes in succession ("). 


An alphanumeric format field may be placed among other fields of the format, For example, 
the statement 


FORMAT (15,8HbFORCEDb=F 10.5) 
can be used to output the line 
22 FORCE = 17.68901 


The separating comma may be omitted alfter an alphanumeric format field. 


5.2.3.4 X Field Specification 


The form of the X field specification is shown in the following, 
EORMAT 


nX 


The X field specification allows blanks to be introduced into an output record or characters 
skipped in an input record. The control character is X, and n is the number of blanks or 
characters skipped. The quantity n must be greater than zero. For example, the statement 


FORMAT (5H STEP 15.10X,3HY =F7.3) 
may be used to output the line 
STEP 28 Y =453.872 


Ten blanks separate the two quantities. 


5.2.3.5 BNand BZ Editing 


The BN and BZ edit descriptors may be used to specify the interpretation of blanks other 

than leading blanks in numeric input fields. At the beginning of execution of each formatted 
input statement, such blank characters are interpreted as zeroes. If a BN edit descriptor is 
encountered in a format specification, all such blank characters in succeeding numeric input 
fields are ignored. The effect of ignoring blanks is to treat the input field as if all the blanks 
had been removed. The remainder of the field is right-justified and padded with leading blanks. 
The only exception is a field of all blanks which would be considered a zero. This treatment 
continues for the duration of the format specification or until a BZ edit descriptor is scanned 
which reverts back to the original mode. 
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The following example right justifies the number 25 in the input field, 
EXAMPLE 


Input Field: tus2i5i_ 4-4 


c 
c A=0. 25 WITH BN 
c 


READ(S, 10034 
106 FORMAT (BN. FS. 2) 


Cc 4=25. 00 WITH EZ 


READ(S, 15034 
156 FORMAT (BZ, F5. 2) 


5.2.3.6 T Field Specification 
For tabular output, the T field specification can be used to designate the position in the output 


record where data transfer will begin. The form of the T field specification is shown in the 
following. 


FORMAT 


Tn 


T is the control character and n designates the character position within the output record 
that will be occupied by the first character of the transfer. For printer output, the first 
character of the output record is reserved for carriage control. Therefore, a one-character 
skew occurs between the character's position in the output record and its position on the 
printed line. 


EXAMPLE 

FORMAT(T20,'NAME') 
This statement results in printing the character string NAME beginning at printer column 
position 19. Also, because of the presence of the carriage control character, the following 
two statements produce identical output. 


FORMAT(T 1,6HDGRADE) 


FORMAT(T2,5HGRADE) 
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Both of these statements print the character string GRADE beginning at the left-most print 
position. 


EXAMPLE 
FORMAT(T 1,5HGRADE) 


The above example causes the format editor to interpret the letter G as a carriage control 
character and output the character string RADE beginning at the left-most print position. 


To create output in many columns, several T field specifications may be included in one 
FORMAT statement. The order that the specifications appear in the FORMAT statement is 
not important; that is, a T70 specification may appear in the statement before a T12 speci- 
fication without affecting the desired appearance of the output. 


EXAMPLE 
FORMAT(T20,'NAME',T40,'AGE',T 1,6HbGRADE) 


The above example prints the character string GRADE beginning at the left-most print 
position, NAME beginning at printed column position 19, and AGE beginning at printed 
column position 39. 


5.2.4 Complex Quantities 


Complex quantities are transmitted as two independent real quantities. The format specifi- 
cation is given as two successive real specifications or one repeated real specification 
using any combination of Fw.d, Ew.d, Dw.d, or Gw.d. 


On input, two successive fields are read and assigned to a complex input/output list element 
as its real and imaginary parts. 


EXAMPLE 
Format External Fields Internal Representation 
F8.5,F8.5 (1234567, 12345.67) 123.45678, 12345.67 
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On output, the constituent parts of a complex value are transmitted under the control of 
repeated or successive field descriptors, 


EXAMPLE 
Format Internal Values External Representation 
2F8.5 1.2345118, 2.345452 1.23451 2.34545 


5.2.5 Repeated Group and Field Specifications 


Repetition of the same field specification may be indicated conveniently in a FORMAT 
statement. In addition, groups of repeated specifications may be represented conveniently. 


5.2.5.1 Repetition of Field Specifications 


Repetition of a field specification may be performed by preceding the control character 
D, E, F, G, 1, L, A, or Z by an unsigned integer giving the number of repetitions desired, 
For example, 

FORMAT (2E 12.4,315) 


is equivalent to 


FORMAT (E12.4,E12.4,15,15,15) 


5.2.5.2 Repetition of Groups 


A group of field specifications may be repeated by enclosing the group in parentheses and 
preceding the group with the repetition number. For example, 


FORMAT (218,2(E 15.5,2F8.3)) 
is equivalent to 
FORMAT (218,E 15.5,2F8.3,E 15.5,2F8.3) 


Up to two levels of parentheses are allowed in group repetition, in addition to those enclosing 
the entire format. For example, the statement 


FORMAT (15,2(E15.5,2(F8.3,F 7. 1))) 


exhibits the maximum level of nesting. The outermost parenthesis pair is termed level Zero, 
the next level one, and the innermost level two. 
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5.2.6 Record Separation Indicator (Slash) 


To specify a group of input/output records with different records having the specifications, 
use a slash "/" to indicate a new record. For example, the statement 


FORMAT (318/15,2F8.4) 
is equivalent to 


FORMAT (318) 
for the first record and 
FORMAT (15,2F8.4) 


for the second record. The separating comma may be omitted when a slash is used, Blank 
records may be written on output or records skipped on input by using consecutive slashes, 


Both the slash and the closing parenthesis at the end of the format indicate the termination 
of a record. If the list of an input/output statement dictates that transmission of data is to 
continue after the closing parenthesis of the format is reached, the format is repeated from 
the last left parenthesis of level one or zero. Level one is used when it exists; otherwise, 
level zero is used. Thus, the statement 

FORMAT (F7.2,(2(E 15.5,E 15.4),17)) 
causes the format 

F7.2,2(E 15.5,E 15.4),17 
to be used on the first record and the format 


2(E 15.5,E 15.4),17 


on succeeding records. 


As a further example, consider 

FORMAT (F7.2/(2(E 15.5,E 15.4),|7)) 
The first record has the format 

F7.2 
and successive records have the format 


2(E 15.5,E 15.4),17 
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5.2.7 Character Format References 


The alphanumeric string comprising a format specification may be stored in a character 
variable or constant. Input/output statements may reference the format by giving the con- 
stant or the variable name rather than the statement number of a FORMAT statement. The 
stored format must have the same form as a FORMAT statement, excluding the word FORMAT. 
The enciosing parentheses must be included. 


As an example, consider the following three sequences which all have the same execution 


effect. 
CHARACTER FORM#iS 
FORM=/(1X, 15, 3X, F9. 2) 
READ(5, FORM)K, X 


READ(S, ‘(1X, 15, 3X, F%. 2) 7K, X 


READ(5, 1000)K, X 
i000 FORMAT(1X, 15, 3X,F9. 2) 


The first example assigns to the character variable FORM a string which contains the desired 
format specifications, and then uses the variable name FORM in the READ statement to 
reference the format. The second example contains the format specifications in a character 
constant within the READ statement. The last example uses a regular FORMAT statement 

to specify the format information, and references the statement label of the FORMAT state- 
ment within the READ statement. 


5.2.8 Assigned Format Labels 


Another way to specify a format reference is to assign the statement label of a FORMAT 
statement to an integer variable in the same way an executable statement label may be 
assigned to an integer variable (see paragraph 4.2.3). After execution of the ASSIGN state- 
ment, the integer variable name may be used to designate the FORMAT statement whose 
label has been assigned to it. 


EXAMPLE 


ASSIGN 1000 TO LABEL 
READ(5, LABEL ?K: X 
1000 FORMAT(1X, 15, 3X. F9. 2) 
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5.2.9 Combinations of Formats 


The different types of field specifications described in the preceding paragraphs may be 
combined and used in a single FORMAT statement. An alphanumeric format field (H field 
specification) may be combined with other field specifications as described in paragraph 


625333. 


5.2.10 Free Field Format 


When performing a formatted READ with FORTRAN-78, the data can be arranged in a free 
field format that terminates each field with a comma. The comma denotes both the end of 
one field and the start of the next field. Data in the field is right-justified. To use this format 
the field specification must account not only for the positions used by the largest number, 
but also for the position used by the comma. For example, to place a 1- to 4~digit positive 
integer into five sequential fields using the free field format, the format statement must be 


changed to 
FORMAT(S5I5) 


to allow for the comma position. The input record can appear as 


3,52,538,6744,5313, 


The last comma in the input record is necessary to maintain the value of the parameter. 
Without that comma, the value is interpreted as a standard I5 format parameter and is 
assigned the value of 53130. The comma can be omitted if a Zero or blank precedes the 


5313 in the input record, i.e., 


»55313 


When the free field format is used with an ACCEPT statement, a carriage return may replace 
the last comma to terminate the data field and the input record. 


5.3 Sequential Access INPUT/OUTPUT 


Sequential access |/O statements allow the programmer to access records within a file in 
sequential order. The characteristics of sequential access files are: 


° If the file contains one or more records, those records exist as a totally ordered set. 


° There exists a unique position of the file, called its initial point. If a file contains 
no records, the unit is positioned (physically or logically) to write, starting at the 
initial point. If the unit is positioned at the initial point and the file contains records, 
the first record of the file is defined as the next record. No record precedes the 


initial point. 
° If a unit is positioned at a point of the file beyond the initial point, a unique pre- 
ceding record and a unique next record are associated with that position. 


° Upon completion of the execution of a WRITE or END FILE statement, there exist 
no records following the records created by that statement. 


° When the next record is transmitted, the position of the records being read or 
written by the unit is changed so that the record just transmitted becomes the 


preceding record, 
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The sequential access I/O statements are READ. WRITE. BACKSPACE. REWIND. and END FILE. 
5.3.1 External File Characteristics 


An external file consists of one or more records, which may be formatted or unformatted. 

A formatted record consists of a string of printable characters. The transfer of such a record 
requires that a format be referenced to supply the necessary positioning and conversion 
specifications, The number of records transferred by the execution of a formatted READ or 
WRITE is dependent upon the list and format specification. An unformatted record consists 
of binary data. When an unformatted or formatted READ statement is executed, the required 
records on the identified unit must be unformatted or formatted records, respectively. 


5.3.2 Internal File Characteristics 


When a formatted READ or WRITE statement specifies an internal file name in place of a 
logical unit number, the I/O transmission is performed internally between memory locations. 
An internal file name is either a character scalar or character unsubscripted array name. If 
a scalar variable is specified, the internal file consists of one record whose length is the de- 
clared length of the character variable. If an array name is used, the internal file consists 

of as many records as there are elements in the array, each record having a length equal 
to the declared character length of the array elements. A multirecord internal file is always 
positioned at the beginning of the first record for each READ or WRITE operation referencing 
it. The 'ERR=label' option may be used in an internal READ or WRITE statement to provide 
error handling for format errors. Internal file 1/O is the functional equivalent of the ENCODE 
and DECODE statements supported under previous releases of FORTRAN. The statements 
ENCODE and DECODE are still recognized by the compiler; however, the buffer variable 
used must be declared to be of type character. 


EXAMPLE 


CHARACTER#&80 INTERN(2) 
WRITE(C INTERN, 1) IANSi, IANS2 
i FORMAT(’ FIRST ANSWER IS’, 15, / 
* ‘ SECOND ANSWER IS’, 15) 


The preceding example will assign the string 'FIRST ANSWER IS' to the first element (record) 
of array INTERN, followed by the five-digit representation of the value in IANS1, followed 
by 60 trailing blanks. The string ‘SECOND ANSWER IS' is assigned to the second element 
(record) of array INTERN, followed by the five-digit representation of the value in IANS2, 
followed by 59 trailing blanks. Note that the terms record and element are equivalent when 
referring to multirecord internal files. 


5.3.3 READ and WRITE Statements 


In general, a READ or WRITE statement provides: 


© Specification of the operation required, by using the verbs READ or WRITE to 
indicate input or output, respectively. 
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© A unit specifier which may be either a logical! unit specifier, or an internal file 
name. A logical unit specifier is used to perform |/O with some external device 
or file and may be an unsigned integer constant or an integer variable name. An 
internal file name is either the name of a character scalar variable, or a character 
array name. 


© Reference to a data format that specifies the types of conversions required between 
the I/O list and the formatted record(s). See paragraph 5.2 for a description of 
format references. 


© A list of the variables whose values are to be transmitted. The values are trans- 


mitted in the order in which information exists on the input record or will exist on 
the output record. 


EXAMPLE 
WRITE(1,5)A,X,K 
The preceding statement specifies that the values of A, X, and K, in that order, are to be 
written on logical 1/O unit 1 according to the format given in the FORMAT statement 
numbered 5. 
EXAMPLE 
BUFF = '00231' 
READ (BUFF,'(15)")J 
The preceding statements specify that variable J will receive a new value by interpreting 


the character information in character variable BUFF with an I5 format. Variable J becomes 
defined with the integer value 231. 


5.3.3.1 Input and Output Lists 


Each READ or WRITE statement includes a list of the names of variables, arrays, and array 
elements, The named elements are assigned values on input and have their values transferred 
on output. The list of a READ or WRITE statement specifies the order of transmission of variable 
values. During input, the new values of listed variables may be used in subscripts for variables 
appearing later in the list. For example 


READ(2,3)L,A(L),B(L + 1) 


reads a new value of L and uses this value in the subscripts of A and B. 
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The transmission of array variables may be controlled by indexing similar to that used in the 

DO statement. The list of controlled variables, followed by index control, is enclosed in 

parentheses and the whole acts as a single element of the list. For example 
READ(7,23)(X(K),K=1,4) 

is equivalent to 


READ(7,23)X(1),X(2),X(3),X(4) 


Like the DO statement, the initial limit, final limit, and increment values may be expressions 
such as 


READ(4,2)N,(GAIN(K),K=1,M+N+1,M) 

The indexing may be nested as in 
READ(2, 13)((MASS(K,L),K=1,5),L=1,4) 

The preceding statement reads in the elements of array MASS in the order 
MASS(1,1),MASS(2, 1) peoeyMASS((5, 1),MASS( 1,2) p200.,MASS(5,4) 


If an entire array is to be transmitted, the indexing may be omitted and only the array identi- 
fier written. The array is transmitted in order of increasing subscripts with the first subscript 
varying most rapidly. Thus, the example above can be written 


READ(2,13)MASS | 


5.3.3.2 Formatted READ Statement 


The purpose of the READ statement is to transfer information to a list of input variables from | 
either an external file or device, or an internal file. 


FORMAT 
READ(u,f[,END=S 4] [,ERR=S2]) )[list] 


In the preceding format, u is one of the following: an unsigned INTEGER*2 constant which 
is the logical unit number of the input device, an INTEGER*2 variable name whose value is 
the logical unit number of the input device, or a CHARACTER variable or array name which 
is the internal input file; f is the format reference to be used; END=S, defines a statement 
number to receive control when an end-of-file is encountered; ERR=S2 defines a statement 
number to receive control when an error is detected; and list is a list of variable names, 
separated by commas, which are assigned the values of the input data. The parameters 
END=S and ERR«S2 are optional and may be in reverse order. The list indicates the number 
of items of data to be placed in memory and the order in which they are read, Any number 
of items may appear in a list. 
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For example, consider the following logical record: 
11-337 2 400 


The following statements can be used in a program to input the information from the above 
logical record. 


READ(5S, 3)N, I, 5,4 
3 FORMAT (I2, 1X. 14, 2X, 11, 3X, 13) 


These statements cause the variable N to be assigned the value 11, |! the value -337, S the 
value 2, and K the value 400. 


Each time the READ statement is executed, the variables obtain new values depending upon 
what information is in the next record read. 


If an input record contains more values than there are items in a list, only enough values 

to fill the list are read. The remaining values are ignored, If a list requires more values than 

are on a record, the next sequential record or records are read until the list has been satisfied. 
The READ (u,f) form may be used with a FORMAT statement to read alphanumeric data and 
modify an existing H-type field (described in paragraph 5.2.3.3). The information read is 
inserted into the memory location of the H-type field. The amount of data to be read in depends 
on the size of the field. 


EXAMPLE 


é FORMAT(C21LHTHIS IS A REPLACEMENT > 


READ(5S, 6) 


These statements cause the 21 characters in the H-type alphanumeric field 
THIS IS A REPLACEMENT 


to be replaced by the first 21 characters from the next record of the file connected to logical 
unit 5. 
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5.3.3.3 Formatted WRITE Statement 


The WRITE statement causes the computer to transfer information from memory to an 
output device, or to an internal file in memory. 


FORMAT 

WRITE(u, fl, ERR=So] )[list] 
In the preceding format, u is either an unsigned INTEGER*2 constant or variable whose value 
is the logical unit number of the output device, or u is a CHARACTER variable or array name 
which is the internal output file; f is a format specifier; ERR=S> defines a statement number 
to receive control when an error is detected; and list is a list of variable names which are 
assigned the values of the output data. The parameters ERR=Sz and list are optional. 
The WRITE(u,f) list form is used to write the data currently stored in memory, corresponding 
to the variable names in the list, onto the file designated by u, using a form specified by the 
FORMAT statement f. lw 


The WRITE(u,f) form is a means of writing alphanumeric data, The FORMAT statement specifies 
the exact data to be written, so the need for an 1/0 list is eliminated, 


EXAMPLE 
WRITE(6, 10) 


ic FORMAT ( 19HTHIS IS A QUOTATION) 


This could also be written 


WRITE (6, 10} 
10 FORMAT( ‘THIS IS & QUOTATION’) 


Two apostrophes together represent one in the format. i 
EXAMPLE 


' 
FORMAT(‘DON"T") 
| 


would be interpreted as (SHDON'T). 
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5.3.3.4 Unformatted READ and WRITE Statements 


FORMATS 
READ(ul, END=S 4] [, ERR=Sq] list 


WRITE(ul, ERR=Sz] Dist 


In the preceding formats, u is an unsigned INTEGER*2 constant, variable, or expression 
assigned the value of the logical unit number used; END=S, defines a statement number to 
receive control when an end-of-file is encountered; ERR=S> defines a statement number 
to receive control when an error is detected; and list is a list of variable names separated 
by commas. Variable names declared to be of type CHARACTER may not appear in an un- 
formatted I/O list. The parameters END=S, and ERR=Sp are optional and may be in reverse 
order. 


The READ(u)list form reads a record of binary information, with no conversion, into memory 
from logical unit u. No FORMAT statement is required. The number of items in the list deter- 
mines the amount of data read. If the record contains more values than are listed, the unread 
items are skipped. The total length of the list of variable names must not be longer than the 
record. 


Omitting the list from the READ statement causes the record on unit u to be skipped and 
no information to be transmitted from the input device. 


The WRITE(u)list form is used to write binary information, with no data conversion, on logical 
unit u. 


5.3.4 File Positioning Statements 
There are three types of mass storage file input/output statements: REWIND, BACKSPACE, 


and END FILE, Use of these statements is limited to magnetic tape, tape cassette, and disk 
applications. 
5.3.4.1 REWIND Statement 
FORMAT 
REWIND u 


In the preceding statement, u is an unsigned INTEGER*2 constant or variable whose value 
is the logical unit number of the mass storage file. 
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This statement directs the designated |/O unit to reposition to the first record of the 
first file. 
EXAMPLES 
REWIND 7 


REWIND KMINO 


5.3.4.2 BACKSPACE Statement 
FORMAT 
BACKSPACE u 
In the preceding format, u is an unsigned INTEGER*2 constant or variable whose value is 
the logical unit number of the mass storage file. 
EXAMPLES 
BACKSPACE 7 


BACKSPACE N 


5.3.4.3 END FILE Statement 


FORMAT 
END FILE u 
In this statement u is an unsigned INTEGER*2 constant or variable whose value is the logical 
unit number of the mass storage file. 
This statement writes an end-of-file mark at the current file position. 
EXAMPLES 
END FILE 4 


END FILE K 
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5.4 Direct Access INPUT/OUTPUT 


Direct access I/O statements allow the programmer to access records within a file in a 
random, rather than sequential, manner. For direct access, a file is viewed as a collection 
of n records, and each record is assigned a unique record number in the range 0 to n-1. 
The programmer must specify in READ, WRITE, and FIND statements not only the unit 
number, as for sequential |/O, but the number of the first record to be read, written, or 
found. There are four direct access I/O statements: READ, WRITE, OPEN, and FIND. 


5.4.1 OPEN Statement 


The OPEN statement may be used to connect an existing file to a unit number, or to create 
a file and connect it to a unit number. The OPEN statement is used for declaring direct 
access files. It is the functional equivalent of the DEFINE FILE statement supported under 
previous releases of FORTRAN. The OPEN statement must logically precede all |/O state- 
ments referring to the defined unit number. Only direct access I/O statements may refer to 
units defined by OPEN statements. The advantage of the OPEN statement is that it provides 
mnemonic keywords for specifying its parameters. 


Additional information describing file creation and physical records is found in Appendix B 
of the FORTRAN~-78 User's Guide for the TAXO Operating Systems. 


FORMAT 


OPEN (u,ACCESS=DIRECT,RECL=1,FORM=f, ASSOC V=i, NUMREC=n) 


All the parameters in the parentheses except u are optional and may appear in any order. 
In the preceding format: 


° 


u is an INTEGER*2 variable or constant that specifies the I/O unit designation. 
The unit number must be the first member of olist. 

° ACCESS=DIRECT specifies the access method, The only access method permitted 

is DIRECT, which is also the default method, 


RECL=l designates the INTEGER*2 variable or constant that specifies the record 
length. The default value is 288 bytes. 


FORM=f specifies the format control for |/O operations. The default value is 
UNKNOWN. If the format type is UNKNOWN or FORMATTED, the record length 
is in bytes, If the format type is UNFORMATTED, the record length is in words, 
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° ASSOCVsi designates the INTEGER*2 variable that assigns the associate 


variable. At the conclusion of each read or write operation, the associate variable 


is set to the record number of the record that immediately follows the last record 


transmitted. At the completion of a find operation, the associate variable is set 
to the record number of the record found, The default is a null variable not 
accessible by the user. 


° NUMRECsn designates the INTEGER*2 variable or constant that defines the number 


of records. The default is 6. 


EXAMPLE 


OPEN (5,ACCESS=DIRECT,RECL=80,NUMREC=500,ASSOCVsID) 


In this example assume that the file associated with unit 5 does not exist. The OPEN state- 


ment will create and then open a relative record file with 500 records 80 bytes long. The 
associate variable is ID and the format type is UNKNOWN, 
5.4.2 READ Statement 
FORMATS 
READ (u'nl, f] Ty END=S,] Ly ERR=So] Jlist 


READ (ul, f] ,REC=n[, END=S 4] Ly ERR=S9] list 


In the preceding statements 
° uis an |/O unit designation, 
° nis an integer expression representing the relative record number. 


° f is an optional format reference. 


© RECsn is an alternate form for specifing the record number.If REC=n is used with 


END=S, and ERR=So, the options may be in any order. 


© END=S, and ERR=S> are optional and their order may be reversed, They indicate 


exit destinations in the event of an end-of-file or an error, respectively. 


° list is a list of variable names separated by commas. 


~ 
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EXAMPLE 


DIMENSION K(6) 
OPEN( 4, ACCESS=DIRECT, NUMREC=500, RECL=680, 
* FORM=FORMATTED, ASSOCV=IF4) 

106 FORMAT (/5(7X, I3)} 

100i READ(4, 100, REC=6)(K(I}, I=1, 6) 

iG02 READ(4‘IF4-3, 100}K 


The READ statement 1001 reads records 7 and 8 from the file associated with unit 4. The 
READ statement specifies record 6, but the format statement contains a leading slash 
which causes a skip to the following record. Two records must be read in order to satisfy 
the I/O list of six elements in the array K. At the completion of the first READ, the 
associated variable IF4 is set to the value 9, the record following the last one read, Thus, 
the second READ, statement 1002, performs exactly as the first READ because the relative 
record expression, IF 4-3, is equal to 6. 
5.4.3 WRITE Statement 
FORMATS 

WRITE (u'nl, f] L, ERR=So] Dlist 


WRITE (ul, f] ,REC=n[, ERR=So )list 


In the preceding statements 
° uis an 1/O unit designation. 
° nis an integer expression representing the relative record number. 
° f is an optional format reference. 


° RECen is an alternate form for specifying the record number. If REC=n is used with 
ERR=S>, the options may be reversed. 


° ERR&So is optional and indicates an exit destination in the event of an error. 


EXAMPLE 


OPEN(S, ACCESS=DIRECT, NUMREC=40, RECL=20, 
* FORM=UNFORMATTED, ASSOCV=IFS> 
I=24 
8& WRITE(S, REC=1+134, B,C 
99 WRITE(S ‘IFS? XY 
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The first WRITE statement, number 88, writes record number 25 and the second WRITE, 
number 99, writes record 26. 


5.4.4 FIND Statement 
The FIND statement locates the next input record while the present record is being processed, 
This increases the execution speed of the program. FIND does not transfer any data to or 
from the file. There is no advantage to preceding a WRITE with a FIND. 
FORMATS 

FIND (u'n) 


FIND (u,REC=n) 


In the preceding statements 
° uis an integer I/O unit designation. 
° nis an integer expression representing the relative record number. 


° RECen is an alternate form for specifing the record number, 


EXAMPLE 


FIND(3‘1526) 


READ(3 1526) x 


5.5 Console DISPLAY INPUT/OUTPUT ¥ 


The ACCEPT and DISPLAY statements allow the user to transfer data to and from the screen 
of video display terminal (VDT) 1/0 devices. 

5.5.1 ACCEPT Statement 

FORMAT 


ACCEPT(u,f,LINE=n 1,POSITION=n2,ERASE,PROMPT, 
GRAPHICS,ECHO,ERRes)list 
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All the parameters in the parentheses, with the exception of u and f, are optional and 
may appear in any order. 


In the above statement 


u is an INTEGER*2 variable or constant that designates the |/O unit number 
assigned to the VDT. 


f is a format reference. 


LINE=n1 designates the INTEGER*2 variable or constant n1 as the line number 
on the VDT screen, in the range 1 to the number of lines on the screen, where 
the first record is accepted. If the value is greater than the number of lines on 
the VDT screen, it is adjusted modulo to the maximum line number, If the value 

is zero or the LINE phrase is not present, the data is accepted from the next line 
at the specified cursor position. However, If the value specified in the POSITION 
phrase is also zero, the data is accepted from the current cursor position with- 
out incrementing the line number. If incrementing to the next line generates a line 
number greater than the maximum number of lines on the screen, an automatic 
erase occurs with reset to line one, 


POSITION=n2 designates the INTEGER*2 variable or constant n2 as the position 
number on the VDT screen, in the range 1 to the number of characters per line, 
where the first record is accepted. If the value is greater than the maximum number 
of characters within a line on the screen, it is adjusted modulo to the maximum 
number of characters. If the value is zero, the data is accepted from the current 
cursor position. If the POSITION phrase is not present, POSITION is set equal to 1. 


ERASE causes the screen to be cleared before input. 


PROMPT causes the input field on the screen to be filled with asterisks before 
input. 


NOTE 


When using PROMPT for alphanumeric input, the system will not allow 
asterisks to be keyed in, 


° 


GRAPHICS allows the control characters >00 to >1F to be input as graphic characters, 
The graphics character value must be accepted into a character variable as shown 
in the following. 


CHARACTER#& GCHARS 
c ACCEPT 8 GRAPHICS CHARACTERS 
ACCEPT (&, ‘(A} ‘, GRAPHICS, PROMPT, ECHO) GCHARS 
stop 
END 
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Refer to the TAXO Manual SCI Commands ect. for further information, Note that 
control keys such as RETURN generate graphic characters in this mode. 


ECHO causes the contents of an item to be displayed on the VDT screen. Conversion 
and justification occur prior to display. If the ECHO phrase is not specified, 
the original input data remains in the field. 


ERR = s indicates a statement number to which control is transferred when an 
error is detected. 


° list is any normal I/O list of variables, arrays, and array elements. 
EXAMPLE 


ACCEPT (10, 99, LINE=3, POSITION=10, PROMPT) ICOUNT 
99 FORMAT(IS) 


This ACCEPT statement would display five asterisks on line 3 in positions 10 through 14, 
position the cursor at the first asterisk, and wait for a five-digit integer to be input. Input 

is terminated by filling the field or pressing the line termination key. When input is terminated, 
the value is transferred to the variable ICOUNT. 


NOTE 


The item is accepted under normal FORTRAN input conversion rules. If the 
number of characters input is less than the field length, the data is right- 
justified within the field. 


To input alphanumeric data, use character variables or arrays, To insure that the input field 
width agrees with the alphanumeric character capacity of the variable type, the following 
formats should be used. 


Variable Type Format Specification 
CHARACTER or CHARACTER*1 Al 
CHARACTER*2 A2 
CHARACTER*3 A3 


CHARACTER#*n An 
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Any normal I/O list can be accepted. Each element in the list is accepted separately. 

The data is transferred to the variables in accordance with the associated format statement, 
One or more fields are contained in each record, Each record is limited to the width of the 
VDT screen. The first record begins at the specified line and position and continues to the 
end of that line. Each succeeding record causes the line number to be incremented and the 
initial position to be reset to the position specified for the first line. 


5.5.2 DISPLAY Statement 
FORMAT 


DISPLAY (u,f,LINE=n 1,POSITIONsn2,GRAPHICS, 
INTENSITY sHIGH/LOW,ERASE,ERRe=s)list 


All the parameters in the parentheses with the exception of u and f are optional and may 
appear in any order. 


In the preceding statement 


° uis an INTEGER*2 variable or constant that designates the I/O unit number 
assigned to the VDT. 


° fis a format reference. 


° LINE=n1 designates the INTEGER*2 variable or constant n1 as the line number 
on the VDT screen, in the range 1 to the number of lines on the screen, where 
the first record is displayed. If the value is greater than the number of lines on 
the VDT screen, it is adjusted modulo to the maximum line number. If the value 
is zero or the LINE phrase is not present, the data is displayed at the next line 
at the specified cursor position. However, if the value specified in the POSITION 
phrase is also zero, the data is displayed at the current cursor position without 
incrementing the line number. If incrementing to the next line generates a line 
number greater than the maximum number of lines on the screen, an automatic 
erase occurs with reset to line one. 


° POSITION=n2 designates the INTEGER*2 variable or constant n2 as the position 
number on the VDT screen, in the range 1 to the number of characters per line, 


where the first record is displayed. If the value is greater than the maximum number 


of characters within a line on the screen, it is adjusted modulo to the maximum 
number of characters. If the value is zero, the data is displayed at the current 
cusor position. If the POSITION phrase is not present, POSITION is set equal to 1. 
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° GRAPHICS allows the control characters >00 to >1F to be output as graphic 
characters. Since a graphics character value cannot be equated to a character 
variable, an EQUIVALENCE statement must be used to assign the graphics character 
value to a character variable as shown in the following. 


CHARACTER#40 LINE 
INTEGER#2 ILINE(20} 
EGUIVALENCE (LINE, ILINE} 
DATA ILINE/20#>16/ 

c WRITE GRAPHIC LINE TO DISPLAY 
DISPLAY (6. ‘(A)’, GRAPHICS}SLINE 
STOP 
END 


Refer to the TAXO Manual SCI Commands ect. for further information. Note that 
control keys such as RETURN generate graphic characters in this mode, 


° INTENSITY=HIGH/LOW allows the data to be displayed in either high or low intensity 
The default is low intensity. For futher information refer to the TAXO Manual SCI 
Commands ect.. 


° ERASE causes the screen to be cleared before output. 


° ERR =s indicates a statement number to which control is transferred when an 
error is detected. 


° list is any normal I/O list of variables, arrays, and array elements, 


To output alphanumeric data, use character variables or arrays. Refer to the ACCEPT statement 
for information on alphanumeric field widths. 


Any normal |/O list can be displayed. The data is transferred from the variables in accordance 
with the associated format statement. One or more fields are contained in each record, The 
first record begins at the specified line and position and continues to the end of that line. 

Each succeeding record causes the line number to be incremented and the initial position 

to be reset to that specified for the first line, i) 


5.6 FORTRAN Unit Numbers 


Input and output unit numbers used in FORTRAN programs can be any INTEGER*2 value 

form1 to 99, These unit numbers must be assigned to actual units before the program is exe- 
cuted, The FORTRAN-~78 User's Guide for the TAXO Operating Systems contains LUNO assign- 
ment under TAXO. 
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When working with FORTRAN I/O operations at a particular unit number, only certain operations 
can follow a given operation. If one logical unit number is used to perform a WRITE, READ 
sequence (i.e., both input and output), a REWIND must be performed after the WRITE command, 
Refer to the FORTRAN-78 User's Guide for the TAXO Operating Systems to validate TAXO 

1/O operations, 
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6 EUNCTIONS AND SUBPROGRAMS 


6.1 General 


With the FORTRAN language it is possible to specify relatively complicated operations, such 

as the determination of the sine of an angle or the printing of variable multiline header 
information, in one program statement. This is accomplished by first programming the operation 
to be performed in a subprogram, and then accessing the code for this operation explicitly 

by a CALL statement or implicitly by a function reference. Several methods and requirements 
exist for using the available FORTRAN facilities; they are discussed in detail in this section, 


6.2 Subprograms 


The following paragraphs define several terms needed to explain arithmetic functions and 
subprograms, and describe the purpose of dummy identifiers. 


FORTRAN subprograms may be internal or external, An internal subprogram is defined within 
the program that calls it and is defined within a single statement, the function definition state- 
ment. Internal subprograms are defined and may be used only within the program containing 
the definition. External subprograms are defined separately from (external to) the program 
that calls them and are complete programs conforming to all the rules of FORTRAN programs, 
They are compiled independently. 


Two types of external subprograms may be defined: function subprograms and subroutine 
subprograms. The use of the declarations FUNCTION and SUBROUTINE in the definition of 
these subprograms is described in paragraphs 6,4 and 6.5. 


Intrinsic functions are external subprograms which are predefined in the FORTRAN language. 
The intrinsic function identifiers, definitions, and types are given in Section 7. These definitions 
are overridden by using the function identifier in any context other than a function reference. 


Any subprogram, internal or external, may call other subprograms: and recursion is allowed. 
If a subprogram calls itself, either directly or indirectly, recursion exists (refer to paragraph 6.9). 


Subprogram definition statements declare certain identifiers to be dummies representing the 
arguments of the subprogram. They are used as ordinary identifiers within the subprogram 
definition and indicate what sort of arguments may appear and how the arguments are used. 
The dummy identifiers are replaced by the actual arguments when the subprogram is executed, 


Dummy identifiers may not appear in COMMON, EQUIVALENCE, or DATA statements. 
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When a subprogram is called, the number of dummy parameters must be the same as the 
number of actual parameters, and the parameters should agree in their order, type, and 
length. For example, if an actual parameter is a 16-bit integer constant, then the correspon- 
ding dummy parameter should be of INTEGER*2 type. 


The number of parameters is checked and must match at compile time and at run time, Other 
attributes are not checked. If the actual and dummy parameters disagree in type, the sub- 
program must be aware of the disagreement. 


If a dummy parameter is an array name, the corresponding actual parameter may be either 
an array name or an array element. 


lf a dummy parameter is assigned a value in the subprogram, the corresponding actual para- 
meter should be a simple variable, array element, or array name. A constant or expression 
should not be used as an actual parameter if the corresponding dummy parameter may be 
assigned a value. 


6.3 Statement Function Definition Statement 


The statement function definition statement defines an internal subprogram. The single state- 
ment contains the entire definition of the function. All statement function definitions must 
precede the first executable (nondeclarative) statement of the program. 


FORMAT 


identifier(Lidentifier 4,identifiero,...] )=expression 


To the left of the equal sign, an identifier followed by zero or more identifiers in parentheses 
appear. The first identifier is the function; the identifiers in parentheses represent arguments 
of the function and are known as dummy identifiers. To the right of the equal sign, the state- 
ment contains an expression whose value is the value of the function. 


The function identifier is the name of the subprogram being defined, A function is single-valued 
and may have any number of arguments, The data type of the function is determined by 
the type of the function identifier. 


The arguments of the function are named by the dummy identifiers, which appear in parenthe- 
ses after the function name. The following rules apply to dummy identifiers. 


° They must be unique only within the individual function definition statement. 


© They may be identical to identifiers of the same type appearing elsewhere in the 
program. 


° They must agree in order, number, type, and length with the actual arguments 
provided at execution time. 
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The portion of the function definition statement to the right of the equal sign must be an 
expression. The expression defines the value of the function. The occurrence of a dummy 
identifier determines how the identifier is handled in the defining expression, Identifiers that 
do not represent arguments are treated as ordinary variables, The following rules apply to 
the defining expression. 


° Dummy identifiers may appear only as scalar variables in the defining expression, 
not as subscripted variables. 


° The defining expression may include references to external functions or other 
previously defined internal functions. 


EXAMPLES 


SSQR CK} =K# (CK+1 2 # (QHK4+1 976 

NOR(T, S}=. NOT. (T. OR-S) 

ACOSH(X 3 =(CEXP(X/A) EXP (-X/A) 972. 
ROOT1 ( }=(—-B+SGRT (B##2—-4, #A#C})/(2. #A) 


In the second example above, the data type of the function NOR must have been previously 
declared to be logical. In the third example above, X is a dummy identifier and A is an ordi- 
nary identifer. At run time the function is evaluated using the current value of the quantity 
represented by A. In the last example above, there are no dummy identifiers (arguments), and 
the statement function only references ordinary program variables. 


6.4 Eunction Subprogram 


A function subprogram is a function which is single-valued and is referenced as a basic element 
in an expression. The subprogram is compiled separately from the main program, yet there may 
be an interchange of data between it and the main program, Its variable names are completely 
independent of the main program and other subprograms. A function subprogram is called 
implicitly by use of a function reference in the main program. 


A function subprogram begins with a FUNCTION declaration and returns contro! to the calling 
program by means of one or more RETURN statements, The function subprogram must terminate 
with an END statement. 


FORMAT 

[type] FUNCTION identifier(Lidentifier 4,identifiers,...] ) 
This statement declares that the program which follows is a function subprogram. The first 
identifier is the name of the function being defined. This identifier must appear as a scalar 


variable and be assigned a value during execution of the subprogram., This value is the function 
value. 
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EXAMPLE 


FUNCTION ROOTI(CA. B,C) 

ROOT1I=(—-B+SGRT (B##2-4. O#A#C))/ (2. OFA) 
RETURN 

END 


Identifiers appearing in the list enclosed in parentheses are dummy identifiers representing 

the function arguments. They must agree in number, order, type, and length with the actual 
arguments given at run time. A function subprogram may have any number of arguments, 
including arrays, within the limitations of available storage area, Function subprogram argu- 
ments may be expressions, array names, or subprogram names, Dummy identifiers may appear 
in the subprogram as scalar identifiers, array identifiers, or subprogram identifiers. 


Dummy identifiers which represent array names must be dimensioned within the subprogram 
by a DIMENSION or type statement. Dimensions given as constants must equal the dimensions 
of the actual arrays given at run time. In a DIMENSION or type statement, dummy identifiers 
may be used to specify variable dimensions for array name arguments, For example, in the 
statement sequence 


FUNCTION TABLE(A, M,N, B, X,Y) 


DIMENSION A(M,N},B(10),C(50) 


the dimensions of array A are specified by the dummies M and N, and the dimension of array 
B is given as a constant, The values given for M and N at run time must be the actual dimensions 
of array A, ° ’ 


Dummy dimensions may be given only for dummy arrays, In the preceding example, the array 
C must be given absolute dimensions since C is not a dummy identifier. 


The type of the function is the type of the identifier which names the function, Alternately, the 
function may be explicitly typed in the FUNCTION statement by replacing the word FUNCTION 
with one of the following identifiers. 


INTEGER FUNCTION 
INTEGER*2 FUNCTION 
INTEGER*4 FUNCTION 

FIXED (scale) FUNCTION 

REAL FUNCTION 

REAL*4 FUNCTION 

REAL*8 FUNCTION 

DOUBLE PRECISION FUNCTION 
COMPLEX FUNCTION 

LOGICAL FUNCTION 


Note that functions of type CHARACTER are not permitted. 
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The symbolic name of the function may not appear in any nonexecutable statement except 
as the name in a FUNCTION statement. 
For example, the statement 

COMPLEX FUNCTION HPRIME(S,N) 
defines a complex function HPRIME. 
EXAMPLES 

FUNCTION MAY (RANGE,XP,YP,ZP) 

REAL FUNCTION COT(ARG) 
The first example defines an INTEGER*2 function subprogram MAY, with four dummy argu- 
ments. The second example explicitly defines the subprogram COT to be REAL. 
6.5 Function Reference 
A mathematical function is a quantity or group of quantities that are defined in terms of 
relations between or operations upon some other quantity or quantities over a specified 
range of values. The value of the function depends on one or more independent variables, 
These independent variables are called arguments. 
EXAMPLE 

f(x) = x24 2x +1 
The function value f(x) is defined in terms of the single argument x. For a given argument 
value such as x = 3, there is an associated function value, in this case f(x) = 16. 
In a program, an arithmetic function is a subprogram that acts upon a number of arguments 
and produces a single numberical quantity as the function value. Function references are 
denoted by the identifier which names the function, followed by an argument list enclosed 
in parentheses. If the user wishes to code a function subprogram with no arguments, the 
argument list consists of a left parenthesis followed by a right parenthesis. 
FORMATS 


identifier (argument 4, argument>,...,argument,, ) 


identifier ( ) 
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An argument may be an expression, an array identifier, or a subprogram identifier. 


A function reference represents a quantity, namely the function value, and acts as a basic 
element. The type of the function value is given by the type of the identifier which names 
the function. The type of the function is independent of the types of its arguments. 


EXAMPLES 
X = COS(THETA) 
| = IZETA(S+SQRT(S)) 


A = FUNC( ) 


NOTE . 


In these examples, a real value X is returned for the function COS, 
an integer value | is returned for IZETA and a real value A is returned 
for FUNC, in accordance with the rules for implicit declarations of 
variables, 


References to logical functions are written in the same way as references to arithmetic 
functions. The identifier used to name the function must be a logical identifier. The function 
arguments may be expressions (logical or arithmetic), array identifiers, or subprogram iden- 
tifiers. 


6.6 Subroutine Subprogram 


Another type of subprogram that is compiled by itself is the subroutine subprogram. A sub- 
routine subprogram is not a function; it may be multivalued and can be referred to only by 
a CALL statement. A subroutine subprogram begins with a SUBROUTINE declaration and 
returns control to the calling program by means of one or more RETURN statements. There 
are two possible forms for the SUBROUTINE declaration, 


FORMATS 

SUBROUTINE identifier 

SUBROUTINE identifier(Lidentifier 4, identifiero,...] ) 
This statement declares that the program that follows is a subroutine subprogram. The first 
identifier is the subroutine name. The identifiers in the list enclosed in parentheses are 
dummy identifiers representing the arguments of the subprogram. These identifiers must 


agree in number, order, type, and length with the actua! arguments given to the subprogram 
at run time. 
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Subroutine subprograms may have expressions, array names, and subprogram names as 


arguments. The dummy identifiers may appear as scalar, array, or Subprogram identifiers 
within the subprogram. 


Dummy identifiers that represent array names must be dimensioned within the subprogram 
by a DIMENSION or type statement. As in the case of a function subprogram, either constants 
or dummy identifiers may be used to specify dimensions in a DIMENSION or type statement. 


A subroutine subprogram may use one or more of its dummy identifiers to represent results, 
The subprogram name is not used for return of results. 


A subroutine subprogram need not have any arguments. 
EXAMPLES 
SUBROUTINE EXIT — 
SUBROUTINE FACTOR (COEF,N,ROOTS) 
SUBROUTINE HEADING ( ) 


SUBROUTINE RESIDU (NUM,D,DEN,M,RES) 


The example in figure 6-1 illustrates the principle features of a subroutine subprogram. The 
arguments of subprogram AVGN include values that are inputs to the subroutine and results 
calculated by the subroutine for output to the main program. The inputs are an array name 
and the number of values to be averaged. The subprogram calculates an average value and 
the number of values that are greater than or equal to Zero. 


To call this subprogram, the main program might contain a statement such as 


SUBROUTINE AVGN(X, J, AVG, IZ) 
SUBROUTINE COMPUTES AVERAGE OF FIRST J ELEMENTS 
AND THE NUMBER OF ELEMENTS THAT ARE GREATER THAN 
OR EGUAL TO ZERO. 
DIMENSION X(100) 
AVG=O. 0 
IZ=0 
IF (J. GT. 106) RETURN 
DO 10 I=1,J 
AVG=AVG+X (I) 
IFCX(1I). GE. 0. O)IZ=I1Z+i 
16 CONTINUE 
XNwJ 
AVGRBAVG/ XN 
RETURN 
END 


aan 


Figure 6-1. Example of a Subroutine 
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This statement calls the subprogram to calculate the average of the first 30 values of an 
array named ALTA. The computed average is stored as the value of YMEAN, and the 

value IZ computed by the subprogram is stored as the value of NPOS. If the main program 
later requires the average of the first 52 values of an array named ARRAY, with the com- 
puted average to be stored in FPT3 and the number of values greater than or equal to zero 
stored in K, it would include this statement: 


CALL AVGN (ARRAY,52,FPT3,K) 


The calling program must provide values for the subprogram's input variables and must 
contain defined variables to store the values calculated by the subprogram. 


The subprogram need not have arguments, as illustrated in the following example. 


EXAMPLE 


SUBROUTINE WARN 
WRITES. 15) 
15 FORMAT(’ & DEFAULT HAS OCCURED ’?> 
RETURN 
END 


6.7 Specification Subprograms For Data Initialization 


A specification subprogram, unlike the function and subroutine subprograms, may not be 
referenced by any other program. This is indicated in the structure of the specification sub- 
program, i.e., it contains no executable statements. The purpose of this type of subprogram 
is to establish initial values of variables in the labeled common areas, 


A specification subprogram begins with a BLOCK DATA declaration; therefore, it is some- 
times referred to as a block data subprogram. Following the opening declaration, there may 
be DIMENSION, type, COMMON, and DATA deciarations as required to specify the values 
to be assigned to the variables. 

The data specification statements DATA and BLOCK DATA are used to specify initial values 
for variables. These values are compiled into the object program. They become the values 
assumed by the variables when execution begins. 

6.7.1 DATA Statement 

FORMAT 


DATA Vv /d /V2/G2/ yore 
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In the preceding format, v is a variable list and d is a data list. 


The variable lists in a DATA statement consist of scalar variables, array names, or array elements 
separated by commas. 


Variables in common may appear in the lists only if the DATA statement occurs in a block 
data subprogram. 


DATA statements may also be used in main programs and in function and subroutine subprograms 
to set values for noncommon variables. When used in this way, the DATA statements must be 
placed before arithmetic function definition statements, if there are any, and must precede 

all executable (nondeclarative) statements. 


The storage area that will be occupied by the initialized data items of each data list must 
match the storage length of the initialization. For example, simple integer data requires one 
word and real number data requires two words. The storage requirements of the seven data 
types are listed in the discussion of types of data in Section 2. 


The following considerations apply to data items in a DATA statement, 


° Data items may be numeric constants or alphanumeric strings. For example 
y 9) 


DATA ALPHA,BETA/5., 16.E-2/ 


specifies the value 5.0 for ALPHA and the value 0.16 for BETA. The type of the 
variable determines the data type of the stored constant. 


Character strings can only be used to initialize items of type CHARACTER, The 

length of the character constant cannot exceed the length of the item it is initializing. 
If the length of the character constant is less than the length of the character variable 
it initializes, blank characters are padded on the end, 


CHARACTER ANSW*4,NAME* 10 
DATA ANSW/'NO'/,NAME/'JOHN'/ 


Character arrays whose elements have odd lengths cannot have their even-numbered 
elements initialized because they begin on odd-byte boundaries, Only individual, 
odd-numbered elements of such arrays may be initialized by the DATA statement. 


CHARACTER*3 TABLE(5),LIST(2) 


DATA TABLE(2),TABLE(4)/'***",'$$$'/,LIST/2*"BAD'/ 
DATA TABLE(1), TABLE(3), TABLE(5)/3*'AOK'/ 


The first data statement above contains illegal initializations that will result in a 
syntax error and the second statement contains legal initializations. 


Any data item may be preceded by an integer and an asterisk. The integer indicates 
the number of times the item is to be repeated. 


DATA A(1),A(2),A(3),A(4) ,A(5)/6 1E2,4% 32E 1/ 
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The preceding statement specifies five values for the array A: the value 6100 for A(1) and 
the value 320 for A(2) through A(5). 

© When an unsubscripted array name is included in the variable list, all elements 
of the array are initialized to the indicated values. The notation is equivalent to 
specifying each element of the array individually. 


DIMENSION A(5) 
DATA A/61E2,4*32E1/ 


is equivalent to 


DATA A(1),A(2),A(3),A(4),A(5)/6 1E2,4*32E 1/ 


6.7.2 BLOCK DATA Statement 
FORMAT 
BLOCK DATA [name] 


This statement declares that the program which follows is a data specification subprogram. 
Data specification for variables in common blocks requires the use of a block data subprogram. 


The first statement of the subprogram must be the BLOCK DATA statement. The keywords 
BLOCK DATA may be followed by an optional subprogram name. If none is given, the default 
is $BLOCK. The user should not have more than one unnamed block data subprogram linked 
in one task. The subprogram may contain only the declarative statements associated with 
the data being defined. 


EXAMPLE 


BLOCK DATA SPEC 

COMMON /R/X, Y/C/Z, WV 

DIMENSION Y(2} 

COMPLEX Z 

DOUBLE PRECISION X 

DATA Y(1),¥(23, Y(3)/1E-1, 2e3E2/ 

DATA X, Z/11. 877D0, (-1. 41421,1. 41421)/ 
END 


Data may be entered into more than one common block in one subprogram, A blank common 
cannot be initialized. Any common block mentioned must be listed sufficiently to define those 
variables being used. However, the first reference to a common block must be listed in full. 
In the example above, W and V are listed in block C although no data values are defined 

for them. 
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6.8 EXTERNAL Statement 


FORMAT 
EXTERNAL identifier ,identifier,...,identifier 


This statement declares that the listed identifiers are external subprograms or subprograms 

that are passed as arguments to other subprograms, Any subprogram name given as an argument 
to another subprogram must appear in an EXTERNAL declaration in the calling program. If an 
intrinsic function appears in an EXTERNAL statement, it is not available for use in the current 
program. 


EXAMPLE 


EXTERNAL SUB1, SUB2 

CALL EVAL(SUB1, 1. 5, ANSWER? 
CALL EVAL(SUB2, . 87, ANSWER) 
END 


SUBROUTINE EVAL (SUB. ARG, ANSWER > 


CALL SUB (ARG, ANSWER) 


RETURN 
END 


6.9 INTRINSIC Statement 


FORMAT 


INTRINSIC identifier,identifier,...,identifier 
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This statement declares that the listed identifiers are the names of intrinsic functions that 
are to be passed as parameters to another program module. The names of intrinsic fuctions 
for type conversion, lexical relationship, and for largest or smallest value must not be used 
as parameters. 


EXAMPLE 


INTRINSIC SIN, COS, ALOG 


6.10 RECURSIVE Statement 
FORMAT 
RECURSIVE [identifier,identifier,...,identifier] 


The RECURSIVE statement is used to declare subprogram modules to be recursive. A recursive 
subprogram is one that calls itself either directly or indirectly. In the case of direct recursion 
where a subprogram contains a statement which references that subprogram, the compiler 

can automatically determine that the recursion exists, and the RECURSIVE statement is not ad 
needed in that program module. The RECURSIVE statement is needed to declare that subprogram 
to be recursive in any program module which references it. 


EXAMPLE 


PROGRAM MAIN 
RECURSIVE IFAC 


I=IFAC (J) 


STOP _ 
END 


INTEGER FUNCTION IFAC(N) 
IF(N. LT. 2) THEN 

IFAC=i 
ELSE 

IFAC=N# IFAC (N-i > 
ENDIF 
RETURN 
END 


In the example above, the main program calls recursive function IFAC, which calculates the 
factorial of its input argument. The RECURSIVE statement is needed in the main program to 
correctly call the recursive function. 
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For indirect recursion, the RECURSIVE statement without an identifier list must be used in 
every subprogram involved in the indirect recursion and the RECURSIVE statement with an 
identifier list must be used in any program module which calls one of the subprograms involved 
in the indirect recursion in order to declare them to be recursive. 


EXAMPLE 
— =O PROGRAM MAIN 
RECURSIVE A,C 
CALL A 
CALL C 


STOP. 
END 


SUBROUTINE A 
RECURSIVE 
RECURSIVE B 
CALL B 

END _ 
SUBROUTINE B 


RECURSIVE 
RECURSIVE C 


CALL C 

END 
SUBROUTINE C 
| RECURSIVE 
RECURSIVE A 


CALL A 


END 
When using recursive subprograms, the user must take great care to guarantee that an infinite 
calling loop does not exist. If this condition does occur, the program will increase in size until 
there is no more system memory available for the dynamic data area stack, At that point 
the task will terminate giving traceback information describing all the recursive calls which 
occurred, 


NOTE. 
Recursive subprograms depend on the dynamic memory allocation services offered when 


running with a standard operating system. For this reason they cannot be used in a 
stand-alone environment. 
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FORTRAN-78 RUN-TIME LIBRARY 


7.1 Library Functions 


The FORTRAN system supplies a library of standard functions that may be referenced from 
any program. These functions are called intrinsic functions and are listed and defined in table 
7-1. The intrinsic functions are referenced in a program in the same manner as normal user- 
supplied functions. The intrinsic functions are known by the compiler and may be treated 

in nonstandard ways (such as inline function expansion). The programmer may replace the 
intrinsic functions by including the name of the new functions in an EXTERNAL statement 

in all calling programs. If an intrinsic function is passed as an argument to a subprogram and 
is not declared as an external function, the function name must be included in an INTRINSIC 
statement. 


Table 7-1. FORTRAN Library Basic Intrinsic Functions 


Type of Type of 
Function Format Parameter Result Definition 
Exponential EXP(a) Real Real e@ 
DEXP(a) Double Double 
CEXP(a) Complex Complex 
Natural Log ALOG(a) Real Real In(a) 
DLOG(a) Double Double 
CLOG(a) Complex Complex 
Common Log ALOG 10(a) Real Real 10949 (a) 
DLOG 10(a) Double Double 
Sine (angles SIN(a) Real Real sin(a) 
are in radians) DSIN(a) Double Double 
CSIN(a) Complex Complex 
Cosine (angles CcOS(a) Real Real cos(a) 
are in radians) DCOS(a) Double Double 
CCOS(a) Complex Complex 
Tangent (angles TAN(a) Real Real tan(a) 
are in radians) DTAN(a) Double Double 
CTAN(a) Complex Complex 
Hyperbolic Sine SINH(a) Real Real sinh(a) 
Hyperbolic COSH(a) Real Real cosh(a) 


Cosine 
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Table 7-1. FORTRAN Library Basic Intrinsic Functions (Continued) 


Type of Type of 
Function Format Parameter Result Definition 
Hyperbolic TANH(a) Real Real tanh(a) 
Tangent 
Arcsine ASIN(a) Real Real arcsin(a) 
(angles are DASIN(a) Double Double 
in radians) ASIN2(a 1,a2) Real Real arcsin(a1/a2) 
DASIN2(a 1,42) Double Double 
Arccosine ACOS(a) Real Real arccos(a) 
(angles are DACOS(a) Double Double 
in radians) ACOS2(a 1,a2) Real Real arccos(a1/a2) 
DACOS2(a 1,a2) Double Double 
Arctangent ATAN(a) Real Real arctan(a) 
(angles are DATAN(a) Double Double 
in radians) ATAN2(a 1,22) Real Real arctan(a1/a2) 
DATAN2(a 1,a2) Double Double 
Square Root SQRT(a) Real Real (a)% 
DSQRT(a) Double Double 
CSQRT(a) Complex Complex 
Absolute Value CABS(a) . Complex Complex Absolute value of 
complex 
Boolean ICR(a 1,a2) Integer Integer Inclusive OR 
Functions* LAND(a 1,a2) Integer Integer Logical AND 
NOT(a) Integer Integer Logical negation 
IEOR(a 1,42) integer Integer Exclusive OR 
Bit Manipulation+ ISHFT(a 1,a2) Integer Integer Shift a1 by a2 bits 


(shift left logical if a2 
is positive; shift right 
logical if a2 is negative) 
IBTEST(a 1,a2) Integer Logical Test the a2th bit of a1 
(bits are numbered 
right to left, 0 to 15) 
IBSET(a 1,a2) Integer Integer Set the a2th bit of a1 
(bits are numbered right 
to lefi, 0 to 15) 
IBCLR(a 1,a2) Integer Integer Clear the a2th bit of a1 
(bits are numbered right 
to left, 0 to 15) 


*Boolean functions IOR, LAND, NOT, and IEOR are expanded in line. 


+Bit manipulation functions ISHFT, IBTEST, IBSET,.and IBCLR will be expanded in line if a2 
is a constant. 
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Table 7-1. FORTRAN Library Basic Intrinsic Functions (Continued) 


Function 


Absolute Value 


Truncation 


Rounding 


Remaindering 


Choose Largest 
Value 


Choose Smallest 
Value 


Type Conversion 


Format 


ABS(a) 
IABS(a) 
LABS(a) 
DABS(a) 


AINT(a) 
INT (a) 
LINT(a) 
IDINT(a) 
DINT(a) 


ANINT(a) 
NINT(a) 


DMOD(a 1,a2) 
AMOD(a 1,22) 
MOD(a 1,42) 

LMOD(a 1,42) 


AMAXO(a 1,€2,004) 
AMAX 1(8 1,22 yee) 
MAXO(a 1,2 ye0e) 

MAX 1(8 1,82 sooo) 

LMAXO(a 15825000) 
LMAX 1(€ 1,82 yooe) 
DMAX0(a 1,82 ,004) 
DMAX 1(8 1,82,.06) 


AMINO(a 1,22 ,+00) 
AMIN 1(8 1,82 ye00) 
MINO(a 1,22 ye00) 

MIN 1(@ 1,42 000) 

LMINO(a 1,225.04) 
LMIN 1(€ 1,842,002) 
DMINO(a 1,22 yes.) 
DMIN 1(2 1,82,000) 


FLOAT(a) 
REAL (a) 
LFLOAT(a) 
IFIX(a) 

LFIX(a) 
SNGL(a) 
DBLE(a) 
CMPLX(a 1,a2) 
ICHAR(c) 


Type of 
Parameter 


Real 
Integer 
Extended 
Double 


Real 
Real 
Real 
Double 
Double 


Real 
Real 


Double 
Real 
Integer 
Extended 


Integer 
Real 
Integer 
Real 
Extended 
Real 
Integer 
Double 


Integer 
Real 
Integer 
Real 
Extended 
Real 
Integer 
Double 


Integer 
Integer 
Extended 
Real 

Real 

Double 

Real 

Real 
Character*1 


Type of 
Result 


Real 
Integer 
Extended 
Double 


Real 
Integer 
Extended 
Integer 
Extended 


Real 
Integer 


Double 
Real 
Integer 
Extended 


Real 
Real 
Integer 
Integer 
Extended 
Extended 
Double 
Double 


Real 

Real 
Integer 
Integer 
Extended 
Extended 
Double 
Double 


Real 
Real 
Real 
Integer 
Extended 
Real 
Double 
Complex 
Integer 


Definition 


lal 


Nearest whole number 


ai(mod a2)=a1-[a1/a2] *a2 
where [a1/a2] is the 
largest integer<=ABS(a 1/a2) 


Max (8 1,2 y.00) 


Min (8 1,825.0.) 


Convers, fr. Integer to Real 
Convers, fr. Integer to Real 
Convers. fr. Extended to Real 
Convers. fr. Real to Integer 
Convers, fr. Real to Extended 
Convers, fr. Double to Rea! 
Convers. fr. Real to Double 
ai + a2*SQRT(-1) 

Position of character within 
ASCIl collating sequence 
(refer to Appendix E) 
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Table 7-1. FORTRAN Library Basic Intrinsic Funktions (Continued) 


Type of Type of 
Function Format Parameter Result Definition 
Transfer of Sign SIGN(a 1,42) Real Real Sign of a2 times a1 
ISIGN(a 1,a2) Integer integer 
LSIGN(a 1,a2) Extended Extended 
DSIGN(a 1,42) Double Double 
Positive DIM(a 1,42) Real Real a1 - min (a1,a2) 
Difference IDIM(a 1,42) Integer Integer 
LDIM(a 1,22) Extended Extended 
Complex Number CONJG(a) Complex Complex Obtain Conjugate of Comple> 
REALCX(a) Complex Real Obtain Real Part of Complex 
AIMAG(a) Complex Real Obtain Imaginary Pag 
Complex 
Lexical Relation LLT(c 1,¢2) Character Logical Portable lexical 
LLE(c 1,c2) Character Logical comparisons based on 
LGT(c 1,¢2) Character Logical ASCII collating order 
LGE(c 1,¢2) Character Logical 


7.2 Library Subprograms 


The FORTRAN library provides standard subprograms to perform a variety of operations, These 
subprograms are: 


Subprogram Operation 

BUFIN Buffered Input 

BUFOUT Buffered Output 

IUNIT Check Buffered I/O Status 

RECSKP Record Skip 

RANF Get Random Number 

RANSET Preset Random Number Generator Seed 
RANCUR Determine Current Random Number Seed 
Svc Generate Supervisor Call 

SVCFUT Generate File Utility Supervisor Call 
LOC Access Absolute Address 

EXTASK Execute Task 

BIDTSK Bid Task 

DLYBID Delayed Bic Task 

ICRU Read CRU Input 

OCRU Output to CRU 

DATIME Obtain Date and Time 

ADATE Obtain ASCII Date 
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Subprogram Operation 

ATIME Obtain ASCII Time 
MDATE Obtain Military Date 
MFLD Mask Field 

NERRST Find I/O Error 


The following paragraphs explain the function and format of these additional library members. 


7.2.1 Buffered 1/0 


Two subroutines are available to perform buffered input and ouput. Any program may use 

these subroutines by referencing them with a CALL statement. Both of these subroutines 

enable the programmer to implement overlapped I/O. That is, by calling one of these subroutines 
the program can place the I/O operation under control of the called subroutine and return 

to program execution while the I/O operation is in progress, If this technique is used, the 
program must ensure that the previous |/O operation is complete before starting another I/O 
operation with the same unit. This is done by calling the function IUNIT to check the status 

of the unit involved, Only when the previous |/O operation is complete can the program proceed 
with the next I/O operation. The following paragraphs describe the buffered I/O subroutines 

and the I/O unit test status function contained in the FORTRAN library. 


7.2.1.1 Buffered Input 


The BUFIN subroutine transmits one physical record from an I/O device to a prescribed buffer 
area in memory. The form of the call to this subroutine follows. 


FORMAT 


CALL BUFIN(UNIT,MODE,BUFFER,NUMBYTES) 


The required paramters for this call are defined as follows. 


° UNIT — an INTEGER*2 variable or constant that designates the unit number assigned 
to the I1/O device. 


© MODE — an INTEGER*2 variable or constant that indicates the transfer mode 
(0 = ASCII, 1 = binary). 


© BUFFER — a variable name that indicates the starting address of the data transfer, 


° NUMBYTES — an INTEGER*2 variable that contains the number of bytes to be 
transferred, upon input, and the number of bytes actually transferred after an JUNIT 
test indicates successful completion of the operation. 


This subroutine retrieves one physical record from I/O unit number UNIT and stores it in memory 
beginning at address BUFFER using the format indicated by the mode indicator MODE. If the 
value of NUMBYTES is less than the number of bytes in the physical record, only NUMBYTES 
bytes are transferred and the remaining bytes in the physical record are ignored, If the value 

of NUMBYTES is greater than the number of bytes in the physical record, the transmission 

stops at the end of the physical record, 


The status of [/O unit number UNIT must be checked with an IUNIT test before proceeding 
with the operation. When the IUNIT test indicates the BUFIN operation is complete, the number 
of bytes read is returned in NUMBYTES. 


EXAMPLE 


10 


20 


Tele 102 
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DIMENSION IARRAY(40} | 
REWIND i5 | 
DO 20 I=1,100 
ICHAR=8C0 | 
CALL BUFIN(15S, 1, IARRAY. ICHAR?) 


. (Any processing not involving UNIT 15 or 
. ITARRAY} 

ISTAT=IUNIT(15) 

IFCISTAT. EG. 1}¢G0 TO 10 

IFCISTAT. NE. 2>STOP 


CONTINUE 
END 


Buffered Output 


The BUFOUT subroutine transfers one physical record from a buffer area in memory to a specified 
1/O device. The form of the call to this subroutine follows. 


FORMAT 


CALL BUFOUT(UNIT,MODE,BUFFER,NUMBYTES) 


The required parameters for this call are defined as follows. 


° 


= 
UNIT — an INTEGER*2 variable or constant that designates the unit number assigned 
to the I/O device. 


MODE — an INTEGER*2 variable or constant that indicates the transfer mode 
(0 = ASCII, 1 = binary). 


BUFFER — a variable name that indicates the starting address of the data transfer. 


NUMBYTES —~ an INTEGER*2 variable or constant that contains the number of 
bytes to be transferred. 
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This subroutine transfers NUMBYTES bytes from the buffer area in memory beginning at 
location BUFFER to I/O unit number UNIT. The data is written in one physical record using 
the format indicated by the mode indicator MODE, 


The status of I/O unit number UNIT must be checked with a IUNIT test before proceeding 
with the operation. 


NOTE 


The file must be terminated with an END FILE statement before the I/O unit may 
be repositioned. 


~ 


EXAMPLE 


DIMENSION IARRAY (40) 
REWIND i5 
BDO 20 I=1. 100 


ICHAR=80 
CALL BUFOUT(15, 1, LARRAY, ICHAR) 


. (Any processing not involving UNIT 15 or 
IARRAY} 


ic ISTAT=IUNIT(1S) 
IFCISTAT. EG. 1360 TO 10 
IFC ISTAT. NE. 2}STOP 

20 CONTINUE 
END FILE 15 
STOP 
END 


7.2.1.3 Check Buffered I/O Status 


The library function IUNIT is used to return the current status of an I/O operation initiated 
by a BUFIN or a BUFOUT operation. 


FORMAT 


IOSTAT = IUNIT(INUM) 


7/8 


INUM is an INTEGER*2 variable or constant that specifies the |/O unit number. The following 
values are returned by the function. 


Result Definition 
1 Operation incomplete 
2 Operation successful 
3 End of file 
4 Error 


7.2.2 Record Skip 


The subroutine RECSKP allows the user to change position within a sequential or relative 
record file by any number of physical records, forward or backward, The form of the call 
to this subroutine follows. 


FORMAT 
CALL RECSKP(UNIT,NUMREC,IEOF) 
The required parameters for this call are defined as follows. 


° UNIT — an INTEGER*2 variable or constant that designates the unit number 
associated with the file. 


NUMREC — an INTEGER*2 variable which contains the number of records to skip 
upon input and the number of records actually skipped upon return from the subroutine. — 


A positive value indicates forward spacing and a negative value indicates backspacing 
within the file. 


IEOF —— an INTEGER*2 variable that is set to one if an end-of-file is encountered 
or to zero if no end-of-file is encountered. 


This subroutine skips the specified number of records forward or backward within a file. The 
process stops if an end-of-file is encountered. The end-of-file is skipped over, the number 

of records actually skipped (not including the end-fo-file) is returned in NUMREC, and IEOF ww 
is set to 1. If an end-of-file is not reached, RECNUM remains equal to its input value and 

IEOF is set to 0. 


7.2.3 Pseudo-Random Number Generator 


The FORTRAN random number generator uses two Seed values to generate a random number, 
The sequence of random numbers produced has a period equal to 98,304, Three subprograms 
are available to control the random number generator, These allow the user ot generate pseudo- 
random numbers, to set the seed values of the generator io any desired values, and to determine 
the current seed values of the generator. 
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7.2.3.1 Get Random Number 


The library function RANF provides the user with the ability to access a number generated 
by the pseudo-random number generator. The function can be used in an assignment statement. 


FORMAT 
X = RANF( ) 


The parentheses are necessary to cause the compiler to recognize RANF as a function, The 
result of this statement is the assignment of a pseudo-random number to the variable X, The 


random number provided is a REAL* 4 number in the range 0 to + 1. 


7.2.3.2 Preset Random Number Generator Seed 


The library subroutine RANSET allows the user to preset the seed values of the random number 
generator to any desired values. One of the seed values should be odd to avoid a degenerate 


sequence of numbers. The form of the call to this subroutine follows. 
FORMAT 
CALL RANSET(I,J) 


The INTEGER*2 arguments | and J can be any integer values. They designate the values 
to which the random number generator seed is set. If both | and J are equal to 0, the seed 


values are not changed. 


7.2.3.3 Current Random Number Seed 


The library subroutine RANCUR allows the user to determine the current seed values of the 
random number generator. The form of the call to this subroutine follows. 


FORMAT 
CALL RANCURI(I,J) 


The current seed values are returned in the INTEGER*2 arguments | and J. 


7.2.4 Generate Supervisor Call 


The library subroutine SVC allows the user to generate a supervisor call by providing the SVC 
control block. The form of the call to this subroutine follows. 


FORMAT 
CALL SVC (SVC BLOCK,IERR,AUTO ABORT FLAG,REGISTER BLOCK) 
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The required parameters are defined as follows. 


SVC BLOCK — the name of the INTEGER*2 array that contains the SVC control 
block. 


JERR — an INTEGER*2 variable that receives a control byte from the SVC (usually 
an error code). 


AUTO ABORT FLAG == a LOGICAL input variable that when true allows errors 
to abort task. When the call is aborted, argument IERR contains a nonzero value 
and the operating system initiates a traceback. 


REGISTER BLOCK — the name of a 10-element INTEGER*2 array that the SVC 
uses for pseudo-registers. 


7.2.5 Generate File Utility Supervisor Call 


The SVCFUT subroutine allows the user to set up a File Utility SVC call block. The form of 
the call to this subroutine follows. 


FORMAT 


CALL SVCFUT(OPCODE,TYPE,UNIT,PCODE,RECLEN,NUMCHAR,FILENAME,NUMREC, 


ALLOC,IERR) 


The required parameters for this call are defined as follows. 


° 


OPCODE — an INTEGER*2 variable or constant that specifies the utility operation 
code. 


TYPE — an INTEGER*2 variable or constant that specifies the file type. 


1 — Sequential file 
2 — Relative record file 


UNIT —— an INTEGER*2 variable or constant that specifies the unit number. 


PCODE — the passcode feature is not implemented by the current operating 
system; therefore, treat this parameter as a dummy entry. 


RECLEN — an INTEGER*2 variable or constant that specifies the record length. 
The number must be even. Under TAXO, a value of Zero indicates the system 
default record length is to be used. Under TX990, the value 128 is always used 
for sequential files, and a vaiue must be specified for relative record files, 


NUMCHAR — an INTEGER*2 variable or constant that specifies the number of 
characters in the file name. The range is 1 to 49 characters. 


7A TRIUMPH-ADLER 


TA1600/Fort.Ref.Man./0982/e 


7/11 
© FILNAME — a CHARACTER variable, array, or constant that contains the pathname. 
Synonyms are not allowed, 
° NUMREC — an INTEGER*4 variable or constant that specifies the number of records 
in the file. Under TAXO, a Zero value indicates that the system default allocation 


size is to be used. The largest allowable number of records is 2**24, 


° ALLOC — an INTEGER*2 variable or constant that specifies the type of file allocation. 
Under TAXO, a Zero indicates a bounded file, and a 1 indicates an expandable file. 


° JERR — an INTEGER*2 variable that receives an error code. 


7.2.6 Absolute Address 


The library function LOC allows the user to access the absolute address of a single argument. 
The function can be invoked with an assignment statement as follows. 


FORMAT 

| = LOC(B) 
The single argument, B, specifies the parameter whose absolute address will be assigned 
to the variable |. The returned value is INTEGER*2, 


7.2.7 Execute Task 


The EXTASK subroutine allows the user to issue an Execute Task supervisor call. The called 
task must use global luno I/O. The form of the call to this subroutine follows. 


FORMAT 


CALL EXTASK(LUN, TASK,STATID,SUSPND,IERR) 


The required parameters for this call are defined as follows. 


© LUN — an INTEGER*2 variable or constant that contains the logical unit number 
associated with the program file in which the task is installed, 


© TASK — a three-element INTEGER*2 array whose first element contains the installed 
task ID of the task to be executed upon input, and which returns the run-time ID 
of the task. The other two elements contain the task parameters. 
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© STATID — an INTEGER*2 variable or constant which represents the station ID 
to be associated with the called task. A value of 0 indicates that the station ID 
is to be the same as that associated with the calling task. A value of >FF indicates 
that no station is to be associated with the task, 


© SUSPND — a LOGICAL variable which indicates whether the calling task is to 
be suspended until the called task terminates. 


° JERR — an INTEGER*2 variable that receives an error code from the SVC upon 
completion of the task. A value of 0 indicates successful completion. 


7.2.8 Bid Task 
The library subroutine, BIDTSK, allows the user to issue a Bid Task SVC. The called task must 
use global luno |/O. The task specified must be installed nonreplicatable on the system program 
file, as described in the FORTRAN-78 User's Guide for the TAXO Operating Systems. 
The form of the call to this subroutine follows. 
FORMAT 

CALL BIDTSK(TASK,IERR) 


The required parameters for this call are defined as follows. 


° TASK — a three-element INTEGER*2 array whose first element contains the task ID. 
The other two elements contain the task parameters. 


° IERR — an INTEGER*2 variable that receives the task state from the SVC on completion | 
of the call. 


7.2.9 Delayed Bid Task 


The library subroutine DLYBID allows the user to issue a Delayed Bid Task SVC. The called 

task must use global luno I/O. The task specified must be installed nonreplicatable on the 
system program file, as described in the FORTRAN-78 User's Guide for the TAXO Operating 
Systems. The form of the call to this subroutine follows. ww 


FORMAT 


CALL DLYBID(TASK,DATE,IERR) 
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The required parameters for this call are defined as follows. 


° 


TASK — a three-element INTEGER*2 array whose first element contains the task 
ID. The other two elements contain the task parameters. 


DATE — a five-element INTEGER*2 array that contains the date and time at which 
to start the task. Element 1 contains the binary value of the year, element 2 contains 
the binary value of the day; element 3 contains the binary value of the hour; 
element 4 contains the binary value of the minute; element 5 contains the binary 
value of the second. 


IERR — an INTEGER*2 variable that receives the completion code from the SVC 
upon completion. 


7.2.10 CRU Input , 


The library function ICRU allows the user to read from 1 to 16 input CRU lines, The function 
can be used in an assignment statement. 


FORMAT . 


| = ICRU(CRU BASE,NUMBER OF LINES) 


The required arguments are defined as follows. 


© CRU BASE — the base address that is used in the data transfer. This address 
is contained in workspace register 12. 


© NUMBER OF LINES — an INTEGER*2 variable or constant that contains the number 
of contiguous CRU lines to be transferred. 


The value read is assigned to the variable |, right-justified. 


7.2.11 CRU Output 


The library subroutine OCRU allows the user to output data to 1 to 16 output CRU lines, The 
subroutine call appears in the following form. 


FORMAT 


CALL OCRU(CRU BASE,NUMBER OF LINES,VALUE) 


7/14 


The required arguments are defined as follows. 


© CRU BASE — the base address that is used in the data transfer. This address 
is contained in workspace register 12. 


© NUMBER OF LINES — an INTEGER*2 variable or constant that contains the number 
of contiguous CRU lines to be transferred. 


© VALUE — a literal value, expression, or other legal argument whose value (right-justified) 
will be output to the selected CRU lines. 


7.2.12 Obtain Date and Time 


The library subroutine DATIME allows the user to fetch the binary values of the date and time 
parameters supplied by the operating system. The form of the call to this subroutine follows. 


FORMAT 

CALL DATIME(TIME) 
The required argument TIME is a five-element INTEGER*2 array. The first element receives 
the binary value for the year, the second element receives the binary value for the day, the 


third element receives the binary value for the hour, the fourth element receives the binary 
value for the minute, and the fifth element receives the binary value for the second. 


7.2.13 Obtain ASCII Date 


The library subroutine ADATE allows the user to access the date information supplied by the 
operating system and store the values as ASCII characters. The form of the call to this 
subroutine follows. 


FORMAT 

CALL ADATE(DATE) 
The required argument DATE is a three-element CHARACTER*2 array. The fist element 
receives two ASCII numbers representing the month (01 through 12), the second element 


receives two ASCIl numbers representing the day (01 through 31), and the third element 
receives two ASCII numbers representing the year. 
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7.2.14 Obtain ASCII Time 


The library subroutine ATIME allows the user to access the time information supplied by 
the operating system and store the values as ASCII characters. The form of the call to 
this subroutine follows. 

FORMAT 


CALL ATIME(TIME) 


The required argument TIME is a three-element CHARACTER*2 array. The first element 
receives two ASCII numbers representing the hour (01 through 24), the second element 
receives two ASCII numbers representing the minute (00 through 59), and the third element 
receives two ASCII numbers representing the second (00 through 59). 


7.2.15 Obtain Military Date 

The library subroutine MDATE allows the user to access the date information supplied by 
the operating system and store the values in military date format, e.g., 04 JUL 79. The 
routine receives the information in integer format and reformats it in ASCII characters that 
correspond to the required format. The form of the call to this subroutine follows. 
FORMAT 


CALL MDATE(DATE) 
The required argument, DATE, is a CHARACTER*9 variable in which the military date is 
stored. 
7.2.16 Mask Field 
The INTEGER*2 function MFLD allows the user to select a bit, or a group of consecutive 
bits from a 16-bit word and to store these bits right-justified in the value of the function. 
The form of the call follows. 
FORMAT 


| = MFLD(IPOS,NUM, WORD, IERR) 
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The required parameters for this call are defined as follows, 
° IPOS — an INTEGER*2 variable or constant that contains the starting bit position. 
This position is counted from zero, The right-most bit of the word will be the zero 
position. The left-most bit of the word is bit 15. 
The following example shows the bit positions within a word. 


(15——98765432 10) 


To select the desired starting bit, count from. right to left. For example; if the user 
wants bits 4 thru 7, IPOS will equal 4 (not 7), and NUM will equal 4. 


° NUM — an INTEGER*2 variable or constant that contains the number of conse- 
cutive bits to be selected. One bit is the smallest number of bits that can be selected 
and 16 is the largest number. 

° WORD — a 16-bit variable or constant from which the bits will be selected, 

° JERR —= an INTEGER*2 variable that receives the completion code. A returned 
value of zero indicates an error-free function call. A minus one (~1) flags an input 
argument error. 

The error flag (IERR) is set to -1 when any of the three following conditions is not met. 

O0<sIPOS<=15 

1<=NUM<=16 

IPOS+NUM<=16 


The error message INVALID ARGUMENT PASSED TO MELD is given the first time such an 
error is found. 


If the variable IWORD contained the hexadecimal value EF and the user wanted three bits 
starting with bit 4, the following would be true. 


IsMFLD(4,3,IWORD,IERR) 


On return from the function I=26. 


7.2.17 Find I/O Error 


The library function NERRST returns the status of an 1/O operation. NERRST may be used 
with any I/O statement containing the ERR=Sz option. 


FORMAT 


IERR = NERRST(IUNIT) 


7A TRIUMPH-ADLER 


TA1600/Fort.Ref.Man./0982/e 


7/17 


IUNIT is an INTEGER*2 variable or constant that specifies the I/O unit number. The 
following values are returned by the function. 


Result Definition 
0 No error 
1 Iiegal unit 
2 Illegal operation 
3 I/O error 
- End-of-file (only if the END=S, option 
has not been selected; if END=S, 
has been selected, the result is 0) 
5 File previously open 
a ie Invalid format 
12 Overflow on input 
13 Field overflow on output 
14 Illegal input character 
15 Numeric or logical mismatch between 
format and list 
EXAMPLE 


CHARACTER#80 BUFFER 
IUNIT=6 


READCIUNIT, ‘(A) ‘, END=999, ERR=200) BUFFER 


20U —« LERR®NERRST(IUNIT) 


WRITE(7, 225) 1IUNIT, IERR 


999  FORMAT(’ I/0-ERROR, 
STOP 
END 


UNIT = 


‘,13, ‘ ERROR = 


‘, 13) 


If an error is detected on the READ statement, NERRST is used to find the I/O error 


number to be used in the error message. 
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APPENDIX A 
CHARACTER STRING MANIPULATION 
A.1 General 
In addition to the basic FORTRAN operations that may be performed on character data, there 
are several useful functions which make it easier to program applications that involve text 
processing. The character string subprogram package was written to provide users with the 
capability to perform these operations with a single subprogram call. Some of the uses of 
this package are as follows. 
°  Intermixing character strings 
° Editing character data for human error 
° Searching for a particular character set 
° Arranging character data in ascending or descending order. 
This package was written with the intention of allowing the user as much freedom as possible 
in moving, comparing, and validating character strings. The following should be kept in mind, 
NOTE 


The bounds of input and output arrays may be exceeded if care is not 
exercised in selecting and maintaining the indices for these arrays. 


The following subprograms are included in the character string subprogram package. 


Subprogram Description 

SUBSTR Transfer character data from one storage 
location to another 

KOMSTR Compare two character strings for lexical 
ordering 

INDEX Search character string for match to se- 

cond character string 

IVERFY Compare character string for permissable 
characters 

LENGTH Find length of character string 

DNCASE Convert character string to lowercase let- 
ters 
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Subprogram Description 
UPCASE Convert character string to uppercase letters 
TRANS Translate character string 
REPEAT Repeat character string in target character 


string 


A.2 SUBSTR Subroutine 


The SUBSTR subroutine gives the user the capability of transferring character data from one | 
storage location to another in units of one character. Its function is similar to that of the character 
assignment statement with the added capability of referencing substrings within a character 


datum. 


The form of the subroutine call follows. 


FORMAT 


CALL SUBSTR(STRA,IPOSA,LENA,STRB,IPOSB,LENB,IERR) 


The required parameters for this call are defined as follows. 


° 


STRA —= a CHARACTER variable or consiant that contains the characters to be 
moved, 


IPOSA — an INTEGER*2 variable that contains the position of the first character 
to be moved. 


LENA — an INTEGER*2 variable that contains the number of characters to move 
from STRA. 


STRB — a CHARACTER variable that receives the characters, 


IPOSB — an INTEGER*2 variable that contains the starting position in which to 
store the characters. Nea! 


LENB — an INTEGER*2 variable that contains the number of characters to store 
in STRB. 


IERR —— an INTEGER*2 variable that receives one of the following error codes. 


0 ——- Normal completion. 
-1— Invalid argument(s) passed to SUBSTR. 
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If LENA is greater then LENB, only LENB characters will be moved. If LENB is greater than 
LENA, the remaining character positions will be filled with blanks. IPOSA, LENA, IPOSB, 
and LENB must be INTEGER*2 expressions, IERR must be an INTEGER*2 variable. STRA is 
either a character variable or constant. STRB is a character variable, 


IPOSA and IPOSB must have values which are greater than zero, and LENA and LENB must 
have values which are greater than or equal to Zero, If these conditions are not met, the 
following actions will be taken. 


° The message INVALID ARGUMENT PASSED TO SUBSTR will be given the firs time 
an error is found. This message will not be given for subsequent errors in the same 
task. 


° The parameter IERR will be set to minus one and then control will return to the 
calling program without any further action. 


It is the responsibility of the user to ensure that no parameters are passed whose values 
will cause the SUBSTR subroutine to access a character variable beyond its declared length, 


EXAMPLE 
Cc 
Cc THIS PROGRAM INSERTS THE CORRECT DATE 
c INTO A MESSAGE 
c 
CHARACTER#& A 
CHARACTER#i8 B 
INTEGER#2 E 
DATA A/ ‘09/11/7787 
DATA B/ ‘TODAY IS ???P?P?P?P?P??I, 
CALL SUBSTR(A, 1,8, B, 10,9,E) 
c B NOW CONTAINS THE STRING ‘TODAY IS 09/11/78 ‘ 
c A IS UNALTERED 
Cc VARIABLE E CONTAINS THE VALUE ZERO 


END 


A.3 KOMSTR Function 


The INTEGER*2 function KOMSTR allows the user to compare two character strings, and 
returns an INTEGER*2 result which reflects the lexical ordering of the two input strings, The 
collating sequence used is the internal ASCII code of the characters in each string, The 
KOMSTR function is similar to the character relational operations available in the FORTRAN 
language, with the added capability of referencing substrings within a character datum. 
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The form of the call follows. 


FORMAT 


| = KOMSTR(STRA,IPOSA,LEN,STRB,IPOSB,IERR) 


The required parameters for this call are defined as follows. 


STRA —— a CHARACTER variable or constant that contains the first character 
string to be compared. 


IPOSA —= an INTEGER*2 variable or constant that specifies the starting character 
position for the compare in the first character string. 


LEN —_ an INTEGER*2 variable or constant that contains the number of characters 
to compare. 


STRB — a CHARACTER variable or constant that contains the second character\y 
string to be compared. 


IPOSB — an INTEGER*2 variable or constant that contains the starting character 
position for the compare in the second character string. 


IERR —— an INTEGER*2 variable that receives one of the following error codes. 


0 —— Normal completion 


-1 — Invalid input parameter(s) 


The input parameters are checked for positive values. If a parameter is negative or Zero, 

IERR is set to -1, and the function return value is set to zero. The message INVALID ARGUMENT 
PASSED TO KOMSTR will be given the first time an error is found. If the input parameters 

STRA and STRB have an address outside the program space, the end vector will be taken. 

This error may be the result of the input parameters being out of sequence. The input parameters 
IPOSA, LEN, and IPOSB are not checked for exceeding variable boundaries. 


The return values may be 


STRA>STRB,KOMSTR = +1 
STRA=STRB,KOMSTR = 0 
STRA<STRB,KOMSTR = -1 


—— 


{ 
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CHARACTER#i2 B, Ax4 


DATA B/‘123456TESTi2°/ 
A= ‘TEST ’ 


| 
EXAMPLE 


COMPARE STRINGS 


Aaa 


LEN=4 
I=KOMSTR(A, 1, LEN. B, 7, IERR} 


CHECK FOR ARGUMENT-ERROR 
IF(IERR.LT.O}GO TO 100 
MAKE USE OF ARITHMETIC IF 


aan aaa 


IF (I}30, 40, 56 


| STOP 


ERROR-HANDLING 


eOOD 


| STOP 
| END 


This example compares the first four characters of A with four characters of B starting with 
character 7. The return value will be zero, indicating that the strings are equal. The error 
flag will be zero. 
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EXAMPLE 


CHARACTER NAME#S, INAME(46)#2 
DATA NAME/ ‘RAFFERTY ‘7 


READ IN EMPLOYEE NAME AND SARALY TO DATE 


READ(5S, SO, END=998) (INAME( I), I=1, 6), SALARY 
FORMAT (4A2, FiG. 2} 


IS IT RAFFERTY ? 


AAAMKNKANANH 
Qo 


J=KOMSTR (NAME, 1,9, INAME, 1, IERR) 
IF CIERR. EG. -1360 TO 999 
IF(J.EG.0)G0 TO 100 

GO TO io 


Cc 
Cc HANDLE THE RAFFERTY CASE 
Cc 
1 


STOP | 


HANDLE ERRORS 


999 


STOP 
END 


This example searches input data for an employee with the name RAFFERTY. 


A.4 INDEX Function 


The INTEGER*2 function INDEX searches a portion of one character string for a match to 
a portion of another character string. When a match is found, the character number at which 


the matching substring begins is determined. This number is returned as the function value. 
If no match is found, 0 is returned. 
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The form of the call follows. 


FORMAT 


M = INDEX(STRA,IPOSA,LENA,STRB,IPOSB,LENB) 


The required parameters for this call are defined as follows. 


° 


STRA — a CHARACTER variable or constant that contains the string to be searched. 


IPOSA — an INTEGER*2 variable or constant that contains the character position 
in STRA at which the search begins. 


LENA — an INTEGER*2 variable or constant that contains the length of the sub- 
string in STRA to be searched. 


STRB — a CHARACTER variable or constant that contains the string to be matched. 


IPOSB — an INTEGER*2 variable or constant that contains the character position 
in STRB at which the substring to be matched begins. - 


LENB — an INTEGER*2 variable or constant that contains the length of the sub- 
string in STRB to be matched, If LENB>LENA, the value 0 will be returned. LENA 
and LENB are in terms of the number of characters. 


The input parameters IPOSA, LENA, IPOSB and LENB are checked for positive values. Upon 
detection of a negative or zero value, the function returns a value of -1. The message INVALID 
ARGUMENT PASSED TO INDEX will be given the first time this occurs. It is the user's respon- 
sibility to check the parameters being passed against their upper bounds. Execution will con- 
tinue with values that are out of range and may lead to unexpected results, 
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EXAMPLE 
CHARACTER#14 STRA, STRB#4 
DATA STRA/ ‘ABCDEFGHIUKLMN‘/, STRB/ ‘FGHI’/ 
c 
C SEARCH FOR A MATCH 
C 
M=INDEX(STRA, 3, 10, STRB, 2, 3) 
C 
C CHECK FOR ARGUMENT-ERROR 
C 
IF(M. EG. -1}00 TO 100 
STOP 
106 : (Error-handlings) 
STOP 
END 


The preceding program searches 10 characters of STRA beginning at position 3 for a match 
to the three characters in STRB which begin at position 2. The substring GHI is matched, 
returning the value 7 to M. 


EXAMPLE 
CHARACTER#4 A(2),B 
DATA A/‘1234%, ‘5678°/,B/‘456X/ 
Cc 
Cc SEARCH FOR A MATCH 
Cc 
IF (M. EQ. -13G0 TO i106 
STOP 
100 : (Error-handlings) 
STOP 


END 
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The preceding program searches four characters of A beginning at position 5 for a match 
to the first three characters in B. No match to the substring 456 is found, The value 0 is 
returned to M. 


A.5 IVERFY Function 
The INTEGER*2 function IVERFY is designed to compare a character string of data with a 
usersupplied character string composed of all permissible characters, Each character is com- 
pared with the list of valid (acceptable) characters. If an invalid character is located within 
the string, a pointer value is set that is the character number within the string found to be 
invalid. If no invalid characters are found, the pointer value will contain a Zero. 
The form of the call follows. 
FORMAT 

M = IVERFY (STRA,IPOSA,LENA,STRB,IPOSB,LENB) 
The required parameters for this call are defined as follows. 


° STRA = a CHARACTER variable or constant that contains the string to be verified, 


° IPOSA — an INTEGER*2 variable or constant that specifies the beginning character 
of the first string to be verified, 


° LENA — an INTEGER*2 variable or constant that specifies the number of characters 
in the first string to be verified, 


° STRB — a CHARACTER variabie or constant that contains the string against which 
the first string will be checked. 


© |IPOSB — an INTEGER*2 variable or constant that specifies the beginning character 
of the second string. 


° LENB — an INTEGER*2 variable or constant that specifies the number of characters 


in the second string against which the first string will be checked, 


If an invalid argument (a negative or zero parameter) is passed for any of the parameters 


IPOSA, LENA, IPOSB or LENB, the IVERFY routine will output the error message INVALID ARGU- 


MENT PASSED TO IVERFY. IVERFY will subsequently return a -1 to the calling program as the 
value of the function. 
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EXAMPLE 


CHARACTER#i2 ABLE, BAKER 

DATA ABLE/‘ 12345. 17890°/ 

DATA BAKER/‘’ 0123456789. ’/ . | 
i 


VERYFY STRING 

M=IVERFY (ABLE, 1, 10, BAKER, 1, 12) 
CHECK FOR ERROR RETURN 

IF(M. EG. -1360 TO 900 

CONTINUE WITH PROGRAM 


aan oan aaQqn 


STOP _ 


ERROR HANDLING 


naan 


STOP 
END 


In the preceding example, the IVERFY routine would return an 8 for the value of the variable 
M, indicating that an invalid character was detected at the eighth location in string ABLE. 
Note that a blank was considered a valid character. 
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aang iohens) aan 


Maan 


In the preceding example, the user wants to verify only the last four characters contained 
in the string MIKE. The IVERFY routine will return a 0 for the value of M, indicating that the 
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CHARACTER#7 MIKE 

CHARACTER#ii VALID 

DATA MIKE/ ‘6020164 ‘/ 

DATA VALID/ ‘0123456789. ‘/ 
VERIFY DATA 
M=IVERFY (MIKE, 4, 4, VALID, 1. 11> 
CHECK FOR INVALID CHARACTER 
IF(M. NE.0>}GO TO S00 

NO ERRORS. PROCESS DATA. 


WRITE . 


ERROR HANDLING ROUTINES 


STOP. 
END 


characters verified without error. 
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C 
c HANDLE COMMANDS | 
c | 
100 

STOP - 
200 

GO TO 10 
300 

60 TO 10 
acu 

6O TO 10 

END - 


The preceding example is a command processor. The LENGTH function is used to break the 
commands down by length before processing them. 

A.7 DNCASE Subroutine 

The DNCASE subroutine converts a character string to lowercase letters. The character string 
is scanned from left to right replacing uppercase letters with lowercase letters. All other 


symbols in the string remain unchanged. 


The form of the call follows. 


FORMAT 


CALL DNCASE(STRING,LENGTH,IERR) 


The required parameters for this call are defined as follows. 
° STRING — a CHARACTER variable that contains the string to be converted. 


° LENGTH — an INTEGER*2 variable or constant that specifies the number of 
characters in the character string. 
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° |ERR — an INTEGER*2 variable that receives one of the following error codes, 
- 0O-— Normal completion 
- -1— Invalid input parameter(s) 


The input parameter LENGTH is checked to make sure it has a positive value, If the value 

of LENGTH is negative, IERR is set to -1. The message INVALID ARGUMENT PASSED TO 
DNCASE is given the first time an error is found. If the input parameters have an address out- 
side the program space, the end vector is taken. This error may result from the input para- 
meters being out of sequence. The input parameter LENGTH is not checked for exceeding 
the variable boundaries, 


EXAMPLE 


CHARACTER#20 MSG 
DATA MSG/ ‘MESSAGE NUMBER ONE. ‘/ 


FIND THE LENGTH OF THE STRING 
LEN=LENGTH(MSG, 20) 

CONVERT TO LOWER CASE LETTERS 
CALL DNCASE (MSG, LEN, NERR > 
CHECK FOR PARAMETER ERRORS 
IF(NERR.LT.O}GO TO 200 


WRITE OUT THE CONVERTED STRING 
‘i 


aan ang aan oan 


WRITE(&, 100)MSC 
106 FORMAT(’ THE MESSAGE IS ‘, A) 


STOP 
ERROR HANDLING 


STOP 
END 


The preceding example finds the length of the character string in MSG and converts the 
character string to lowercase letters. 


EXAMPLE 


aan AaIaQgranndkn 
Qo 


KONO OAD 


150. 


166 
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CHARACTER#2 BUFFER(40} 
READ IN A LINE OF TEXT 


READ(S, 50, END=100) BUFFER 
FORMAT (40A2?) 


CONVERT TEXT TO LOWER CASE 
CALL DNCASE (BUFFER, 80, IERR) 
WRITE TEXT TO TEMP FILE 


WRITE (4, SOC) BUFFER 
GO TO ic 


READ TEMP FILE BACK IN AND CONVERT 
THE FIRST LETTER OF EACH SENTENCE TO 


UPPER CASE 


END FILE 4 

REWIND 4 

REWIND 5 

FLAG=1 
READ(4, 50, END=999) BUFFER 


FIND FIRST NON BLANK IN THE LINE 


ISTART=1 

LEN®LENGTH( BUFFER, 60) 

IF (LEN. EG. 0)G0 TO 400 

J=INDEX (BUFFER, ISTART, 1,‘ ‘, 1,1} 
IF (J. EG. 0)GO TO 165 
ISTART=ISTART+1 

GO TO i160 


CAPITALIZE FIRST LETTER OF A LINE 
IF (FLAG. EG. 0>}GO TO 250 


CALL UPCASE (BUFFER, ISTART, IERR?> 
FLAG=0 
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FIND END OF SENTENCE 


YWaan 


So I=INDEX (BUFFER. ISTART, Si-ISTART, ”. ‘,1,3) 
IF(I. EG@.0)G0 TO 406 
I=I+3 
IF(I. LE. LEN)GO TO 300 
FLAG~1 
GO TO 400 


CAPITALIZE FIRST LETTER OF A SENTENCE 


Qaan 


OO N=I/2 
L=2-(I-N*2} 
IF(L. EG. 1} N=N+i 
CALL UPCASE(BUFFER(N),L, IERR} 
ISTART=I+1 
IFCISTART. GT. 80}G0 TO 400 
LEN=LEN-I 
GO TO 250 


WRITE UPDATED FILE 


CO WRITE(S, SO) BUFFER 
GO TO 150 

999 END FILE 5 

STOP 

END 


The preceding example first converts the uppercase text in an input file to lowercase text. 
Then the first character in each sentence of the text is converted to uppercase. Note that 
this example uses several character string functions. 


A.8 UPCASE Subroutine 


The UPCASE subroutine converts a character string to uppercase letters. The character string 
is scanned from left to right, replacing lowercase letters with uppercase letters. All other 
symbols in the string remain unchanged. 

The form of the call follows. 


FORMAT 


CALL UPCASE(STRING,LENGTH,IERR) 
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The required parameters for this call are defined as follows: 
° STRING — a CHARACTER variable that contains the string to be converted. 


° LENGTH — an INTEGER*2 variable or constant that specifies the number of 
characters in the character string. 


° JERR — an INTEGER*2 variable that receives one of the following error codes, 
- O-— Normal completion 
- -1— Invalid input parameter(s) 


The input parameter LENGTH is checked to make sure it has a positive value. If the value 
of the input parameter LENGTH is negative, IERR is set to -1. The message INVALID ARGU- 
MENT PASSED TO UPCASE is given the first time an error is found. If the input parameters 
have an address outside the program space, the end vector is taken. This error may result 
from the input parameters being out of sequence. The input parameter LENGTH is not check 
for exceeding variable boundaries. 


EXAMPLE 


CHARACTER#20 MSG 
DATA MSG/ ‘message entry #1'/ 


FIND THE LENGTH OF THE STRING 
LEN=LENGTH(MSG, 20) 

CONVERT TO UPPER CASE LETTERS 
CALL UPCASE (MSG, LEN, IERR} 
CHECK FOR PARAMETER ERROR 

IF (IERR. EG. -1)60 TO 100 


WRITE OUT THE CONVERTED STRING 


aan aan aaa ANH 


WRITE(&, SO}MSG 
FORMAT(’ THE MESSAGE IS ‘, A) 


15) 
Qo 


STOP 


TA1600/Fort.Ref.Man./0982/e 


a 


A/19 


ERROR HANDLING 


STOP 


The preceding example finds the length of the character string in MSG and converts the 
character string to uppercase letters. 


EXAMPLE 


aaa aagagaMKKaAng 
oo 


CHARACTER#4 BUFF (20) 
LB=8O 


READ IN A LINE OF TEXT 


READ(S, 50, END=999 ) BUFF 
FORMAT (20A4) 


CONVERT TEXT TO UPPER CASE 
CALL UPCASE (BUFF, LB, IERR>} 
WRITE TEXT TO TEMP FILE 


WRITE (4, 50} BUFF 
cO TO 10 


STOP 
END 


The preceding example converts the lowercase text in an input file to uppercase text and 
stores the results in another file. 
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A.? TRANS Subroutine 


The TRANS subroutine translates a character string. The original character string is scanned 
from left to right comparing each character with the first translate string. If a match is found, 
the corresponding character in the second translate string is used to replace the character 
in the original character string. 
The form of the call follows. 
FORMAT 

CALL TRANS(STRING,LENGTH,TSTRA,TLEN,TSTRB,IERR) 
The required parameters for this call are defined as follows. 


° STRING — a CHARACTER variable that contains the string to be translated. 


° LENGTH — an INTEGER*2 variable or constant that specifies the number of 
characters in the character string being translated. ww 


° TSTRA —a CHARACTER variable or constant that contains the comparison trans- 
late string. 


° TLEN ——an INTEGER*2 variable or constant that specifies the number of characters 
in the translation character strings. 


° TSTRB — a CHARACTER variable or constant that contains the replacement 
character string. 


° JERR — an INTEGER*2 variable that receives one of the following error codes. 
0 — Normal completion 
-1 — Invalid input parameter(s) 


The input parameters are checked for certain boundary conditions. If the value of the input 
parameter LENGTH or TLEN is negative, IERR is set to -1. The message INVALID ARGUMENT 
PASSED TO TRANS is given the first time an error is found. If the input parameters have an Wy 
address outside the program space, the end vector is taken. This error may result from the 

input parameters being out of sequence. The input parameters LENGTH or TLEN are not checked 
for exceeding variable boundaries. 
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CHARACTER BUFFER#30, TA#4, TB¥4 

DATA BUFFER/ ‘THE PROGRAM CLEVELA] IS ‘’OK’‘! ‘/ 
DATA TA/ ‘EIS SE SA 

BATA TB/‘C3". °/ 


WRITE OUT INITIAL BUFFER 


WRITE (6, 5O}BUFFER 
FORMAT( ‘1iTHE INITIAL BUFFER IS: “1 A} 


CHANGE C TO ¢ 
2 TO } 
e TO st] 
re . 


CALL TRANS(BUFFER, 30, TA, 4, TB, IERR) 
WRITE OUT THE TRANSLATED BUFFER 
WRITE (46, LOG) BUFFER, IERR 

FORMAT ( ‘OTHE NEW BUFFER IS: ‘,&,/' TERR = ‘, 12) 


STOP 
END 


The preceding example translates the character string in BUFFER according to the translate 
strings TA and TB with the following results. 


THE PROGRAM (LEVELA) IS ""OK", 


FORMAT 


A.10 REPEAT Subroutine 


The REPEAT subroutine allows the user to repeat a character string in a target character 
string. 


The form of the call follows. 


CALL REPEAT(STRA,LENA,STRB,LENB,RNUM,IERR) 
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The required parameters for this call are defined as follows. 
° STRA — a CHARACTER variable that contains the target character string. 


° LENA — an INTEGER*2 variable or constant that specifies the length of the target 
character string. 


° STRB — a CHARACTER variable or constant that contains the character string 
to be repeated. 


° LENB — an INTEGER*2 variable or constant that specifies the length of the 
character string to be repeated. 


° RNUM — The INTEGER*2 variable or constant that specifies the number of times 
to repeat the character string. 


° JERR — an INTEGER*2 variable that receives one of the following error codes. 
0 —— Normal completion 


-1— Invalid input parameter(s) 


The input parameters are checked for certain boundary conditions. If the value of the input 
parameter LENA, LENB, or RNUM is negative, IERR is set to -1. Also if the length of the 
character string to be repeated multiplied by the repetition number is greater than the length 
of the target string, IERR is set to -1. The first time an error is found, the message INVALID 
ARGUMENT PASSED TO REPEAT will be given. If the input parameters have an address out- 
side the program space, the end vector will be taken. This error may result from the input 
parameters being out of sequence. The input parameters LENA or LENB are not checked for 
exceeding variable boundaries. 
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EXAMPLE 
CHARACTER#2 BUFFER (20), STARS(20), STAR#4 
DATA STAR/ ‘He / 
LEN=80 


INITIALIZE PRINT BUFFER 
CALL REPEAT(BUFFER. SO. ‘ ‘, 1,80, IERR) 


CHECK FOR PARAMETER ERROR 


aaa ann 


IF(IERR. LT.0)GD TO 500 

CALL REPEAT(STARS, 80, STAR, 4, 20, IERR) 
IF(IERR.LT.0)@0 TOSSO 

STOP 


HANDLE ERRORS 


Maan 


STOP 
END 


The preceding example uses the REPEAT subroutine to initialize buffers used within the 
program. 
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EXAMPLE 
CHARACTER#10 NUMBER, NEWNUM 
Cc 
Cc ACCEPT A NUMBER THAT IS LEFT JUSTIFIED 
c 


1c DISPLAY (5, 15, LINE=S5, ERASE} 
a ) FORMAT ( ‘ENTER NUMBER: vi 
ACCEPT (5S, 25, LINE=5, POSITION=16, PROMPT }NUMBER 
25 FORMAT (A10} 


Cc 
Cc RIGHT JUSTIFY THE NUMBER 
Cc 
L=LENGTH(NUMBER, 10} 
IF(L.EG@.0>GO TO 999 
J=10-L 
K=J+1 
CALL SUBSTR (NUMBER, 1, L, NEWNUM, K. L, IERR) 
CALL REPEAT(NEWNUM, 10, ‘ ‘, 1, J, IERR? 
GO TO i0 
999 STOP 
END 


The preceding example accepts a number which is left-justified in the field and right justifies 
it before processing the number. 
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APPENDIX B 
ISA EXTENSIONS. 
B.1 General 


The FORTRAN library includes a set of subroutines that satisfy the extensions to the FORTRAN 
language recommended by the Instruments Society of America (ISA Extensions S61.1-1975 
and S61.2-1976). These subroutines permit interface with executive programs, process input 
and output functions, provide access to time and date information, and provide procedures 

for file access and control of file contention. The following paragraphs describe the function 
and calling sequence of each of the ISA extension functions provided in the FORTRAN-78 
library. 


The error codes returned by the ISA extensions may be found in table B-1. The value zero 
indicates a normal completion of the ISA extension. 


Unit numbers used in the ISA extensions are associated with logical unit numbers generated 
internally to FORTRAN. This mapping is transparent to the user. The unit number is maintained 
in the file control block (FCB) while the logical unit number (LUNO) is stored in the physical 


record block (PRB) for a given file. Therefore, when debugging a program, remember that a 
unit number used in an ISA extension will not match the LUNO assigned to the file. 


Table B-1. FORTRAN Error Codes 
Code Error Condition Modules 
FORTRAN internal Error 


>FO00 FCB table full. The maximum 
is 20 OPENW 


ISA Extension Errors 


>F 100 The specified task ID does not START/TRNON 
exist. 
>F101 The time specified illegal, START/TRNON 
>F102 The system is unable to bid START/TRNON 
the task. 
>F103 The delay is too long. The WAIT 
maximum delay is 1638.375 
seconds. 
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Table B-1. FORTRAN Error Codes (Continued) 


Code 


>F104 


>F105 


>F 106 


>F107 


>F108 


>F109 


>F10A 


>F10B 


>F10C 


>F10D 


Error Condition 


Illegal time unit specified. The 
range is 0 to 3. 


The number of bit groups is zero. 


The CRU address is invalid. The 
CRU address must be less 
than > 1FFe,. 


The number of bits is invalid. The 
range is 1 to 16, 


Operation incomplete, hardware 
failed to respond. 


Invalid channel number, The range 
is 0 to 63. 


LUNO requested is invalid. The 
range is 1 to 99, 


File previously opened and not 
closed. 


Access method is invalid. The 
range is 1 to 4. 


The number of characters in the 
file pathname is out of range. The 
range is 1 to 49, 


The ISA extension subroutines are listed in the following. 


Subroutine 


START 
TRNON 


WAIT 
DIW 


DOLW 
DOMW 


Description 


Modules 
WAIT 
DIW/DOLW/DOMW 
DIW/DOLW/DOMW/ 
AISQW/AOW/AIRDW 
DIW/DOLW/DOMW | 
AISQW/AOW/ AIRD\) 
AIRDW 

| 
OPENW/CLOSEW/RDRW/ 
WRTRW/MODAPW | 


OPENW 


OPENW 


OPENW 


Specify time delay before execution of task 
Specify time of day to start execution of 


task 


Specify time delay before continuing exe- 


cution of program 


Input bits from CRU interface and store in 


array 
Output bits to CRU 


Output momentary digital output signal bit 


groups to CRU 
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Subroutine Description 

DATE Determine calendar date 

TIME Determine time of day 

AISQW Analog input in sequence 
AIRDW Analog input in random order 
AOW Analog output 

CFILW Create file 

DFILW Delete file 

OPENW Open unit number with file name 
CLOSW Close unit number with file name 
RDRW Read relative records 

WRTRW Write relative records 

MODAPW Change access method of opened file 


B.2 Start a Program 


The START subroutine allows the user to specify a time delay before beginning execution 

of a specified task. On TAXO, the specified task must be installed nonreplicatable on the 
system program file, as described in the FORTRAN-78 User's Guide for the TAXO Operating 
Systems. The called task must use global LUNO I/O. The form of the call follows. 


FORMAT 


CALL START(TASK,DELAY,UNITS,IERR) 


The required parameters for this call are defined as follows. 


° 


TASK — a three-element INTEGER*2 array. The first element contains the task 
ID of the task to be started. The last two elements contain the four bytes of task 
parameters, i.e. terminal I/O, data, etc. 

DELAY —- an INTEGER*2 variable or constant that indicates the number of units 
that the task will be delayed before starting. A negative or zero value results in 
an immediate start. 


UNITS — an INTEGER*2 variable or constant that specifies the units to be used 
in measuring the time delay. The following values indicate the units. 


— 0 = 8.3-millisecond units (system clock), minimum of six units 


— 1 = one-millisecond units (delay is computed in 50 millisecond intervals), 
minimum of 50 units 


— 2 = one-second units 
— 3 = one-minute units 


IERR —- an INTEGER*2 variable that receives an error code. 
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B.3 Starta Program at a Specified Time 


The TRNON subroutine allows the user to specify a time of day at which to start execution 
of a specified task. The specified task must be installed nonreplicatable on the system pro- 
gram file, as described in the FORTRAN-78 User's Guide for the TAXO Operating Systems, 
The called task must use global LUNO I/O. The form of the call to this subroutine follows. 
FORMAT 

CALL TRNON(TASK,TIME,IERR) 

The required parameters for this call are defined as follows. 

° TASK — a three-element INTEGER*2 array. The first element contains the task 
ID of the task to be started. The last two elements contain the four bytes of task 
parameters. 

° TIME — a three-element INTEGER*2 array. The three elements of the array contain 
the hour, minute, and second values that define the military (24-hour) time at why 
the program will be started. 

° JERR — an INTEGER*2 variable that receives an error code. 

NOTE 


The TRNON subroutine is not supported on TXDS or TX5DS systems 


B.4 Delay Continuation of a Program 


The WAIT subroutine allows the user to specify a time delay before continuing with the exe- 
cution of a program sequence. The form of the call to this subroutine follows. 


FORMAT 


CALL WAIT(DELAY,UNITS,IERR) 
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The required parameters for this call are defined as follows. 


° 


DELAY — an INTEGER*2 variable or constant that indicates the number of units 
that the program sequence will be delayed. A negative or zero value results in 
no program delay. 


UNITS —— an INTEGER*2 variable or constant that specifies the units to be used 
in measuring the time delay. The following values indicate the units, 


— 0 = 8.3-millisecond units (system clock), minimum of six units 


— 1+ one-millisecond units (delay is computed in 50 millisecond intervals), mini- 
mum of 50 units. 


— 2 «= one-second units 
— 3 = one-minute units 


IERR —— an INTEGER*2 variable that receives an error code. 


B.5 Digital Input 


The DIW subroutine allows the user to input sets of bits from the CRU interface and store 
those bits in a specified array. The form of the call to this subroutine follows. 


FORMAT 


CALL DIW(NUMBER,SPECIFICATION, TARGET, IERR) 


The required parameters for this call are defined as follows. 


NUMBER — an INTEGER*2 variable or constant that indicates the number of 
bit groups to be input from the CRU. The value must be 1 to 16. 


SPECIFICATION — a 2 by "number" INTEGER*2 array. The first element of each 
pair is the decimal CRU base as it would be used in register 12, and the second 
element is the number of bits to transfer. 


TARGET — an INTEGER*2 array that is "number" elements long into which the 
bit groups are stored right-justified. 


IERR —— an INTEGER*2 variable that receives an error code. 
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B.6 Latched Digital Output 


The DOLW subroutine allows the user to output bits of information to the CRU. The bits can 
be latched in either the set or the reset condition. The form of the call to this subroutine 


follows. 


FORMAT 


CALL DOLW(NUMBER,SPECIFICATION,SOURCE,MASK,IERR) 


The required parameters for this call are defined as follows. | 


° 


NUMBER — an INTEGER*2 variable or constant that indicates the number of bit 
groups output to the CRU. The value must be 1 to 16. 


SPECIFICATION —— a 2 by "number" INTEGER*2 array. The first element of each 
pair is the decimal CRU base as it would be used in register 12, and the second 
element is the number of bits to transfer. 


SOURCE — an INTEGER*2 array that contains the output bit groups that are 
fetched. The data in the array must be right-justified. 


MASK — this parameter is ignored. MASK is provided for compatibility with other 
ISA libraries. 


IERR —— an INTEGER*2 variable that receives an error code. 


NOTE 


CRU input and output data is bit-reversed. For a detailed description 
of the CRU refer to the Assembly Language Programmer's Guide. 


B.7 Momentary Digital Output 


The DOMW subroutine allows the user to output bit groups to the CRU. The bit groups consist 
of momentary digital output signals. All specified bits are reset after a delay of TIME. The 
form of the call to this subroutine follows. 


FORMAT 


CALL DOMW(NUMBER,SPECIFICATION,SOURCE,TIME,IERR) 


4 
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The required parameters, except TIME, have the same definitions as those for the latched 
digital output subroutine, DOLW. TIME specifies a time delay before continuing the execution 
of a program sequence. The parameter TIME is a three-element INTEGER*2 array whose 
elements correspond to the parameters of the WAIT routine (refer to paragraph B.4). 


B.8 Obtain DATE 


The DATE subroutine allows the user to determine the correct calendar date. The form of 
the call to this subroutine follows. 


FORMAT 

CALL DATE(I) 
The required argument | is an INTEGER*2 array containing three elements. The first element 
receives an integer representation of the year (A.D.), the second element receives an integer 


representation of the month, and the third element receives an integer representation of the 
day. 


B.9 Obtain TIME 


The TIME subroutine allows the user to determine the correct time of day. The form of the 
call to this subroutine follows. 


FORMAT 
CALL TIME(J) 


The required argument J is an INTEGER*2 array containing three elements. The first element 
receives the hour, the second element receives the minute, and the third element receives 
the second, 


B.10 Analog Data Handling 


The 990 computer family provides for analog data handling by an optimal analog-to-digital 
(A/D) conversion module and an optimal digital-to-analog (D/A) conversion module. These 
hardware modules can be controlled from FORTRAN by the following ISA subroutines. 


Subroutine Purpose 

AISQW Analog Input in Sequence 
AIRDW Analog Input in Random Order 
AOW Analog Output 
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The AISQW subroutine samples points in an order determined by the A/D hardware module. 
Each call samples consecutive channels beginning with channel 0 on a single A/D module. 

A maximum of 64 channels can be sampled since an A/D module contains no more than 64 
channels. The subroutine samples a channel and stores a 12-bit integer in the least significant 
bits of the input array. Any scaling of values is left to the user. The choice of voltage or 
current range and code scheme (straight binary or two's complement) is a hardware switch 
selectable option left to the user. 


The AIRDW subroutine samples points in an order determined by the user. Each call samples 
any number of channels from any A/D module sequence in any order. The user specifies the 
number of sample points desired and provides an INTEGER*2 array to receive the converted 
values. In addition, each point must be tagged with a CRU A/D module address and a channel 
number. The subroutine samples the specified channel on the specified module and stores a 
12-bit integer in the 12 least significant bits of the input array, Any scaling of values is left 

to the user. The choice of voltage or current range and code scheme (straight binary or two's 
complement) is a hardware switch option left to the user. If two's complement is switch-selec- 
ted, the converted value is-not sign-extended to 16 bits; it remains a 12-bit positive integer. 
To sign-extend the value, the following statement can be used, 


IF (IBTEST(IVALUE,11)) THEN IVALUE = IOR(IVALUE,4ZFO000) 
In the preceding statement, IVALUE contains the 12-bit value. 


The AOW subroutine converts digital values in an order specified by the user. The digital-to- 
analog hardware module can be ordered with one to four channels. Each call converted any 
number of values through any combination of D/A modules using any combination of channels. 
The user specifies the number of values to convert and provides an INTEGER*2 array containing 
the values. Each value must be tagged with a CRU D/A module address and a set of channel 
enable bits. The subroutine transfers the 12 least significant bits of a value to the D/A hard- 
ware module and channels. All scaling is left to the user. The choice of voltage or current 

range and code scheme is a hardware switch option left to the user. 


All three subroutines have a similar calling sequence. The form of the call for the AISQW sub- 
routine follows. 


FORMAT 
CALL AISQW(POINTS,CRUADDR,INVECTOR,IERR) 
The required parameters for this call are defined as follows. 


° POINTS — an INTEGER*2 variable or constant that specifies the number of con- 
versions to perform. 


© CRUADDR — an INTEGER*2 variable or constant that indicates the CRU address 
for the chassis slot containing the A/D module. This is the value stamped on the 
chassis for the slot. 
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INVECTOR —— an INTEGER*2 array that receives the converted values. 


IERR —— an INTEGER*2 variable that receives an error code. 


The form of the call for the AIRDW subroutine follows. 


FORMAT 


CALL AIRDW(POINTS,ADDRCHAN, INVECTOR,IERR) 


The required parameters for this call are defined as follows. 


° 


POINTS —- an INTEGER*2 variable or constant that specifies the number of con- 
versions to perform. 


ADDRCHAN — a two-dimensional INTEGER*2 array that tags the source for each 
conversion. The dimensions should be 2 and the number of points. For each two- 
word entry, the first word contains the A/D module CRU address and the second 
contains the channel number to use on the A/D module. For example, for point J, 
the entry is 


ADDRCHAN(1,J)=CRU address of A/D module 
ADDRCHAN(2,J)=Channel number 


Note that channel numbers are in the range 0 to 63. 
INVECTOR — an INTEGER*2 array that receives the converted values. 


IERR ——_ an INTEGER*2 variable that receives an error code. 


The form of the call for the AOW subroutine follows. 


FORMAT 


CALL AOW(POINTS,ADDRCHAN,OUTVECTOR, IERR) 


The required parameters for this call are defined as follows. 


° 


POINTS — an INTEGER*2 variable or constant that specifies the number of con- 
versions to perform. 


° 


° 
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ADDRCHAN — a two-dimensional INTEGER*2 array that tags the destination for 
each conversion. The dimensions should be 2 and the number of points. For each 
two-word entry, the first word contains the D/A CRU module address and the 
second contains the channel enable mask. The channel enable mask is stored in 
the four least significant bits of the array word. The bits are defined as follows. 


Bit Number 


15 


14 


13 


12 


Purpose 


1 = Enable conversion on channel 0 


0 = Inhibit conversion on channel 0 


1 = Enable conversion on channel 1 
0 = Inhibit conversion on channel 1 


1 = Enable conversion on channel 2 
0 = Inhibit conversion on channel 2 


1 = Enable conversion on channel 3 
0 = Inhibit conversion on channel 3 


OUTVECTOR — an INTEGER*2 array that contains the digital values to convert. 


IERR —— an INTEGER*2 variable that receives an error code. 


B.11 Create File 


The CFILW subroutine allows the user to create a specific file. The form of this call follows. 


FORMAT 


CALL CFILW(TYPE,NUMCHAR,FILENAME,RECLEN,NUMREC,IERR) 


The required parameters for this call are defined as follows. 


° 


TYPE — an INTEGER*2 variable or constant that specifies the file type. 


1 — Sequential file 


2 — Relative record file 


NUMCHAR — an INTEGER*2 variable or constant that specifies the number of 
characters in the file name. The range is from 1 to 49 characters. 


FILENAME —— a CHARACTER variable, array, or constant that contains the file 
pathname. Synonyms are not allowed. 
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© RECLEN — an INTEGER*2 variable or constant that specifies the record length 
in bytes, The number must be even, Under TAXO, a value of zero indicates the 
system default record length is to be used. 

° NUMREC — an INTEGER*4 variable or constant that specifies the number of 
records in the file. Under TAXO, a Zero value indicates an expandable file using 
the system default for the allocation size. A positive value specifies a bounded 
file, and the largest allowable number of records is 2**24, The file type is non- 
contiguous and expandable. 


° |IERR —— an INTEGER*2 variable that receives an error code. 


B.12 Delete File 
The DFILW subroutine allows the user to delete a specified file. The form of this call follows. 
FORMAT 


CALL DFILW(NUMCHAR, FILENAME, IERR) 


The required parameters for this call are defined as follows. 


° NUMCHAR — an INTEGER*2 variable or constant that specifies the number of 
characters in the file name. The range is from 1 to 49 characters. 


° FILENAME — a CHARACTER variable, array, or constant that contains the file 
pathname. Synonyms are not allowed. 


° JERR — an INTEGER*2 variable that receives an error code. 


B.13 OPEN/CLOSE Unit Number 


The subroutine OPENW/CLOSEW allows the opening and closing of a unit number with an 
associated file name. A call to OPENW does not create a nonexistent file. The form of these 
calls follows. 
FORMATS 

CALL OPENW(UNIT,NUMCHAR,FILENAME,ACCESS,IERR) 


CALL CLOSEW(UNIT,IERR) 
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The required parameters for these calls are defined as follows, 


° 


UNIT — an INTEGER*2 variable or constant that specifies the unit number. 


NUMCHAR — an INTEGER*2 variable or constant that specifies the number of 
characters in the file name. The range is from 1 to 49 characters. 


FILENAME == a CHARACTER variable, array, or constant that contains the file 
pathname. Synonyms are not allowed, 


ACCESS — an INTEGER*2 variable or constant that specifies the access method. 


Relative Record File Sequential File 

1 — Read only 1 — Read only 

2 — Shared 2 — Read/Rewrite 
3 —— Exclusive all 3 —— Exclusive 

4 —— Exclusive write 4 —— Exclusive 


IERR —— an INTEGER*2 variable that receives an error code. IERR may be omitted 
and automatic error termination with traceback will result if an error occurs. If IERR 
is included, it is the user's responsibility to code an error-handling routine to handle 
values returned in IERR. A file must have been previously opened with ISA Extension 
OPENW before it can be closed or have automatic error termination with traceback. 


| B.14 Read/Write Relative Records 


This subroutine allows the user to read and write relative records. The file to be read or 


written must have been previously opened by the ISA Extension OPENW and must be closed 
by the ISA Extension CLOSEW. The forms of these calls follow. 


FORMATS 


CALL RDRWUNIT,RECNUM,BUF,WRDCNT,IERR) 


CALL WRTRW(UNIT,RECNUM,BUF,WRDCNT,IERR) 
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The required parameters for this call are defined as follows. 


UNIT —— an INTEGER*2 variable or constant that specifies the unit number. 


RECNUM — an INTEGER*4 variable or constant that specifies the record 
number to be read or written, The maximum value is 2**24, 


BUF — a variable or array of any data type into which information is read or 
written. 


WRDCNT — an INTEGER*2 variable or constant that specifies the maximum 
number of words to be read or written. 


IERR —— an INTEGER*2 variable that receives an error code. IERR may be 
omitted and automatic error termination with traceback will result if an error 
condition occurs. 


B.15 Change File Access Method 


The subroutine MODAPW allows the user to change the access method of a previously 
opened file. The form of the call follows. 


FORMAT 


CALL MODAPW(UNIT,ACCESS,IERR) 


The required parameters for this call are defined as follows, 


° 


° 


° 


UNIT —- an INTEGER*2 variable or constant that specifies the unit number. 


ACCESS —— an INTEGER*2 variable or constant that specifies the access method. 


Relative Record File Sequential File 

1 — Read only 1 — Read only 

2 — Shared 2— Read/Rewrite 
3 — Exclusive all 3 — Exclusive 

4 — Exclusive write 4 — Exclusive 


IERR —— an INTEGER*2 variable that receives an error code, IERR may be 
omitted and automatic error termination with traceback results if an error occurs. 
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To obtain more representable values, i.e. more precision, floating point is necessary. For 
example, a variable declared 


REAL*4 PIG 
would cover the range 
-7.237 E 75 <= PIG <= -5.397 E -79 


and 


and 
5.397 E -79 <= PIG <= 7,237 E 75 


On the number line, this is represented as: 


eee tr 
~7.237E75 -5.397E-79 0 5.397E-79 7.23E75 
and in expanded view 


<~ t t hH t 


5+ 2° 16-6 -4999999404 5 5000000596 .5+2*°16**-6 


Note that the range is larger than that of fixed-point numbers and that the number of values 
is much larger. There are about 4.3 billion (2**32) distinct values resulting in much better 
precision. Also note there is a set of values about Zero that are not representable, namely 
any X, such that 


ABS(X) <5.397 E -79, X not equal zero. 


D.3 Fixed-Point Number Operating Characteristics 


The user of fixed-point numbers must be aware of their operating characteristics, since loss 
of precision is far easier than with other data types such as INTEGER or REAL. The following 
example demonstrates the dangers involved. 


EXAMPLE 


FIXED (-7) A 

FIXED (4) B 

FIXED (0) C 

4=45G6-7 

B=27G4 

C=A+B 

WRITE(6, 800034, B.C 
8006 FORMAT(1X, SF20. 8) 

END 
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Table E~1. ASCIi Character Set (Continued) 
Hexadecirna! Decima! Holierith Code 

Value Value Character (Punched Card) 
22 34 - 8-7 

23 35 # 8-3 

24 36 $ 11-8-3 
25 37 % 0-8-4 
26 38 & 12 

27 39 ‘ 8-5 

28 40 ( ’ 12-8-5 
29 41 ) 11-8-5 
2A 42 ™ 11-8-4 
2B 43 + 12-8-6 
2c 44 - 0-8-3 
20 45 - 11 

2E 46 : 12-8-3 
2F 47 / 0-1 

30 48 0 0 

31 49 1 1 

32 50 2 2 

33 51 3 3 

34 52 4 4 

35 53 5 5 

36 54 6 6 

37 55 7 7 

38 56 8 8 

39 57 9 9 

3A 58 : 8-2 
3B 59 : 11-8-6 
3C 60 - 12-8-4 
3D 61 = 8-6 
3E 62 > 0-8-6 
3F 63 ? 0-8-7 
40 64 @ 8-4 
41 65 A 12-1 
42 66 B 12-2 
43 67 Cc 12-3 
44 68 D 12-4 
45 69 E 12-5 
46 70 F 12-6 
47 71 G 12-7 
48 72 H 12-8 
49 73 | 12-9 
4A 74 J 14-1 
4B 75 K 11-2 
4C 76 L 11-3 
4D 77 M 11-4 
4E 78 N 11-5 
4F 79 e) 11-6 
56 80 P 11-7 
51 81 Q 11:8 
52 82 R 11-9 
53 83 Ss 0-2 


> 


» vieericiwan./U982/e 
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Table E-1. 


Hexadecimal 
Value 
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ASCII Character Set (Continued) 


Decimal 
Value 


Character 


2. 
—“--*NSXEK CH MHRA DVODZH—- Ke Fa HrTOAAO TMH 1 Se rmrnN<xs<ecad 


DEL 


Hollerith Code. _:«.~ 
(Punched Card) 


0-3 
0-4 
0-5 
0-6 
0-7 
0-8 
0-9 
12-8-2 
0-8-2 
11-8-2 
11-8-7 
0-8-5 . 
8-1 3. 
12-0-1 ; 
12-0-2 

12-0-3 

12-0-4 

12-0-5 

12-0-6 

12-0-7 

12-0-8 

12-0-9 

12-11-1 

12-11-2 

12-11-3 ‘ 
12-11-4 ° 
12-11-5 at 
12-11-6 oe 
12-11-7 - 
12-11-8 

12-11-9 

11-0-2 “ 
11-0-3 2 
11-0-4 

11-0-5 

11-0-6 

11-0-7 

11-0-8 

11-0-9 

12-0 

12-11 

11-0 

11-0-1 

12-9-7 


——————————————EEEEEeEeEeEeEeeeeEeEeee 


———eeee 
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To obtain more representable values, i.e. more precision, floating point is necessary. For 
example, a variable declared 


REAL*4 PIG 
would cover the range 
-7.237 E 75 <= PIG <= -5.397 E -79 


and 


and 
5.397 E -79 <= PIG <= 7.237 E 75 


On the number line, this is represented as: 


a 
-7.237E75 -5.397E-79 0 5.397E-79 7.23E75 
and in expanded view 


—— 
5 +2*16-6 .4999999404 5 5000000596 .5+2*16**-6 


Note that the range is larger than that of fixed-point numbers and that the number of values 
is much larger. There are about 4.3 billion (2**32) distinct values resulting in much better 
precision. Also note there is a set of values about Zero that are not representable, namely 
any X, such that 


ABS(X) <5.397 E -79, X not equal zero. 


D.3 Fixed-Point Number Operating Characteristics 


The user of fixed-point numbers must be aware of their operating characteristics, since loss 
of precision is far easier than with other data types such as INTEGER or REAL. The following 
example demonstrates the dangers involved. 


EXAMPLE 


FIXED (-7) A 

FIXED (4) B 

FIXED (0) C 

4=450-7 

B=2704 

C=A+B 

WRITE(6, 860034, B,C 
800G FORMAT(1X, SF20. 6) 

END 
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Execution of this program results in 
5760.00000000 B=1.68750000 C=1665.00000000 


At first glance, the value of C in the preceding example appears to be incorrect. The expected 
value for C is 5761.68750000, which is the value produced under floating-point arithmetic. 
Actually, the most significant bit was lost because of the alignment of fixed-point numbers 
before arithmetic operations. The FORTRAN-78 compiler generates a command to shift the 
contents of A 11 bits to the left and to align the binary point of A and B before generating 

the add A and B command. This shift command loses the most significant bit of A. The follo- 
wing illustrates the process. 


Original value of A: 0000 0000 0010 1101 
Shifted value of A: 0110 1000 0000 0000 
Original value of B: 0000 0000 0001 1011 
Result of addition: 0110 1000 0001 1011 


By converting the result from binary to decimal, the correct answer is 1665.6875. To avoid 
such a loss of precision, table D-1 gives the smallest and largest magnitude numbers repre- ww 
sentable with a given scale factor. 


Table D-1. Fixed Point Scale Factor Vs. Decimal Magnitude 


Scale Smaliest Decimal Number Largest Decimal Number 
nec 4454661 28730D-0F 1525832 34012127D-04 
o 30 .VELS2257 440-09 3051 6464480242540-04 
@ 29 1184626451 49D-08 "BL dese s604S50S0-04 
G Se 2725 29027SD-08 » 1220665872097021-02 
27 . 74505S05970-05 » 24413317441 94030-03 
o 24 » 147011411 9D-07 RSS 2L634S823806D-02 
a 35 » PYROSTSSESD—-07 9765326974677 461 20-03 
G 34 » 9404644 78D-07 1952065295 255220-02 
m@ Do LLiSZOPVZE9ED-04 > 3904130770710450-02 
mh 2D .2324125771D-04 »7S1 2241581 42090D-02 
i 24 1476937 15520-04 1 1546245231428418D-01 
bt ZO P67 431 44D-06 13124 704632546836D-01 
@ 19 > 190734S4233D-05 » 42498092651 34720-01 
& 16 LSE1449724660-08 124996185 3027240 00 
@ 17 "342929453 10-05 TSASSS2 3706054490 O0 
mM 14 » LS2RS7S9060-04 479984741 2109370 OO 
Q@ 15 » 3051 7572120-04 9SSS494ER24218750 OO 
G14 41025 154625D-04 LippPoaee44R43750 O1 
i 13 112207031 250-03 "RSSSE7 7029687500 Of 
ew a2 »24414046250D-02 799 S9s7500D O1 
Q 11 , 488:26:1 2500D-03 115% 187500 O2 
@ 10 » P745625000D-03 «3199 3750000 OZ 
mS "49521 25000n-02 , Gere 7=S0000D OZ 
a * 200L250000L-02 SLS7ePs SOOOOL! OF 
mF "351 2500000D-02 tate a SOOOOOL! OF 
a & .15625000000-01 ~511?84375000000D O2 
a § so 12S5OO000OD-O1 © {OBB248 750000000 04 
fi 64 ie SEOOOOOOOD—O 1 "2047 927500000000 04 
a 3 ,L250000000D OO 40758 7SQQOOO0000B 04 
rm = a SEGOOOOOOOD oOo PHLPLTEOOOOOOOOOD 04 
re) ¥ SSOOOOOOOOOL OO ~ 1422250 OHOOOOOOOOL OS 
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Table D-1. 


Smallest Decimal Number 


D/5 


Fixed Point Scale Factor Vs. Decimal Magnitude (Continued) 


a ———— 


Largest Decimal Number 


» 10000000001 01 »327670000000000D O5 
a =e ° 8600000000 O1 -4552400000000000 OF 
@ =2 - 40000000000 O1 ~1310689000000000D 06 
q -3 ~ 20000000000 O1 »262136000000000D 04 
@ -4 - 14000000002 O2 »524272000000000D 06 
Q -5 ~SZ200000000D O2 .104854400000000D 07 
Q@ -& - 64000000002 02 » 2Z09708S00000000D 07 
Q -7 - 12800000001 OF -419417400000000D 07 
QO -8 ~ 25600000001 O3 »S29S352000000000 07 
a -9 - 51200000001 OF -1477670400000001 O8 
G-10 » 10240000001 04 > 335534020000000D 0s 
D-11 »2048000000D 04 -6710681460000000D O08 
G-12 » 40960000000 04 ~124213632000000D 0o¢ 
f-13 21220000000 04 » 2624272464000000D OF 
G-14 ~1638400000B 05 - 5268545230Q00000 OF 
Q-15 ~ 32762000000 OF »107370205400000D 10 
G=-16 - 65534600000 OF »214741511200000D 106 
-17 -13107200001D 04 e4294S3622400000D 10 
15 - 26214400001 04 ~S58947244500000D 10 
n-19 ~S242eS0000D 04 ~171793442940000D 11 
B-ZO »1042574000L 07 -3435846297920000D 11 
-21 ~ 20971520000 07 -6871723795240000D 11 
G-22 ~ 41943040001 07 ~137434757148000D 12 
-25 ~EBRES40S000L 07  .27486951S3346000R 12 
Gi- 24 ~1477721400D OF »5497390366720000 12 
A-25 ~SPSS44S2000 O8 -109947S07234400D 13 
Zé -67108846400D 08 »21989546144468500D 13 
-27 »1342177280D OF 439791 2293376002 13 
H-22 26543545401 OF ~87958245246475200D 12 
H-29 ~SaéS70P7120D oF «175916491735040D 14 
I-30 ~1073741824D 10 eSTLSS29E3470080D 14 
B-Z1 ~21474834648D 10 2 703665966940160D 14 

D.3.1 Assignment Statement a 


In the assignment statement 


X2Y 
X has a scale factor of sx, and Y has a scale factor of sy. 
If sx=sy, no shifting precedes the MOV command which stores the integer portion of Y in X. 


Or if sx<sy, a shift right arithmetic (SRA) of Y sy-sx bit places precedes the MOV. This shift 
loses sy-sx of the least significant bits, 


Or, if sx>sy, a shift left arithmetic (SLA) of Y sx-sy bit places precedes the MOV. This shift 


loses sx-sy bits from the most significant bits. If these bits were leading zeros, no information 
is lost. 


A point to remember in each of these cases is that the resulting value retains the scale factor 
of sx. An example of a common source of precision loss is 


FIXED (0) X 
X = 54Q3 
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Since the scale of the right-hand side (3) is greater than the scale of the left-hand side (0), 
the right-hand side is right-shifted three bits. Unfortunately, this results in a loss of precision, 
since 54 (110110 in binary) right-shifted three bits is 6 (110 in binary). 


To avoid problems with assignment statements, be sure to scale the receiving item with a 


Scale factor equal to or larger than the scale factor of the source item. 


D.3.2 Addition and Subtraction 


Addition and subtraction are treated identically for fixed-point operations; therefore, all 
references in this paragraph to addition apply equally to subtraction. In the expression 


X+Y 
X is FIXED (sx), and Y is FIXED (sy). 
If sx=sy, no operations precede the normal addition. 


Otherwise, the operand whose scale factor is smaller is snifted left ABS(sx-sy) bit places. 
This aligns the X and Y binary points, but loses sx-sy leading bits from the shifted operand. 
The scale factor of the result is the MAX(sx,sy). 


D.3.3 Multiplication 


Multiplication of two fixed-point numbers requires no shiiting before the normal integer multi- 
plication. The only source of precision loss is a product overflowing 15 binary digits. The 
scale factor of a product is the sum of the scale factors of the operands. 


D.3.4 Division 


Division has the same restrictions as multiplication. The only difference is that the resulting 
scale factor is the difference of the operand scale factors, 


D.4 Scaling to Avoid Loss of Precision 


The fixed-point number iQs represents the number i/2**s. This is analogous to the floating- 
point number jEt representing the number j*10**t. The difference is that the fixed-point 
number is divided by the scale factor while the floating-point number is multiplied by the 
scale factor. This is the most important point to remember in scaling fixed-point numbers. 


There are severa! guidelines to follow in scaling fixed point numbers to avoid loss of precision. 
© Try to keep the scale factors of the right-hand and left-hand sides of assignments 
nearly equal. For example, as the results of a calculation grow larger, scale tempo- 
raries up to larger values; if they grow smaller, scale temporaries down. 


© Avoid using the shift functions in fixed-point expressions. 


© Analyze and arrange expressions so that intermediate values have nearly the same 
scale factors. 


ne 
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APPENDIX E 


ASCII/HOLLERITH CHARACTER CODES 


The ASCII character set used by FORTRAN-78 is shown in table E-1. The table shows the 
hexadecimal value and decimal value of the ASCII code for each character, the character 
itself and the corresponding Hollerith (punched card) code. 


Table E-1. ASCII Character Set , HS, 


Hexadecimal Decimal Hollerith Code. 
Value Value Character (Punched Card) 
00 i) NUL 12-0-9-8-1 
01 1 *SOH: “FS 12-9-1 - 
02 2 STX 12-9-2 
03 3 7 ETX+' 12:9-3. 
04 4 EOT.--:-. 9-7 
05 5  ENQ' - = 0-9-8-5 
06 6 ACK 0-9-8-6 
07 7 BEL 0-9-8-7 
08 8 BS 11-96. , 
09 9 HT 12-9-5 - 
0A 10. a ou iia 0-9-5 
OB 1 ee eae Spee 
0G 2 Ag VERS dng oo, DORE Cw geet nee 
oD 13 - ee SS et greg ten ak ee 
) 14 so 12-9-8-6 
OF 15 Ss! 12-9-8-7 
10 16 DLE 12-11-9-8-1 © © 
1 7 DC1 11-9-1 
12 18 38 TS on! Bee te ar 11-9-2 a : 
13 99 ise -152 ‘<DC3 = 11-9-3 Vice? Gan 
14 20 - DCs 9-8-4 
15 21 NAK 9-8-5 
16 22 SYN. - 9-2. 
17 23 =: 0-9-6 
18 Be tees, CAN . 11-9-8 
‘s 19 a ee EM . 11-9-8-1 
> 1A ee. aa 9-8-7 
Ey 1B es .., £SS , 0-9-7 
S 1c 28 Seg 11-9-8-4 
3 1D ye GS 11-9-8-5 
o 1E 30° ** 50 38 ge 11-9-8-6 
= 1F 31 US 11-9-8-7 
2 20 32°24: SPACE None 
_ 21 <\ a a. i 128-7 
6 
~ 
=) . 
3 
< ih on ; 
i—_ ee I SQ SS yes 
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Table E~1. ASCIi Character Set (Continued) 


Hexadecirma! Decima! Holierith Code 
Value Value Character (Punched Card) 
22 34 a 8-7 
23 35 # 8-3 
24 36 $ 11-8-3 
25 37 % 0-8-4 
26 38 & 12 
27 39 : 8-5 
28 40 ( ‘ 12-8-5 
29 41 ) 11-8-5 
2A 42 : 11-8-4 
2B 43 + 12-8-6 
2C 44 ; 0-8-3 
2D 45 ~ 11 
2E 46 , 12-8-3 
2F 47 / 0-1 
30 48 0 0 
31 49 41 1 
32 50 2 2 
33 51 3 3 
34 52 4 4 
35 > 53 5 5 
36 54 6 6 
37 55 7 7 
38 56 8 8 
39 57 9 9 
3A 58 : 8-2 
3B 59 x 11-8-6 
3C 60 < 12-8-4 
3D 61 = 8-6 
3E 62 > 0-8-6 
3F 63 ? 0-8-7 
40 64 @ 8-4 
41 65 A 12-1 
42 66 B 12-2 
43 67 Cc 12-3 
44 68 D 12-4 
45 69 E 12-5 
46 70 F 126 
47 71 G 12-7 
48 72 H 12-8 
49 73 | 12-9 
4A 74 J 11-1 
4B 75 K 11-2 
4C 76 L 11-3 
4D 77 M 11-4 
4E 78 N 11-5 
4F 79 ie) 11-6 
56 80 P 11-7 
51 81 Q 11:8 
52 82 R 11-9 
53 83 S$ 0-2 
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Table E-1. ASCII Character Set (Continued) 


Hollerith Code». «:«..~ 


xadecimal Decimal 
™ Value Value Character (Punched Card) 
54 84 T 0-3 
55 85 U 0-4 
56 86 V 0-5 
57 87 Ww 0-6 
58 88 xX 0-7 
59 89 Y 0-8 
5A 90 z 0-9 
5B 91 [ 12-8-2 
5C 92 \ 0-8-2 
5D 93 ] 11-8-2 
5— 94 A 11-8-7 ‘3 
5F 95 0-8-5 Oh 
60 96 \ 8-1 * 
61 97 a 12-0-1 
62 98 b 12-0-2 
63 99 re 12-0-3 
64 100 d 12-0-4 
65 101 e 12-0-5 
66 102 f 12-0-6 
67 103 9 12-0-7 
68 104 h 12-0-8 
69 105 i 12-09 
6A 106 j 412-11-1 
68 107 x 12-11-2 
6C 108 | * 42-44-3 , 
6D 109 m 12-11-4 ‘ 
6c 110 n 12-11-5 i. 
6F 111 0 12-11-6 we 
70 112 p 42-11-7 
71 113 q 12-11-8 
72 114 r 12-11-9 
73 115 s 11-0-2 
- 116 t 11-0-3 
75 117 u 11-064 
76 118 v 11-05 
77 119 w 110-6 
78 120 x 11-0-7 
= 121 y 11-08 
7A 122 z 11-0-9 
7B 123 { 12-0 
7C 124 12-14 
7D 125 \ 11-0 
7E 126 ~ 11-0-1 
7F 127 DEL 12-9- 


TA1600/Fort.Ref.Man./0982/e 


7 Al TRIUMPH-ADLER 


A 


TA1600/Fort.Ref.Man./0982/e 


INDEX 


A Field Specification 
Absolute Address Function 
ACCEPT Statement 


I/| 


5-11 
7-11 


4-17; 5-30 


Access Input/Output, Sequential 5-19 
Access Method Subroutine, Change 

File B-13 
ADATE Subroutine 7-14 
Addition Module: 

Addition of Fixed-Point Numbers D=6 
AIRDW Subroutine B-7 
AISQW Subroutine B-7 
Alphanumeric Data Specification 5-10 
Analog Data Handling B=7 
Analog Input in Random Order 

Subroutine B-7 
Analog Input in Sequence Subroutine B-7 
Analog Output Subroutine B-7 
Analog-to-Digital Conversion Module B-7 
AOW Subroutine B-7 
Arguments; 

Mathematical Function 6-5 
Subprogram Dummy 2-22 
Arithmetic: 

Assignment Statement 3-11 
Expression Evaluation 3=4 
Expressions 3-1; 3-3 
IF Statement 4-5 
Operators 3-1; 3-2 
Relational: 

Operations 3-7 
Operators 3-7 
Array: 

Assumed-Size Dummy 2-17 
Bounds Checking 2-18 
Array Element 2-15 
Arrays 2-15; 6—4 
Memory Representation of 2-16 
ASCIl Character Codes Appendix E 
ASCII Date Subroutine 7-14 
ASCII Time Subroutine 7-15 
ASSIGN Statement 41; 4-2 
Assigned: 

Format Labels 5-18 
GO TO Statement 4e2 


Assignment: 


Statement 

Arithmetic 

Character 

Logical 

Type Conversion 
Assumed-Size Dummy Array 
ATIME Subroutine 


B Compiler Option 
BACKSPACE Statement 
Bid Task: 

Subroutine 

Delayed 

BIDTSK Subroutine 
Blank Common 
BLOCK DATA; 
Statement 
Subprogram 

Block IF Statement 
BN Edit Descriptor 
Bounds Checking, Array 
BOUNDS: 

OFF Statement 

ON Statement 
Buffered: 

Input Subroutine 

1/0 

1/O Status Subroutine 
Output Subroutine 
BUFIN Subroutine 
BUFOUT Subroutine 
BZ Edit Descriptor 


CALL Statement 

Carriage Control 
Characters 

CFILW Subroutine 

Change File Access Method 
Subroutine 

Character Codes: 

ASCIl 

Hollerith 

Character Constants 
Character Data, Transferring 
CHARACTER: 

Assignment Statement 
Carriage Control 

Data Type 

Formats, Referencing 
Relational: 

Operations 


3-11; D=5 
3-11 
3-16 
3-11 

T3=12 
2-17 
7=15 


2-18 
5-1; 5=26 


7-12 
7-12 
7-12 
2-21 


6-8; 6-10 
6-10 
4a7s hn 2 
5-13 
2-18 


2-18 
2-18 


7-5 
7-5 
7-7 
7-6 
7-5 
7-6 
5-13 


1-8; 4-15; 6-6; 7-5 


5-10 
T5-10 
B-10 


B13 


Appendix E 
Appendix E 
2-12 
A-2 


3-16 
T5-10 
1-9 
5-18 


3-8 


‘Operators - = 
Set, FORTRAN 
Statement 
Nees String: 
apts . Length 
* Repeat 
Subprogram Package 
Compare 
Convert to Lowercase _ 
_. Convert to Uppercase. 
‘e- ‘Search 
“” Translate 
Check Buffered 1/0 Status Function 
Close Unit Number-Subroutine 
CLOSEW Subroutine 
‘Comment Lines 


’ 


! 
or my 


‘ 


T1-3; 1-43 1-8 | 


a >. ‘Common: 
3-3-2 EOE 2=27° 
ach ~ Labeled . 20; 2-23 - 
¢ 2 e_- Statement 2-17; 2-18; 2-19, 
Abas Compare Character Strings _ An3y Ando ® 
a Compilation, Conditional T1-3; 1-9 
Aa Compiler Option: , 
ga sn2 © 2-18 . 
OF spate as eae 1-9" 
baa Complex Data Specification 5-15 | 
oe Complex Numbers. 2-12" 
sie COMPLEX Statement ~ 2-2; 2-4 
ibe Computed GO TO Statement - 4-2. | 
mag Conditional Compilation T1=3; 1-9°* © 
e- Console Display — Bog 5-30... 
le Constants .. 2-1; 2-2 | 
gs Character . 2-12 
hs Hexadecimal 2-8 
7 Logical 2-13 
or on Continuation Line T1-3; 1-3 
ye CONTINUE Statement 4=8; 4 5 
OC a3 Control Statements 44 P 
is we CoMversion: 
— Assignment Statement Type T3=12° 
; Character String to Lowercase A-14 
or Character String to Uppercase A=17 
3% Double-Precision Format 5-7 
, Fixed-Point Format 5-6. 
a Floating-Point Format 5=6 
- Hexadecimal Format 5-8 
me Input/Output Type 5-t 
Integer Format 5-3 


Module: 

Analog-to-Digital Be 
Digital-to-Analog Be 
Floating-Point to Integer 4 
Integer to Floating-Point Fe! 
Rea! Format 5-3; 5+ 
Single-Precision: ] 
COPY Statement 1-1 
Create File Subrouiine B-7 
CRU Input Function 7-" 
CRU Output Subroutine 7-1 
Current Random’ Number Seed 

Subroutine 7- 
D Field Specification 5-2; 5-7; 54 
‘Data Specification: - . 
Alphanumeric : 5-" 
Complex 5-1 
Literal 5-" 
Logical ww 5. 
Numerical 5. 
DATA Statement 1-43 6+ 
Data Type, CHARACTER 1+ 
Data Types 2- 
Date and Time Subroutine 7-1 
DATE Subroutine - 
Military 7-7 
DATIME Subroutine 7-) 
Declaration i 2- 
Definition Statement: 1 
Function 6+ 
Statement Function 6" 
Delay Continuation Subroutine B. 
Delayed Bid Task Subroutine 7-) 
Delete File Subroutine Be 
DFILW Subroutine | Be 
Digital Input Subroutine By 
Digital Output Subroutine: 

Latched 3 BY 
Momentary Be 
Digital-to~Analog Conversion Module By 
DIMENSION Statement 2-16; 2-20; 6—4; 6: 
Dimensioned. Variables 21 
Direct Access Input/Output 5-| 
Display: 
High/Low Intensity oT 
Statement 5-! 
Division Module: 
Division of Fixed-Point Numbers D 
DIW Subroutine B 
DLYBID Subroutine 7=| 
DNCASE Subroutine Ae! 
Loop 1-9; 4G: nd 


| 


i/il 


Nested stead 4-12 Logical ugly 3-7} 3-9 
Statement ~~. he Bs Gad Within Parentheses, Evaluating * * Pane 3-6 
DOLW Subroutine = ss. | B-6 _  EXTASK Subroutine ~ (xa LS 7-11 
DOMW Subroutine _..- viet: B-6 External: We siti 
DOUBLE PRECISION Statement 2-2; 2-4 . File Characteristics a1 Ng 5-20 
Double-Precision: a ae ag Statement aly 2-3; 6-11 
Format Conversion we 5-7 Subprogram 66s 6-11 
Numbers 7 erie 2-11... Subprograms lace S 6-1 
Memory Representation 7 sty ijt Eta 
of 2-12 F Compiler Option the ae - = 1-9 
Dummy: 2 rd i : F Field Specification © a 5325 5-3; 5-8 
Arguments, Subprogram. 2-22 : Factors, Scale jmp a7 ne Das 
Array, Assumed-Size _ Fait Field, Source Record _ 3 pan : T1-3 
identifiers 6-1; 6-2. Field Specification: °°) 
i 4 A la tecniteende 5-11 

E Field Specification  ,. ... 5=2;5-5;5-8 .. D.- "Seah breeder aes On 7s 5-8 
Edit Deseripterss "S55 SSS a oor oe . napa dhs 5-55 5-8 
BN ~ eal besek 5-13. F ——— 52} es 
BZ Bp ae Batt nc naiacgi ome So” 
ELSE Block case, SeBrestB- te, thes: i sey ‘ate 5-12 
ELSE IF Block a « ‘eewee st ae Asad etaaa,s Sate 
ELSE: sa Site A ,.; ., Repeated © YON oe Ron a 5-16 
IF Statement-- =... ,.. 3 a T wat. ¥en 5-14 
Statement _— iP i hn? x tein AS oy te 5-13 
END FILE Statement . “ie 5-1; 5-26 © Zz ~ 5-2; 5-8 
END: aren ere a File Access Method Subroutine ae 
iF Statement meesiretrreee 7-1) os Change ~*~ : ed B-13 
Statement He 4-18. . Eile Characteristics: are anal. oc 
End-of-File Record” Racer 1-2 “"* External ‘emai? a tal ae - tees 5-20 
EQUIVALENCE Statement. "=. 2-18; 2-23 ._.- Internal © "7 5-20 
Civer Cadag: (ht ret aes ale Sequential. aia dean wate cottccten ek} 5-19 
ISA Extension |... Wee ee TB-{. . File Name} 'Interrial © soe nea 5-20 
Error Function, Find. Vo. : pene 7-16 ° File Positioning Statements - amet ne 5-25 
Evaluating Expretnane-) Within sie bs " i File Utility Supervisor Call 3 = ° - 
Parentheses Dos Sa. Subroutine, Generate ~~ ~> 7-10 
Evaluation: oe he ‘a yee. = Bind: Rig 
Arithmetic Expression. ae: apt 3-4 °°. Character String Length“. ~ 2) >< A=12 
Exponent... “t an 1-9. ~ I/O Error Function Bre = Yep 7-16 
Logical Expression——~=~ Ba aig 3-10 Statement siesta alae: 3 : ne 5-30 
Execute Task Subroutine... ........ Fa44 FIXED Statement ss. Feiaeeeeeaniea 2-2; 2-6 
Exponent Evaluation - . | en dit? 1-9 Fixed-Point: * Pong eee 

ai Exponentiation, ~ j=. j=, ¢¢ Se5 . Format Conversion . ss ae eas lenis 5-6 

3 Expression Evaluation: ,., . Se Integers ~ ~~ ee 2-13 

2 Arithmetic KEM Soeyek 2416 3-4 Memory Representation. Es ae 2-14 

J Logical ee 3-10 | Number: thas 

c Expressions tee gee ee 3-1 Addition of P itcce ces aa oe D-6 

£ Arithmetic aera 2g 3-1; 3-3 Applications , a ee oe per De1 

> , - . Division of thea tiies D6 

ag al ie eae Multiplication of Bev eS Ve Dé 

é Aer 

~~ 1 Sp F 

S * YOR SA. Gg 

mu =4 ii 

< ~ es? 

= s0C—m 


\ ‘ 


we Operating Characteristics 
Precision 
ee Subtraction of _., 
25 Floating-Point; 
{2 Division: 
Format Conversion 


g2* Ot a3 
rt £3 


_” | if) = <2 -Multiplication: 
a Number Precision 
a Format: 
a Conversion: .. - 
slik Double-Precision 
aris Fixed-Point 
Floating-Point 

5 Hexadecimal thos 
Sage: Integer * 2B eens 
aaa 3 Real ming re 
“ Free Field .. (a8 ebook 
i Format Labels, Assigned. 
ois Format Specification, Stored ~ 
ay 5 FORMAT Statement. _ . 

' Formats, Rasesanerng ¢ Character 
Sao Gwo-sk Formatted: ie 
; READ Statement 
eo Records er 

j WRITE Statement, . - _-. 
— FORTRANS - 5 - -..0 0). - 
Character Set , 
i Language oe 


“77 1° ' + 9*~Braction Addition Module, ~~ - 


; Program 


Statement Numbers 
Statements 1-2; 
FORTRAN-66 VS FORTRAN-78 
Free Field Format s 
Function 

Absolute Address 7 
Arguments, Mathematical 

Check Buffered |/O Status 

CRU Input 7 
Definition: 

Statement 

Statement, Statement 

Find I/O Error 7 
Get Random Number 

ICRU 7 
INDEX 
Intrinsic 
JUNIT 
IVERFY 
LENGTH a 
Loc 7 
Mask Field 7 
Mathematical . 

MFLD = ies 

NERRST 

RANF 

Reference, Mathematical 

Statement 

Subprogram .. ... - 

G Field Specification 

Generate: . 

File Utility Supervisor 

Call Subroutine 

Supervisor Call Subroutine 
Generator, Pseudo-Random Number 
Get Random Number Function 
GO TO Statement 

Assigned 

Computed .-...... 
Unconditional . 


6-1; 6-11; 


— 
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Graphics TENTS 5-30; 5.33 
Group Specification, Repaiiay 52 5-16; 5-16 
: Te: Soa 
H Field Specification * AC PG. 3 5-12 
Hexadecimal: Met oS eae 
Constants “Gas D> 22-8 
Format Conversion **-" > © 5-8 
High/Low Intensity Display 5-33 
Hollerith Character Codes” Appendix E 
| Field Specification oorg | it 5-2; 5-3 
ICRU Function - ae ee 7513 
Identification Field ~ ~ 12° 8i- 4-3 
Identifier = ae 2-1; 2-2 
Identifiers, Dummy’ Pe, 6=1;-6-2 
IF Block SS he Ts He 12 
IF Level cui 3 -aQ 2 
IF Statement mc - hed 
Arithmetic Wierd 5° hed 
Block TPA 4-74 bn t2 
ELSE 2.3 * che 
END 4734-13 
Logical sn 4-6 
IMPLICIT Statement ~~!" "5 2~33-2-6 
INDEX Function “3 52 Ab 
Indexing 1/O Lists Thad 2522 
Lists ALES TSI IE Aa. apd :,» 524 
Subroutine, Buffered “5°! *5°: 7-5 
Input/Output:Consolé’ Display 5-30 
Direct Access **" 5-27 
Sequential Access =< fae. ne 5-19 
Statements “°° FS be ee 3 5-1 
Type Conversion-' © ~°'* > 5-1 
Instrument Society” of *” 
America(ISA)": ” t 1-1 
Integer 7S 2-7 
Conversion: ‘'S ~~. > 
Fixed-Point < 2-13 
Memory Representation’ ‘of 2-14 
Format Conversion’ “°° "°° 5-3 
Statement 2-2; 2-3 


Integer to: 


~ WV 


Single-Precision: 


’ Intensity Display, High/Lewe' 5-33 
Internal: 1B! GS 3 

File Characteristics iS MOMS atte 2 5-20 
File Name eee 5-20 
Subprogram 6-1; 6-2 
Intrinsic; $2439 

Function’ “°° < if54 ts 16— 13 Ge 115 FI 
Statement ERLE ge ett 2-35 6-11; 7-1 
IUNIT Function ~~=!°2"~ at 7-7 
ISA 15 -Mg= 1-1 
Extension Error Codes: - : TB~1 
Extensions SITET 21 Lok B-1 
IVERFY Function wed? a A-9 
Buffered Bete, 7-5 
Error Function, Find 7-16 
Lists, Indexing Bt* 5-22 
I/O Status sulrossling, " %ay= 

Check Bulferéd ty es 7-7 

BUT Ba ten 193 
KOMSTR Subréutine me TAN 3 A-3 


2-20; 2-23 


Subroutine B-6 
LENGTH Function ” A-12 
Library sopprgeais. 7-4 
Lists: 

Indexing I/O Bee 5.2 5-22 
Input 5-21 
Output 5-21 
READ 5-21 
WRITE 5-21 
Literal Data Specification 5-10 
LOC Function 7-11 
Logical: 

Assignment Statement 3-11 
Constants 2-13 
Data Specification 5-9 
Expression Evaluation 3-10 
Expressions 3-1; 3-73 3-9 
IF Statement 4mé 
Operations 3-9 
Operators 3-1; 3-9 
Statement 2-2; 2-5 
Loop, DO 1-93 4-8; 4-9 
Lowercase, Convert Character 

String to A-14 
Machine Code 1-1 
Mask Field Function 7-15 


Mathematical Function: ... 


Military Date Subroutine” 


Relational 


_ Arguments 
. Reference 2 6=5 

“ MDATE Subroutine: 2. 2 > 5. 7.4% 
Memory Representation of: rae eS PP aad? 
Arrays ot pie. Ok 2-16 - 
Double-Precision Numbers eat ga2 2e42: :; 
Fixed-Point Integers .. ae op 2-14 
Real Numbers ee ee 2-10 : 
MFLD Function a ee 7=15.- 


2% 


.MODAPW Subroutine- ... B-13 
"' Module: ; -enee ee. eid wc 
Analog-to-Digital Conversion 5k et ae Be7 . 
Digital-to~Analog Conversion __ ae? ee 
Subroutines: ~ seal stivied “ae 
Momentary Digital Output, > 3. -. on™ 
Subroutine oe val ht B-6 
rt “Multiplication Module: oprtti + los 
, Multiplication of Fixed-Point ZAZA ge 
Numbers niet D=6.. 
; ~ NERRST Function pions 7=16 
a Nested DO sre iers 4m 12 - 
Normalization sion. 2-10. 
_ Numerical Data Specitications =: Ag 52... 
~ Object Code Benim 1-1 
OCRU Subroutine eee 3 7-13 
OPEN: _— ee Saat 
Statement ES 5=27- 
_ Unit Number Subroutine ruth : Bei1. 
- OPENW Subroutine TET a Be44 
: Operator Precedence 3-25 3~10} T3-1 
._ Operators ) 364: : 
Arithmetic - B=13 3-2 
_ 3, Relational 3e7° 5 
‘” Character Relational 43 ey 3-8 
Logical = Tom lg BSS. 


7 5315 An 


esaerenrc 1. 


Output: 
Lists 
Subroutine: 
Analog 
Buffered 


PAUSE Statement 
Power-of-Ten Form 
Precision: 

Fixed-Point Number 
Floating-Point Number 
Preset Random Number Generator 
Seed Subroutine: ~= 
Program: 

Preparation 

Size nese SS : 
Programming: cachaaine 
PROM Programming: 


Pseudo-Random Number Generator 


RANCUR Subroutine: °°: 
Random Number Function 
Random Number Generator Seed 
Subroutine, Preset 

Random Number Seed 
Subroutine, Current «=. 
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a_i nS ee 


eS 
\ {Pon 


RANF Function as 2 7-9 SAVE Steel a 2-25 
RANSET Subroutine hare =9 Scale Factors ha! ne — . 5-8; D-5 
RDRW Subroutine erie B=12 Sealing “ea 33 che D-1; D-6 
READ: ms ek eat Search Character Stfings ne a ue A-6 
ists ee 5~21 Sequential:.~. ee 

Relative Record Subroutine B=12 Access Input/Output ss, gin 5-19 
Statement Meme * 0 caer 5-1; 5-20 - File Characteristics “°° "| a oe 5-19 
Formatted a 5~22 Single-Precision Real subd ner 2-10 
Relative Record _ isis Beis 5-28 Single-Precision: ee ee 
Unformatted - >. a 5. 5-25 Floating-Point: ee Gaia 

Real Format Cenversion. -.0 <9: 5-33 5-5" Source Code '~' oye pei 1-1 
Real Numbers ~~ 34 7. c>e tae9)% 7 oo Source Record Fielas’*"°"" Y le? T1-3 
Memory Representation. rep C 38% 2-10 Source oe wits Ordering | ses 1-4 
Single~Precision rik 2-106 Specification Subp prog ran igh ne igh 6-8 
REAL Statement s813 Vat a 2-2; 2-4 Specified-Time’E écat on. : oe 

Record Separation Indicator = ye 5-17 Subroutine a se 2.5 Bal 
Record Skip:Subroutinie © 9 °'2"~ 7-8 START Subroutine © * ‘#5 Ss B-3 
Records: a 23 eons Woks aes Statement: pa 
Formatted :eaicuei no. a~touses 5-20 ACCEPT ,, 4-17; 5-30 
Uniormatied 5-20 Arithmetic IF as ee 4-5 
RECSKP Subroutine. © © = 17°85 7-8 ASSIGN W8TaLes cic eS aly 4a2 
Recursion o't:.t.+) adi. toute 6~1422 Assigned GO TO Boe Ege 4a? 
RECURSIVEs#12 ¢¢ tea0iust truscar Assignment 3-11; De5 
Statement sab) (a curd. teh 6e42° ~~ Arithmetic WAIeLS VAs S 314 
Subprogram 3:32 (fr us 8a 6<32-. Character aerd.e) 3-16 
Reference, Usthouatical Funauha” re Logical OAR STO 3-11 
Referencing Character Formats: *== 5-18: BACKSPACE =~" GBC ISONOTUE 4. 96 
Relational: BLOCK DATA 6-8; 6-10 
Operations: “ Block IF abO5.409.0G. 9 antag 
Arithmetic 327° BOUNDS: hemes Aso 
Character 3-8 OFF Ses2 and 
Operators 3-1; 4=6-4 ON saci th ae #0 2-18 
Arithmetic i ay Baad CALL afafte tau “448; 4-15; 6-6: 7. 5 
Character 3-8-3 CHARACTER IP NAGEE TE rie 2-5 
Relative Record: sr" 3° GOMMON 28 eli '2= 18; 2-1 
READ Statement 5-28, COMPLEX 1°32 be 9 : 
Subroutine: --< Computed GO TO Lee AN ge 
Read B12. CONTINUE F° Toi gigs 41g 
Write B-42¢ Control "iQiesis > T8Itetelt. 
WRITE Statement 5-29. > -< COPY .5 ee 
Repeat: 3+. ** -EDATA artcitcis ty 
Character String Aqu21 DIMENSION 2 fide 
Subroutine A-21 DISPLAY tas alt 
Repeated: DO- P = 
Field Specification 5-16; 5-16 DOUBLE PRECISION rs a 
Group Specification 5-163 5-16 ELSE — 
RETURN Statement hn 16 IF ~ 
REWIND Statement 5-1; 5-25 END one 


ey. Pt Of S 
“ F ' 


en ¢ & oF ft 


eS «ig 


; 
: u 


to? 


1/Vill 
END FILE 7 5-1;'5426 
END IF 4-73 412 
EQUIVALENCE _ 2-18; 2-23 
EXTERNAL (72:54! 73° afer7 1-4; 2-3: 6-11 
Field a ZF T1-3 
File Positioning =<" -2° 2-05 B25 
FIND eg : BS30 
_ FIXED ; Ae +] 2-2; 2-6 
“RORMAT Cos S97 7 Tang Tee W555. 
Formatted: PUG GL Dita ie 
READ= => af FING Py yp 3 : 5-22 
WRITE “""'S Siosig.-S cise 5-24 
FORTRAN ate TEN 122$ 1. 3 
FUNCTION heal 2-22;-6-3 
Definition eed ee oes a Se 
Definition Statement™ pest aa “2 6a2 
GO TC _ A ae ie ET 41 
IF Pw aree g ~~ Hog 
IMPLICIT meme gtd 2-3; 2-6 
Input/Output arS.o 2a y - =" +501 
INTEGER —— 2-23: 2.9 
INTRINSIC =) ee OOS "faye ae 6-115 7-1 
LOGICAL Seedies shir gtnanes 2-2; 2-5 
Numbers 
OPEN es 
Ordering,-Source = ous 
PAUSE 
READ 
Relative Record j saan! 
REAL Si BSA 2:2: 2-4 
RECURSIVE 1-4; 6-12 
RETURN as ae 4-16 
REWIND © "S860. 2 eit 5-1; 5-25 
SAVE  Qa23 
Statement Function > Befinitign ' 6=2 
STOP 4-17 
SUBROUTINE 2-22 
Type rrwsr are. 17; 6-3; 6-6 
~Gonversion, Assignment - T3-12 
Unconditional GO TO -2 <8 4-4 
Unformatted: “9% © 58” 
READ 5-25 
WRITE 5-25 
WRITE 5-1; 5-20 
Relative Record 5-29 
Status Subroutine; Read - G-3 
STOP Statement — 4-17 
Stored Format Specification 5-18 
String Length, Find Character A-12 
String, Repeat Character A-21 
String: 
Compare Character A-3; A-9 
Convert to Lowercase A-14 
Convert to Uppercase A-17 
Search Character A-6 
Translate A-20 
Subprogram: 
Block Data 6-10 
Dummy Arguments 2-22 


External 

Function 

Internal 

Library 

Recursive 

Specification 

Subroutine 

Subroutine: 

ADATE 

AIRDW 

AISQW 

Analog Input in Random Order 
Analog Input in Sequence 
Analog Output 


Buffered: 

Input 

Output 

BUFIN Si, a. eae 
BUFOUT 

CFILW 

Change File Access Method 
Close Unit Number 
CLOSEW 

Create File 

CRU Output 

Current Random Number Seed 
DATE here 
DATIME 

Delay Continuation. © 
Delayed Bid Task 

Delete File 

DFILW 

Digital Input 


Execute Task 

EXTASK 

Generate: 

File Utility Supervisor Call 
Supervisor Call 

KOMSTR 

Latched Digital Output 
MDATE 

MODAPW 

Momentary Digital Output 
Obtain: 

ASCII Date 

ASCIl Time 

Date 

Date and Time 

Military Date 

Time 

OCRU 
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Open Unit Number 
OPENW 

Preset Random Number 
Generator Seed 
RANCUR 

RANSET 

RDRW 

Read: 

Relative Record 
Record Skip 
RECSKP 


REPEAT “© «sr + 43 . 
Specified-Time: Execution ).:~; = 
START sere 


Statement 
Subprogram 
Subprograms 
SUBSTR 


TIME 

Time-Delayed Execution 
TRANS 

TRNON 


UPCASE  s-visiw. 2retoR 3 = 
WAIT a0 “lh 


Write: 
Relative Record 
WRTRW 


Subscripts< .s¢°:s4 eas 2 4 


SUBSTR Subroutine 
Numbers 


Supervisor Call Subroutine: -. > si 


Generate d22? vit 


Generate File Utility 


SAP ath 
»...:T Field Specification 3-22; < 


SVC Subroutine Sph ayt 
SVCFUT Subroutine ss 


Tabular Output 
TIME Subroutine - 5: -¢jtj2-54 35 
ASCIl ‘ oA 
Date and ~ 
Time-Delayed Execution Subequtine 
TRANS Subroutine rer 
Transferring Character Date ~ 
Translate Character String --; 


/P 
2 


nY 
a iG 
4 AL go oe Lh peg ft 


TRNON Subroutine BACT ICs 
Type: CHttyns 
Conversion: eovintet? 
Assignment Statement: . ‘ou ated 
Input/Output ee 

Statement 2 
Statements Li tthad 
Types, Data an ae 


Bai vi 
., cUnconditional GO TO Statement,: 


Unformatted: sts 
READ Statement ony 
Records stort!” 
WRITE Statement YES 
Unit Number Subroutines ..,: Petre) 
Close SAAS 
»_Qpen “arse 
Unit Numbers... . = silatee 
UPCASE Subroutine aig 
Variable psig 
Variables, Dimensioned -.: >; 
zy .42 
WAIT--Subroutine 4 5*532 
WRITE Lists sirens 
Write: tad ced 


7; tRelative Record Soienutea 
Statementeonsirns  -yers on ~Se 1g 5-293 5-20 


Formatted =~ ~5 tes :.3- 


Relative Record; -~ “2 pats 
' Unformatted ean 
WRTRW Subroutine 


X Field Specification... = > 


oR bF ht Tere 
Z Field. Specification .- 
$74 Ee} 
125 Bas e Ay 
2€3. 182° x = ne 
Pa gl. Sore 
Seay 
nie 
2rire, = 
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